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<7^(Xg)  Variance  of  variable  7.^ 
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T  Shearing  stress 

Tq  Critical  value  of  t:  at  v/hich  g(i:) 

changes  its  characteristics  in  a 
bi-functional  constitutive  rela- 
tionship 

Tj^l  Stress  tensor 

T  Maxiimjim  shear  stress  in  viscometer, 

P6^/2L 

T.^  Shear  stress  at  the  v;all 

T(Tq),  T'<Tj^),  Y(0)   Experimentally  derived  functions 

defined  by  equation  (5«^) 

9  e/2 

|i  Coordinate  of  "latitude"  on  Poincare 

sphere 

X  Extinction  angle:  angle  bett/een  prin- 

cipal optic  axis  and  principal 
direction  of  polarizer  (or  analyzer) 
when  polarizer  and  analyser  are 
crossed.  Of  the  four  angles  thus 
defined,  the  extinction  angle  is 
the  only  one  which  is  less  than  r:/A- 
and  positive. 

>jf  Orientation  angle  which  ni  makes  with 

principal  flow  axis 

yoo  Limiting  value  of  orientation  angle 

as  shear  rate  increases  without 

limit 

'■5"  Effective  orientation  angle 

v-w-  Effective  orientation  angle  at  wall 

'-^  Circular  frequency  of  light 

I,  n,  m  Invariants  of  gj  • 

J.  J 
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C-rad-aate  Council  of  the  University  of  Florida  in 
Partial  Palfill^ent  of  the  Requirenents  for 
the  Degree  of  Doctor  of  Philosophy 

RliEOLOGY  A:-ID  STREA!4ING  BIREFRINGENCE 
OF  AW  AI'IISOTROPIC  FLUID 

3y 

<J.  Eric  Schonblom 

March,  197^ 

Chairinan:  Dr.  E,  Riino  Llndgren 

Major  Departraent:  Engineering  Science,  Mechanics,  and 

Aerospace  Engineering 

The  Intrinsic  viscosity  and  birefringence  of  an 
aqueous  solution  of  Milling  Yello\f  MGoj  a  cciTjrr.ercial 
organic  dye,  are  obtained  expGrimentall3.^,  Each  property 
is  measured  in  a  flow  v/here  the  velocity  is  dependent  upon 
tiro  spatial  coordinates.   It  is  shov/n  that  the  rheological 
and  optical  properties  thus  obtained  nay  be  used  to  cornpars 
hypothetical  velocity  distribu.tions  in  steady  three- 
dimensional  flow's. 

The  rheological  investigation  eriploys  a  Hoppler  Rheo- 
Visconater  in  which  a  ball  slides  without  rotating  through 
the  fluid  within  a  closely  fitted  cylinder.  This  instrument 
has  previously  been  considered  unsuited  for  the  deter- 
nination  of  basic  rheological  constants.  By   modeling  tlie 
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flov/  past  the  ball  on  steady  flow  in  an  eccentric  annulus, 
it  is  shovm  that  the  distribution  of  shear  rates  can  be 
integrated  to  obtain  a  unique  relationship  between  the 
shear  rate  and  the  shear  stress  for  the  fluid.  The 
analysis  is  valid  for  all  fluids  and  can  be  extended 
v/ithout  difficulty  to  viscometers  in  which  the  tightly 
fitted  ball  is  replaced  by  a  cylinder. 

Values  for  the  birefringence  (naxinrom  difference  in 
refractive  index  betvreen  the  principal  optic  axes)  of 
Milling  Yellow  have  been  previously  reported.  The  present 
study  shov/s  that  the  previous  data  exhibit  a  linear  rela- 
tionship between  the  square  of  the  birefringence  and  the 
shear  rate.  An  analysis  demonstrates  that  as  a  result,  in 
a  square  pipe,  the  fringe  order  at  the  wall,  squared,  sho-old 
vary  linearly  with  the  mass  flow  rate  through  the  pipe. 
This  expectation  is  confirned  experimentally,  and  the  bire- 
fringence is  calculated  from  the  data. 

\'lhen   birefringent  fluids  are  observed  in  flov;  oetween 
two  polarizers  a  fringe  pattern  is  seen.   Such  patterns 
have  been  used  to  obtain  pressure  and  velocity  distri- 
butions in  two-diriensioiial  flovrs  and  to  estimate  lift  and 
drag  coofficisnts.  '..ihen  conducted  vrith  Killing  Yellow, 
such  studies  have  been  limited  to  the  extremely  low  flow 
rates  at  vrhich  Hilling  Yellow'  s  birefringence  and  shear 
stress  vary  linearly  with  the  shear  rateo 


The  results  of  the  present  study  extend  these  methods 
to  include  steady  three-dinensional  flows  in  wiiich  velocity 
variations  a-long  the  light  paths  are  permissible.   Further, 
shear  rates  for  ^ihich  the  birefringence  and  shear  stress 
vary  non-linear ly  are  no  longer  excluded.  Although  the 
direct  determination  of  velocity  distributions  from  fringe 
patterns  remains  impractical,  the  pattern  v;hich  corresponds 
to  any  assumed  velocity  distribution  may  be  computed  and 
compared  with  the  fringe  pattern  obtained  experimentally. 
The  method  by  which  fringe  patterns  m.ay  be  calculated  once 
the  velocity  distribution  has  been  assumed  is  outlined 
schematically. 

A  hypothetical  distribution  of  shear  rates  for  Milling 
Yellow  flowing  in  a  rectangular  conduit  has  not  been 
attempted  for  the  rheological  relationship  obtained  with 
the  Hoppler  Rheo-Viscom.eter ;  however,  the  application  of 
other  constitutive  relationships,  notably  that  for  a  pov/er- 
law  fluid,  is  considered  briefly. 
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CHAPTER  OriE 
DTTRODUCTIOII 

The  velocity  distribution  of  anisotropic  liquids 
flowing  steadily  in  rectangular  pipes  can  be  constructed  in 
certain  cases  from  a  knowledge  of  the  optical  and  rheo- 
logical  properties  of  the  fluid.  Specifically,  if  the 
material  is  birefringent,  so  that  the  refractive  index 
bears  a  directional  dependence  upon  the  shear  rate,  then 
any  hypothetical  velocity  distribution  may  be  confirmed 
or  denied  by  observing  the  fringe  pattern  which  results 
when  the  flow  is  observed  between  crossed  polarizers.  The 
successive  steps  in  such  an  evaluation  are  as  follows: 

DeterEination  of  constitutive  relationships.  -  Con- 
stitutive relationsxhips  must  be  provided  which  describe 
the  optical  and  Theological  propertiec  of  the  material. 
Determination  of  ^hgaiL^iaJ:e__distiiiJgutij3n  ~  Based 
upon  the  rheological  properties  of  the  m.edium,  a  compat- 
ible distribution  of  shear  rates  for  steady  flow  in  rectan- 
gular pipes  must  be  calculated.  Depending  upon  the  complex- 
ity of  the  rheological  relationship,  the  m^athematical 
solution  of  this  boundary  value  problem  rn^ay  be  exact  or 
approximate , 


± 


Integration  to  obtain  fringe  patterns.  —  Once  the 
distribution  of  sliear  rates  is  knovm,   integration  of  the 
dependent  optical  properties  along  each  light  path  v;ill 
determine  the  relative  intensity  of  the  emergent  light 
beam.  Fringe  patterns  thus  obtained  may  be  compared  with 
experimental  data  to  evaluate  the  relationships  derived 
in  the  previous  steps. 

Scope  of  Dissertation 

This  dissertation  is  concerned  primarily  with  the 
first  of  the  thiree  steps  just  listed  and  with  the 
properties  of  a  single  birefringent  medium:  an  aqueous 
solution  of  a  coimmercial  organic  dye,  Milling  Yellow  IIGS. 
In  previous  investigations  the  optical  and  rheological 
properties  of  Milling  Yellov;  solutions  (referred  to  here- 
after as  simply  "Milling  Yellow")  have  been  measured  in 
viscoraetric*  flows  using  a  concentric  cylinder  polariscope 
and  a  capillary  viscometer  respectively.  The  present 
study  utilises  non-viscometric  flovrs  to  measure  the  optical 
properties  at  the  wall  of  a  nearly  square  conduit  and  the 
rheological  properties  within  the  eccentric  annulus  of  a 
sliding  ball  viscometer.   Since  neither  of  these  meas- 
urements seem,  to  have  been  employed  previously,  analyses 
are  provided  bo  support  the  present  applications. 


*A  flow  is  viscometric  for  the  purpose  of  this  dissertation 
if  the  velocity  field  has  the  form  ui  =  0,  U2  =  0,  U3  =  u(x) 
vjhere  x  is  a  single  spatial  coordinate.  For  a  more  general 
definition,  see  Coleman,  Markovitz,  and  Noll  (i960). 


The  determination  of  the  distribution  of  shear  rates 
for  Milling  Yellov;  flowing  in  a  rectangular  conduit  has 
not  been  attempted  for  the  measured  rheological  rela- 
tionships; hov/ever,  the  application  of  other  constitutive 
relationships,  notably  that  of  a  so-called  power- law  fluid, 
is  considered  briefly. 

An  optical  analysis  is  performed  to  demonstrate  the 
means  by  which  the  resultant  fringe  pattern  may  be  obtained 
once  the  preceding  steps  have  been  accomplished.   It  is 
shovrn  that  if  the  optical  properties  do  not  change  along 
a  given  light  path  through  the  flowing  medium  the  optical 
relationship  simplifies  to  a  familiar  result  from  two- 
dimensional  optical  stress  analysis. 


CI1A.PTER  TVJO 

RESULTS  OF  PREVIOUS  IF/EST IGAT lOtIS 
AilD  THEIR  ItlPLIC AT  IONS 


This  assessment  of  the  present  state-of-the-art 
is  in  three  parts.  The  first  is  devoted  to  studies  in 
which  birefringent  liquids  have  been  used  to  obtain 
information  concerning  velocity  fields.  Emphasis  is 
laid  upon  those  studies  in  v;hich  Milling  Yellow  was  the 
birefringent  mediiLm,  The  second  part  is  concerned  with 
the  physical  properties  of  Milling  Yellow  and  includes 
a  discussion  of  continuum  mechanics  and  m.odel  construction 
as  they  relate  to  Milling  Yellow' s  rheology.   The  final 
part  describes  previous  investigations  of  velocity 
distributions  in  rectangular  pipes. 

In  Chapter  Five,  preceding  the  analysis  of  the 
Hoppler  Rheo-Viscometer ,  is  a  review  of  the  rolling  ball 
viscom.eter,  the  falling  cylinder  viscometer,  and  the  ball 
and  tube  flow  meter,  subjects  too  specific  for  inclusion 
in  this  more  general  chapter* 


Birefr indent  Flo\-/  Fields 

The  first  reports  of  streaming  birefringence  are 

those  of  Mach  (I873)  and  Maxv/ell  (I873)  a  century  ago. 

Said  MaxT.;ell  (l873j  p.  46): 

I  am  not  av;are  that  this  method  of  rendering 
visible  the  state  of  strain  of  a  viscous  liquid  has 
been  hitherto  employed. 

Although  nany  theories  have  arisen  from  this  humble 
beginning,  the  emplojinent  of  birefringence  for  the  quanti- 
tative investigation  of  flov;  fields  has  remained  scant  to 
the  present  day. 

The  most  popular  media  for  these  studies  have  been 
suspensions  of  colloidal  bentonite  and  solutions  of  organic 
dyes,  notably  Milling  yellow.  Dewey  (1941)  observed  t\TO- 
dimensional  flow  patterns  with  bentonite  and  concluded  that 
quantitative  velocity  gradients  could  be  obtained  from  such 
data.  Similar  studies  by  Weller  (194?)  were  hampered  by 
the  high  viscosity  of  the  polymeric  medium  vrhich  he  used. 
Vinogradov's  (1950)  work  with  colloids  provided  pictures 
of  two-dimensional  flows  around  circular  obstacles  and 
suggested  applications  for  lubrication  theory.  Rosenberg 
(1952) J  another  user  of  bentonite,  described  the  optical 
properties  of  his  medium,  recommended  suitable  concen- 
trations and  colloidal  dimensions,  and  suggested  methods 
of  using  two-dimensional  models  to  calculate  pressure 
distributions,  lift  and  drag  coefficients,  velocity  distri- 
butions and  streamlines.  He  concluded  that  applications 


to  turbulence  would  remain  qualitative  reflecting  Binnie's 
(1945)  experience  with  dilute  solutions  of  bensopurpurin. 
Later  Lindgren  (1953  et  sea.)  and  Way land  (1955)  used 
bentonite  to  visualize  turbulence  but  not  for  the  purpose 
of  computing  the  velocity  field* 

All  of  the  early  quantitative  studies  v^ere  hampered 
either  by  high  viscosities  (as  in  VJeller's  case)  or  by 
marginal  birefringence  (with  bentonite) .  These  consid- 
erations prompted  Jury  in  1950  ^o  suggest  to  Fields  the 
investigation  of  various  organic  dyes  for  their  feasibility 
as  birefringent  media.  Fields  (1952)  concluded  that  the 
most  likely  candidate  for  such  use  was  an  aqueous  solution 
of  commercial  Milling  Yellow.  A  preliminary  study  of  its 
usage  by  Peebles,  Garber,  and  Jiory  (1953)  ratified  this 
conclusion  and  sparked  some  independent  studies  by  other 
investigators.  Although  the  1953  report  did  not  attempt 
a  quantitative  evaluation  of  the  flovf  fields  "v/hich  it 
investigated,  it  did  include  photographs  of  the  two- 
dimensional  flow  patterns,  details  for  the  preparation  of 
the  dye  solution,  determination  of  the  density  (I.005  gm/cis^) 
and  a  plot  of  apparent  viscosity  versus  temperature  for  foui* 
different  dye  concentrations.  The  latter  information  is 
discussed  in  more  detail  in  the  next  section. 

The  first  of  the  independent  studies  using  Milling 
Yellow  \vas  completed  by  Hargrove  and  Thurstone  (1957)  who 
observed  flow  through  an  orifice.  Although  this  study  was 


7 


not  quantitative,  the  usefulness  of  the  medium  prompted 
its  further  use  by  Thurstone  for  the  investigation  of 
v/ave  propagation:  Thurstone  (196l);  Thurstone  and  Scl-irag 
(1962,  196^);  Cerf  and  Thurstone  (1964).  From  these 
studies  emerged  numerical  values  for  the  viscoelastic 
properties  of  Milling  Yellow  which  correlated  well  with 
the  measured  wavelengths  and  propagation  velocities  of 
small  amplitude  waves.  Thurstone  (I96I)  also  replicated 
the  earlier  density  measurement. 

Other  independent  studies  vrere  conducted  by  Swanson, 
ocheuner,  and  Ousterhout  (1965)  and  Swanson  and  Ousterhout 
(1965).   They  assumed  a  linear  relationship  between  bire- 
fringence and  shear  rate  and  demonstrated  the  means  by 
which  a  two-dimensional  flow  field  could  be  calculated  from 
such  data.  They  also  described  a  flow  tunnel  built  for 
this  purpose. 

While  these  independent  studies  were  under  way, 
Peebles  and  his  students,  particulaj-ly  Prados,  continued 
the  original  work  at  the  University  of  Tennessee,  Prados 
(1957)  and  Prados  and  Peebles  (1959)  obtained  velocity 
profiles  for  two-dimensional  flow  in  straight  channels, 
in  converging  and  diverging  channels,  and  in  a  straight 
chan-nel  around  a  cylinder.   Bogue  and  Peebles  (I962) 
suggested  a  technique  for  obtaining  velocity  profiles  from 
isochiTomatic*  fringes  only.  Their  technique  was  applied 


*Isochromatic:  optical  response  (colored  in  white  light) 
which  depends  only  upon  the  birefringence,  An. 
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to  data  obtained  in  a  converging  channel  by  Liu,   Liu  and 
Peebles  (I963)  indicated  that  Hilling  Yellow  can  be  used 
to  describe  two-dimensiona-l  flov;s  in  converging  and 
diverging  charjiels.  in  free  jets,  and  in  wall  jets. 
Building  upon  work  by  Bogue  and  Peebles  (1962)  and  Eirsch 
(1964),  Peebles  and  Liu  (1965)  describe  in  detail  the 
numerical  technique  by  which  two-dimensional  velocity 
profiles  may  be  obtained  from  an  isochromatic  pattern  using 
a  laminar  expanding  jet  as  the  experimental  configuration. 

An  important  restriction  upon  this  research  was  that 
quantitative  evaluation  of  flow  fields  using  Milling  Yeliov/ 
appeared  possible  only  at  extremely-  lox;  flox-/  rates  where 
the  optical  and  riieological  properties  vary  linearly  x-zith 
the  shear  rate  as  predicted  hy   the  theories  described  later 
in  this  chapter.  Further,  only  txvo-drlmensional  config- 
urations, in  which  the  variation  in  fluid  velocity  along  a 
given  light  path  csji  be  neglected,  were  considered 
tractable  for  analysis.   V/hen  these  restrictions  \vere 
observed,  velocity  fields  could  be  calculated  irith  an 
5-''erage  error  of  about  13  percent  according  to  Peebles 
and  Liu  (1965). 

It  should  be  mentioned  that  most  of  the  flow  field 
investigations  were  accompanied  b3''  concurrent  examinations 
of  the  optical  and  rheological  properties  of  the  medium. 
Those  findings  v;hich  rela.te  to  Killing  Yellow  are  reviewed 
in  the  section  v/hicli  follows. 

Excluding  poly!..eric  media  in  which  elastic  properties 


predominate  (extruded  polyethylene  vas  observed  by  Wales, 
1969,  through  windows  set  in  the  long  walls  of  a  capillary 
slit,  for  example)  one  three-dimensional  flow  field  has 
been  examined  quantitatively.  Diirelli  and  Horsard  (1972) 
photographed  flow  around  a  cylinder  in  a  rectangular 
channel  v;ith  an  aspect  ratio  of  C.7^;  that  is,  the  light 
path  was  actually  shorter  than  the  charjiel  width.  This 
arrangement  violated  the  requirement  set  by  the  Tennessee 
studies  for  two-dimensional  flows  based  on  Purday's  (19'^9) 
estimate  that  the  light  path  must  be  5  to  10  times  the 
width  of  the  chan-nel.  Durelli  and  Llorgard  chose  to  treat 
the  flow  as  two-dirf.ensional  and  calculated  average  veloc- 
ities along  each  light  path.  This  assumption  yielded  good 
agreement  with  local  velocities  obtained  by  averaging 
speeds  measured  from,  streak  photographs  of  hydrogen  bubbles 
at  three  locations  in  the  sam.e  channel. 

From  this  review  of  flow  analyses  using  birefringent 
madia  it  is  evident  that  there  is  a  need  for  a  teclinique  by 
wnich  three-dir.ensionjal  flows  can  be  considered.   It  would 
be  helpful  if  the  restriction  to  extrem.ely  low  flow  rates 
could  be  relaxed  or  eliminated c 
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Physical  Properties  of  Milling  Yellov/ 

General  Properties 

The  meaiiira  for  the  present  investigation  was  obtained 
by  dissolving  in  water  a  corarnercial  dye  designated  by  the 
Society  of  Dyers  and  Colour ists  (1971)  as  Colour  Index 
Acid  Yellow  44,  The  coorrion  name  is  Milling  Yellow.  The 
trade  name  for  the  commercial  product  supplied  by  the 
Keystone  Aniline  and  Chemical  Company,  Incorporated, 
Chicago,  is  Milling  Yellow  NGS.   It  is  this  product  which 
was  used  in  the  current  investigations  and  unless  otherwise 
indicated,  the  term  Milling  Yellow  in  this  dissertation 
will  refer  to  solutions  of  this  commercial  product  rather 
than  the  pure  dye stuff, 

Sv;anson  and  Green  (1969)  provide  a  number  of  details 
concerning  the  ph^^sical  chemistry  of  Milling  Yellow.  They 
give  the  structural  formula  as 

CH,  CH 


1 


3  V^^3 


I 


GOH   CHo  CKo    COH 

!^     ) /-. ^.v 


b         NaS03  NaS03        0 
and  state  that  the  birefringence  is  due  to  a  solid  phase 
precipitated  from  solution.  They  describe  this  solid 
phase  as  consisting  of  transparent,  rhombic  crystals  with 
an  aspect  ratio  of  5*7   and  "strong  inherent  polarization." 


11 


Swanson  and  Green  state,  and  the  supplier  confirms,  that 
the  presence  of  impurities,  notably  NaCl  and  HaaSO^  with 
some  sodiiiin  acetate,  may  constitute  more  than  30  percent 
of  the  commercial  product. 

In  dilute  solution  (less  than  1  percent)  Milling 
Yellow  is  lemon  yellow  and  highly  transparent,  Swanson 
and  Green  obtained  birefringence  with  pure  dye  solutions 
having  concentrations  as  low  as  0.1  percent  by  salting 
the  solution  I'Tith  electrolytes,  but  such  solutions  \/ere 
highly  unstable. 

At  higher  concentrations  Milling  Yellow  is  orange 
and  deeply  colored.  The  preparation  of  the  medium  used 
in  this  dissertation  v/as  basically  that  described  by 
Peebles,  Garber,  and  Jury  (1953).  The  dye  v;as  mixed  v/ith 
water  at  a  weight  concentration  below  that  desired  and 
heated  to  just  under  IOC  C,  At  this  temperature  water 
was  evaporated  until  the  desired  concentration  was 
obtained.  Further  details  are  provided  in  Appendix  B, 

Although  earlier  investigators  follow  Peebles,  Garber, 
and  Jury  in  suggesting  dilution  of  the  concentrated  medium 
to  obtain  the  desired  level  of  birefringence  for  a  given 
experiment,  this  dissertation  conc-ors  with  Sv;anson  who 
cautioned  against  dilution  since  the  equilibrium  of  the 
medi^om  is  disturbed  v/hen  distilled  v/ater  is  added,  and, 
on  occasion,  sedimentation  may  result. 

It  was  observed  that  v;hen  Milling  Yellov;  is  suddenly 
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diluted  to  about  1  percent  there  is  a  short  period  during 
which  significant  birefringence  remains  in  the  dispersed 
mixture,  but  the  viscosity  approaches  that  of  water. 
Investigators  who  are  willing  to  tolerate  rapid  changes 
in  the  optical  properties  of  the  medium  may  find  this 
imstable  dispersion  a  useful  medium  for  the  observation 
of  qualitative  phenomena.  Although  the  birefringence 
soon  disappears,  the  color  remains  the  deep  orange  vrhich 
characterizes  solutions  concentrated  by  heating. 

There  is  disagreement  concerning  the  stability  of 
Killing  Yellow  preparations.  Peebles,  Garber,  and  Jury 
(1953)  detected  no  qualitative  differences  in  the  observed 
optical  properties  of  their  1.5  to  1.8  percent  medium 
over  a  period  of  10  months,  nor  was  there  a  perceptible 
darkening  after  more  than  a  week's  continuous  contact  with 
iron  pipe,  steel,  copper,  brass,  or  rubber.  Prados  and 
Peebles  (1959)  did  report  darkening  of  their  1.3  percent 
solution  within  two  weeks  of  preparation  except  for  small 
samples  stored  in  glass  bottles  which  remained  unchanged 
after  three  months*  Peebles,  Prados,  and  Honeycutt  (1965) 
emphasized  that  small  changes  in  concentration  due  to  evap- 
oration have  a  marked  influence  upon  both  the  optical  and 
rheological  properties  of  the  mediiom. 

During  the  current  investigation  it  was  found  that  a 
significant  concentration  gradient  may  develop  between  the 
surface  and  the  bottom  of  the  storage  container  due  to 
evaporation  within  the  container  followed  by  draining  of 
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condensate  from  the  lid.  The  refltocing  action,  unless 
controlled  by  floating  plastic  sheeting  on  the  fluid 
surface,  seems  to  lead  to  sedimentation. 

From  the  preceding  paragraphs  it  is  clear  that  the 
preparation  of  standardized  birefringent  media  having 
specified  properties  is  not  practical  due  to  variability 
in  the  commercial  dyestuff,  apparent  instability,  and 
marked  variation  in  properties  arising  from  evaporation. 
The  recommended  procedure  is  to  measure  all  significant 
properties  at  the  time  of  each  use.  This  has  been  done 
in  the  present  stud;/  and  in  every  previous  quantitative 
investigation  using  Milling  Yellow. 
Optical  Properties 

In  non-steady  shearing  flow  the  optical  properties 
of  Milling  Yellow  have  both  in-phase  and  out-of-phase 
components  which  have  been  studied  by  Thurstone  and  Schrag 
(1962,  1964)  and  Cerf  and  Thurstone  (1964).  They  found 
that  the  birefringence  is  highly  dependent  upon  strain  as 
veil  as  strain  rate  particularly  below  room  temperature 
at  oscillatory  rates  less  than  1  Hz.  Thus,  as  Harris 
,1970)  concludes,  the  analysis  of  unsteady  flows  is  not 
possible  in  the  general  case.   In  the  present  study  only 
steady  state  conditions  are  considered  in  the  measurement 
of  the  optical  properties  and  the  out-of-phase  components 
are  neglected. 

Tv;o  in-phase  optical  characteristics,  the  birefringence 
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and the  extinction  angle,  serve  to  define  streaming  bire- 
fringence in  steady  flow.   Besides  Harris  (1970),  Jerrard 
(1959)  and  Peterlin  (1956)  have  described  these  character- 
istics in  useful  review  articles. 
Theor7v"  of  birefringence 

Based  upon  earlier  work  by  Jeffery  (1922),  Boeder 
(1932)5  Peterlin  and  Stuart  (1939),  and  Snellman  and 
Bjornstahl  (1941),  the  birefringence  of  a  suspension  of 
rigid  non- interacting  ellipsoids  has  been  calculated  by 
Scheraga,  Edsall,  and  Gadd  (1951).  This  body  of  theory, 
which  predicts  a  linear  dependence  of  birefringence  upon 
shear  rate  at  low  flow  rates,  has  been  applied  to  Killing 
Yellow  by  Cerf  and  Thurstone  (1964)  for  the  assessment 
of  small  amplitude  oscillations  and  by  Peebles,  Prados, 
and  Honeycutt  (1965)  and  Swanson  and  Green  (1969)  to 
estimate  particle  size.  The  theory  fails  when  there  are 
interactions  between  particles  or  v/hen  the  particle 
dimensions  exceed  the  upper  limit  of  10"  meters  set  by 
Peterlin  and  St^aart  (1939)  and  Snellman  and  Boornstahl 
(1941) „   Little  is  knovm  about  the  microstructure  of  bire- 
fringent  solutions  of  Milling  Yellow.   Cerf  and  Thurstone 
(1964)  observed  crystals  betx>reen  1  and  2  microns  in  length 
under  the  electron  microscope  but  do  not  report  how  the 
solution  was  prepared  for  viev;ing  in  a  vacuum.*  Recent 


*At  my  request  I^.  Franlc  Hearne  made  a  microscopic  obser- 
vation of  a  2.8  percent  Killing  Yellow  solution  under  an 
oil- immersion  magnification  of  XIOOO.  He  observed  no 
crystals  but  did  obtain  stress  birefringence  in  the  clear 
mediiom  by  pressing  upon  the  cover  glass. 
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descriptions  of  lyotropic  nesophases,  such  as  those  given 
by  Eartshome  and  Stnart  (1970)  suggest  a  viable  alter- 
native to  the  usual  assumption  that  the  properties  of 
Milling  Yellow  are  due  to  a  crystalline  precipitate  of  the 
type  described  by  Swanson  and  Green  (1969).  The  hydro- 
dynamics of  such  mesophases  requires  further  investigation 
before  the  applicability  of  suspension  theory  can  be 
assessed.  i 

Experimental  ireasureaent  of  birefringence 

The  introductory  studies  of  Fields  (1952)  aid  Peebles, 
Garber,  and  Jury  (1953)  provided  qualitative  inforna-tion 
about  the  optical  properties  of  Milling  Yellow,   The  first 
quantitative  evaluation  was  reported  by  Prados  (1957)  and 
Peebles  and  Prados  (1959)  who  calibrated  their  solution 
(roughly  lo3  percent  dye)  in  simple  shearing  flow  usiLng  a 
concentric  cylinder  polariscope.  They  also  verified  that 
the  measureiaent  of  the  distance  between  fringes  ±n  parallel 
channel  flow  could  be  used  as  an  alternative  method  of 
calibration.  The  latter  method  is  still  in  use  (e.g.  Dur- 
elli  and  Norgard,  1972).    Peebles  and  Prados  obtained  a 
nearly  linear  relationship  between  birefringence  and  rate 
of  deformation  for  rates  up  to  19  sec"''.  A  marked  temper- 
ature dependence  was  observed.  The  birefringence  at 
24.75*  C  was  11  percent  higher  than  at  24.95"  C  and  36 
percent  higher  than  at  25.20°  C,  The  extinction  angle 
. was  measured  only  at  25°  C  and  dropped  monotonically  from 
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45°  at  negligible  rates  of  shear  to  about  28°  at  20  sac"''. 

Tile  data  of  Thiorstone  and  Schrag  (1962)  are  of 
limited  usefulness  to  the  present  study  since  they  consid- 
ered oscillatory,  rather  than  steady  shear  flows.  Their 
1.72  percent  solution  showed  a  progressive,  20- fold 
reduction  in  the  optical  coefficient  (loosely,  the  bire- 
fringence) as  the  temperature  was  raised  from  12°  C  to 
42°  C.  V'hen  the  temperature  was  held  constant  at  23°  C, 
the  coefficient  showed  little  change  as  oscillatory 
frequencies  increased  from  10"^  to  1  Hz,  but  dropped 
rapidly  with  further  increases. 

The  range  of  dependence  of  birefringence  upon  shear 
rates  was  extended  by  Hirsch  (1964)  in  his  study  of 
diverging  ducts,  but  the  most  extensive  study  was  reported 
by  Peebles,  Prados,  and  Honeycutt  (1965)  ^■''ho  again  used 
a  concentric  cylinder  polar iscope  for  their  measurements. 
For  shear  rates  ranging  up  to  2500  sec""*  and  concentrations 
between  1.248  and  1.455  percent  by  weight,  they  found 
increasing  non-linearity  as  shear  rates  increased,  although 
"hay  identify  a  possible  "second  range  of  linearity"  at 
the  highest  shear  rates.  All  of  their  data  were  taken  at 
25°  C.  Although  the  concentrations,  which  were  measured 
very  accurately  by  evaporating  samples  to  dryness,  are 
reported  to  four  significant  figures,  dispersion  results 
when  a  correlation  is  attempted  between  concentration  ajid 
the  shear  rate  required  to  obtain  a  given  fringe  order  in 
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the  polariscope.  Figure  1  shows  this  dispersion,  some  of 
which  may  be  due  to  inaccuracies  in  replotting»  The 
remaining  variability  can  be  attributed  to  the  use  of 
the  coEinercial  dyestuff  and  to  the  difficulty  of  preparing 
a  standardized  medium  as  alluded  to  earlier.  In  any  case 
it  is  clear  that  birefringence  increases  markedly  with 
dye  concentration. 

Peebles,  Prados,  and  Honeycutt  also  measured 
extinction  angles  over  the  same  range  of  concentrations. 
They  found  that  the  more  concentrated  solutions  e:diibit 
asymptotic  values  for  the  extinction  angle  in  the  vicinity 
of  20°  as  the  shear  rate  increases  above  40  sec"""  .   At 
lower  concentrations  a  similar  asymptote  is  reached,  but 
at  higher  shear  rates. 

No  expressions  for  either  the  birefringence  or  the 
extinction  angle  are  advanced  by  the  authors  to  represent 
their  findings.  Consequently,  for  the  purpose  of  the 
present  study  it  has  been  necessarj'"  to  construct  empirical 
relationships  which  describe  the  data  of  Peebles,  Prados, 
and  Honeycutt,  This  has  been  done  in  Chapter  Four. 

Swanson  and  Green  (19^9)  were  concerned  only  v;ith 
the  minimum  concentration  at  which  birefringence  could  be 
observed  and  not  with  its  magnitude.  They  hypothesize  that 
the  variability  of  Mil].ing  Yellow  preparations  is  due  to 
a  small  fraction  of  the  dissolved  material,  as  little  as 
0,04  percent,  which  exists  in  suspension.   Their  hypothesis 
is  not  evaluated  in  this  dissertation. 
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FIGu'RE  1.   Dispersion  in  concentration  data  of  Peebles, 
Prados,  and  Honey cutt  (1965) 
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Rheologlcal  Properties 

The  Theological  behavior  of  nedia  such  as  Milling 
Yellow  may  be  described  rigorously  in  terms  of  continuuia 
mechanics,  theoretically,  but  v;ith  less  rigor  in  terms 
of  hydrodynamic  models,  or  empirically,  based  upon  e:rper- 
imental  evidence.  Each  of  these  methods  is  discussed 
in  turn. 
Continuum  mechanics 

V/ith  the  advent  of  liquid  crystals  as  a  practical 

media  for  electronic  display  devices  (Caulfield  and  Soref , 

1971)  is  one  report  among  many),  there  has  been  a  great 

increase  in  publications  relating  to  the  constitutive 

behavior  of  anisotropic  materials.  Hot  all  such  reports 

have  been  useful.  One  reviewer,  Kisiel  (1968,  p.  104-3) 

spoke  for  many  when  commenting  upon  a  study  which  shall 

remain  nameless: 

This  investigation  belongs  to  a  class,  abundant 
at  present,  of  papers  dealing  with  very  general  problems 
with  limited  applicability  to  the  solving  of  practical 
qiaestions. 

An  overview  cf  the  current  state  of  the  art  indicates 
that  rigorous  application  of  the  continuum  mechanics  of 
anisotropic  media  is  limited  to  visconetric  flows  of  the 
simple  type  in  which  the  velocit^r  com.ponents  are  given  by 

Ui   =  0;   U2  =  0;   U3  =  u(x) 
where  x  is  a  single  spatial  coordinate.  Further,  numer- 
ical solutions  are  possible  only  for  the  very  small  class 


20 


Ox  substances,  notably  p-azoxyanisole,  for  which  some,  at 
least,  of  the  necessary  constitutive  constants  have  been 
measured  and  published.  Neither  of  these  conditions  is 
satisfied  in  the  present  dissertation;  hence,  the  dis- 
cussion of  anisotropic  continuum  mechanics  x>7hich  follows 
is  succinct  and  selective. 

Oldroyd  (1950)  established  the  general  procedure  by 
which  constitutive  equations  must  be  constructed  if  the 
necessary  conditions  for  tensor  invariance  were  to  be 
preserved.  Noll  (1958)  introduced  the  concept  of  a 
"simple  fluid":  one  in  v;hich  the  properties  are  completely 
defined  by  the  temperature  and  the  strain  history.  The 
viscometry  of  simple,  non- Newtonian*  fluids  was  examined 
by  Coleman,  Markovitz,  and  Noll  (I966)  in  a  general 
treatise  which  includes  a  bibliography  of  over  350  refer- 
ences spanning  the  period  from  I687  to  1965. 

Specific  constitutive  relationships  for  anisotropic 
fluids  were  formulated  by  Ericksen  (1960a  et  sea.)  and 
Leslie  (1966  et  seq,)  who  postulate  that  at  each  point  in 
the  continuum  there  is  a  preferred  direction  characterized 
by  a  unit  vector,  or  "director,"  dJ.  On  the  basis  of  this 
hypothesis  it  was  found  that,  in  general,  the  constitutive 
stress  tensor  is  non- symmetric,  and  seven  or  more  consti- 
tutive constants  are  required.  The  theory  has  been  applied 
with  some  success  by  Atkin  and  Leslie  (1970)  and  Tseng, 


*Non-Newtonian:  a  substance  is  Newtonian  if  and  only  if  the 
shear  stress  is  directly  proportional  to  the  shear  rate. 
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Silver,  and  Finlayson  (1972)  to  certain  specialized  flows, 

A  more  general  and  even  less  tractable  theory  has 
been  developed  by  Eringen  (1964  et  seg^, )  who  postulates 
a  micromotion  of  the  material  points  which  define  the 
continuvuHo  Associated  with  this  micromotion  are  corres- 
ponding m.icromoments  and  microinertia.  A  linear  consti- 
tutive relationship  is  obtained  at  the  price  of  additional 
unlcnoxm  constitutive  constants.  The  current  literature  is 
replete  with  argument  concerning  the  existence  of  the 
various  constants,  with  their  signs,  and  with  the  rela- 
tionships, frequently  in  the  form  of  inequalities,  among 
tiiem,  Truesdell  (1965)  has  pointed  out  that  the  complexity 
of  modern  continuum  mechanics  is  a  reflection  of  nature 
and  requires  no  apology,  but  a  hopeful  reading  of  the  most 
recent  review  of  anisotropic  continuum  mechanics  by  Ariman, 
Turk,  and  Sylvester  (1973)  leads  only  to  the  conclusion 
that  the  theory  is  not  yet  useful. 
Model  construction 

As  an  alternative  to  the  utilization  of  rigorous,  but 
complex  constitutive  relationships  for  a  continuum,  many 
authors  have  elected  to  model  anisotropic  behavior  in  terms 
of  the  effect  which  the  presence  of  microscopic  particles 
in  a  Nevrtonian  medivun  ha-S  upon  the  macroscopic  properties 
of  the  m.ixture.  The  success  of  such  theories,  of  which 
Einstein's  (1906)  calculation  of  the  intrinsic  viscosity 
of  a  suspension  of  rigid  spheres  is  the  classic  example, 
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has  led  not  only  to  the  analysis  of  particles  whose  shape 
is  less  well  defined,  as  in  colloidal  suspensions,  but 
also  to  inferences  about  the  microstructure  when  an  ill- 
defined  or  poorly  imderstood  medium  is  found  to  obey  the 
predictions  of  a  particular  theory. 

The  most  influential  body  of  analysis  has  grown  from 
Jeffery's  (1922)  solution  for  the  periodic  motion  of 
rigid  ellipsoids  suspended  in  a  viscous  fluid  undergoing 
uniform  shearing  motion.  Jeffery's  solution  was  open- 
ended,  consisting  of  an  infinite  seb  of  permissible  orbits, 
Other  authors,  notably  Peter lin  (1538),  calculated  the 
distribution  of  orbits  which  would  result  from  pertur- 
bations of  the  particles  due  to  Brovmian  notion  as 
expressed  by  the  rotational  diffusivity  constant.   Inte- 
gration of  such  distributions  leads  to  an  estimate  of  the 
viscosity,  Kuhn  and  Kuhn  (19'^5)?  Scheraga  (1955)  j  sua. 
Leal  and  Hinch  (1971)  are  among  those  v;ho  have  performed 
this  integration. 

Cylindrical  particles  have  been  treated  by  Boeder 
(1932),  who  replaced  the  cylinders  by  ellipsoids  of  high 
axial  ratio,  Burgers  (1938),  who  obtained  the  torques  due 
to  shears  for  true  cylinders,  and  Broersma  (I960),  who 
included  end  affects.   Still  later  Bretherton  (1962) 
demonstrated  that  any  rigid  particle  having  an  axis  of 
revolution  can  be  replaced  by  an  ellipsoid  of  appropriate 
dimensions  and  incorporated  into  the  general  theory. 
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A  common  assumption  of  these  theories  is  that  there 
is  no  interaction  between  the  particles.  V/hen  interaction 
is  permitted,  as  in  Ziegel  (1970)  or  Batchelor  (1971),  the 
analysis  is  greatly  complicated. 

Although  rigid  spheres,  ellipsoids,  and  rods  have 
served  as  the  primary  models  for  the  analysis  of  non-linear 
rheological  behavior,  other  shapes  also  play  an  important 
role.  A  sampling  of  investigations  which  have  served  as 
alternate  models  might  include  the  work  of  Taylor  (193^) 
on  drops,  Debye  (194-6)  on  swarms  and  porous  spheres,  Kuhn 
and  Kuhn  (194-3)  and  Kiriarood  and  Riseman  (1948,  194-9)  on 
chains  and  necklaces,  Simha  (1950)  on  dumbbells,  and 
Frohlich  and  Sack  (194-6)  on  elastic  spheres. 

The  practical  value  of  these  theories  is  that  they 
permit  the  replacement  of  a  complex  constitutive  rela- 
tionship v;ith  many  unknovm  constants  by  a  relatively  simple 
constitutive  equation;  however,  the  coefficients  of  this 
•  equation  will  e:«ihibit  an  involved  (though  theoretically 
explicit)  dependence  upon  the  various  material  parameters, 
and  these  parameters  may  prove  as  difficult  to  measure  as 
zr,e   constitutive  constants  vmich  they  replace.  An  examiple 
is  the  rigorous,  three-constant  constitutive  equation 

'^ij  =  -P^ij  -  ^^ij  -  •i^6n+P)gikgkj  +  ^e^  -r-i 

where  the  constants  are  the  intrinsic  viscosity  /r,  and  the 
primary  and  secondary  normal  stress  functions  G^^and  p. 
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For  the  model  of  rigid  ellipsoidal  particles  in  suspension, 
re  has   been  calculated  by  Saito  (1951)  and  Schoraga  (1955), 
and  the  stress  functions  have  been  obtained  by  Giesekus 
(1962) »  The   study  by  Scheraga  tabulates  its  results  in 
terns  of  the  rotational  diffusivity  constant  and  the  ratio 
of  the  lengths  of  the  major  and  minor  axes  of  the 
ellipsoid.  In  practice  it  has  been  coamoner  to  infer 
these  properties  from  the  macroscopic  properties  rather 
txhan  the  reverse.  Thus  the  model,  even  when  it  is  valid, 
may  not  be  predictive. 
Experimental  measurements  of  rheclogj'" 

Amenemhet's  (154-0  B.C.)  boastful  account  of  his  water- 
clock,  which  was  capable  of  compensating  for  seasonal  vari- 
ations in  viscosit}''  (due  to  temperature  changes),  begins 
the  \iTltteTi   record  of  vhaology.     It  is  clear  that  Amenemhet 
did  his  work  without  the  benefit  of  continuum  mechanics. 

The  present  imowledge  of  the  rheology  of  Milling 
Yellow  is  also  founded  upon  experiment.  Although  the 
empirical  relationships  which  describe  these  data  may 
violate  conditions  of  invariance  prescribed  by  continuum 
mechanics  and  include  constants  which  cannot  be  obtained 
from  model  construction,  they  may  be  employed  with  care 
provided  that  the  flows  to  v/hich  they  are  applied  do  not 
differ  too  greatly  from  those  in  v;hich  the  experimental 
data  v;ere  obtained,  A  more  lengthy  discussion  of  the 
feasibility  of  employing  empirical  relationships  to 
describe  the  rheology  of  Milling  Yellov/  vrlll  be  found  at 
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the  begiimrng  of  Chapter  Five. 

VJhen  the  feasibilitj'-  of  Milling  Yellov;  as  a  bire- 
fringent  nediiim  was  established  by  Peebles,  Garber,  and 
Jury  (1953) >  measurenents  of  the  apparent  viscosity  were 
made  in  a  rolling  ball  viscometer  at  various  temperaturea 
For  solutions  varying  in  concentration  from  1.46  to  2.02 
percent,  a  sharp  exponential  rise  in  viscosity  was 
observed  as  the  temperature  decreased.   In  the  2  percent 
solution,  the  viscosity  doubled  as  the  result  of  a  two- 
degree  temperature  drop.  Above  a  certain  critical  temper- 
ature the  optical  activity  of  the  solution  ceased  and  the 
viscosity  approached  that  of  water*   It  was  recognized 
that  the  apparent  viscosity  had  a  shear  rate  dependence 
v/hich  was  not  obtained  from  the  measurements. 

Prados  (1957)  and  Peebles  and  Prados  (1959)  measured 

the  viscosity  of  a  1.3  percent  solution  but  did  not  report 

the  results.  For  shear  stresses  less  than  5  d:^Tie/cm^, 

Prados  assumed  that  the  viscosity  was  constant,  citing 

Honeycutt  and  Peebles  (1955)  as  his  author itj"-  that  Milling 

Zellow  solutions: 

...exhibit  marked  non-Newtonian  behavior  \-Then 
subjected  to  shearing  stresses  greater  than  five  to  ten 
dynQS   per  square  centimeter.   (Prados,  1957?  pp.  55-56) 

Thur stone  (I96I)  measured  the  acoustic  impedance  of 

a  1.39  percent  solution  in  a  circular  tube  whose  base  vms 

excited  by  lo\-;- amplitude,  axial,  oscillatory  vibrations. 

At  an  uncontrolled  temperature  between  22°  and  26°  C,  he 
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fo-ond  that  Milling  Yellow  exhibited  viscoelastic  properties. 
That  is,  the  local  stresses  were  a  function  of  both  the 
shear  and  the  shear  rate.  As  the  freqiiency  increased 
frod  10  to  300  Hz,  the  viscous  term  of  the  complex  viscos- 
ity coefficient  dropped  from  66  to  13  centipoise,  the 
elastic  ccnponent  having  about  the  same  magnitude  as  the 
viscous  component  over  this  range • 

Thur stone  and  Sclirag  (1964)  and  Cerf  and  Thur stone 

(1964)  did  not  report  the  viscous  and  elastic  terms 
separately.  Thurstone  and  Sciirag  found  that  the  complex 
viscosity  coefficient  is  approximately  the  same  for  both 
axial  and  transverse  oscillations  of  the  medium.  Cerf 
and  Thurstone  found  that  elastic  forces  predominate  at 
frequencies  below  0.3  Hz,  but  that  viscous  forces  are 
dominant  above  10  Hz  for  oscillatory  shear  -v/aves«  At 
very  high  frequencies  a  limiting  viscosit^r  of  45  centi- 
poise was  obtained  for  their  1.73  percent  solution  at  25"°  C. 

In  stead;/  flows  the  most  extensive  examination  of 
Milling  Yellow  rheology  is  Peebles,  Prados,  and  Honeycutt 

(1965)  who  measured  apparent  viscosities  with  a  capillary 
viscometer  at  25°  C  over  a  concentration  range  of  1.25 

to  1.50  percent.  For  each  sample  thej   obtained  smooth 
mono tonic  curves  for  calculated  values  of  wall  shear 
stress  versus  shear  rate  with  linear  responses  when  the 
shear  rate  exceeded  about  2500  sec"^  «  Ifaezi   the  wall  shear 
stress  corresponding  to  a  given  shear  rate  is  replotted 


27 


versus  concentration,  there  is  a  significant  scatter  of 
the  data,  Just  as  a  similar  replotting  of  the  optical 
data  (Figure  1)  also  resulted  in  dispersion.  This 
confirms  a  difficulty  experienced  in  all  investigations 
including  the  current  one:  specification  of  the  commer- 
clal  dye  concentration  is  insufficient  to  define  the 
properties  of  the  medium  even  at  a  fixed  temperature. 

An  important  result  of  tiie  investigation  of  Peebles, 
Prados,  and  Honey cutt  was  the  demonstration  that  plots  of 
apparent  viscosity  versus  wall  shear  stress  are  independent 
of  the  diameter  (and  L^/D  ratio)  of  the  capillary  in  which 
the  measurement  is  made.  As  Skelland  (1967,  pp.  32-39), 
among  others,  has  pointed  out,  this  coincidence  of  curves 
indicates  the  absence  of  inlet  effects  of  the  type 
described  by  Maude  and  Vihitmore  (1956)  or  of  wall  effects 
such  as  slippage  or  the  radial  migration  of  microscopic 
elements  as  measured  by  Goldsmith  and  Mason  (1961,  1962, 
I964-)  and  Gauthier,  Goldsmith,  and  Mason  (1971).   In  the 
absence  of  a  reliable  explanation  of  Milling  Yellow's 
exceptional  properties,  the  elimination  of  such  effects 
from  consideration  is  welcome, 

Peebles,  Prados,  and  Honeycutt  (1965)  conclude  that 
Hilling  Yellow  is  well  represented,  though  not  tmiquely, 
by  the  Pov/ell-Eyring  equation: 

P-     =     V-QO    +  -^(Po  -  >ioo)  sinh-'(g/Ap) 
where  y  is  the  viscosity  at  shear  rate  g,  and  ]1q,  }1qo  ,   and 
Ap  are  constants.  The  authors  provide  straight-line  plots 
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of  these  constants  versus  concentration,  and  for  the  four 
concentrations  plotted  the  agreement  is  excellent.  Based 
upon  these  plots,  the  following  numerical  relations  can 
be  obtained: 

logio  Po   =   10.72  q  -   12.92, 

logio  poo  =   1.33  q  -   1.63, 

logic  A^  =  -10.37  q  +  16.16. 
The  units  of  p,  7105  and  jioo    are  centipoise,  shear  rates 
g  and  A  are  in  sec"*,  and  q  is  the  weight  percent  of 
Milling  Yellow  in  the  solution. 

The  data  of  Peebles,  Prados,  and  Honeycutt  are  not 
compelling  with  regard  to  the  prediction  of  the  Powell- 
Eyring  equation  that  the  viscosity  vrill  approach  a 
constant  value  at  low  shear  rates.  Since  most  of  the 
quantitative  studies  reported  in  the  literature  were 
conducted  at  very  low  shear  rates,  the  absence  of 
conclusive  data  in  this  range  is  of  major  concern, 
Fortunately,  Peebles  and  Liu  (1965)  have  provided  a  plot 
from  Hirsch's  (1964)  dissertation  which  indicates  clearly 
that  the  viscosity  does  approach  a  constant  value  for 
shear  rates  below  about  5  sec""* .  These  data  were  obtained 
in  a  capillary  viscometer  and  replicate  closely  data  from 
a  '"Rotovisco"  instrument  when  the  two  methods  are  compared 
at  shear  rates  around  50  sec"''.  The  latter  instrument 
shows  the  upper  range  of  linear  responses.   It  should  be 
recalled  that  Prados  (1957)  m.easured,  though  he  did  not 
report,  constant  viscosities  at  loxj'er  shear  rates. 
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Velocity  DlstrlbTition  in  Rectangular  Conduits 
V/henever  a  differentiable  expression  for  the  velocity 
distribution  is  inioi-m,  the  shear-rate  distribution  is 
defined  by  direct  differentiation.   Once  the  shear-rate 
distribution  is  known,  the  birefringence  and  orientation 
angle  (or  extinction  angle)  can  be  calculated.  In  the 
present  dissertation  the  shear-rate  distribution  is 
required  in  a  rectangular  conduit.  This  distribution 
has  not  been  calculated  for  a  fluid  with  Milling  Yellow's 
rheological  properties.  The  review  which  follows  includes 
those  studies  which  show  a  potential  usefulness  in  the 
development  of  such  a  distribution. 
Ilei^onian  Flow  in  Cylinders 

The  determination  of  velocity  distributions  in  pipes 
dates  from  the  experimental  studies  of  Hagen  (1839)  and 
Poiseuille  (1840).   In  modern  derivations  the  equation 
which  bears  their  names 

\vhere  u(r)  is  the  speed  at  a  distance  r  from  the  center- 
line,  R  is  the  pipe  radius,  jUjj  is  the  viscosity,  and  dp/dZ 
is  the  pressure  gradient,  is  obtained  directly  from  the 
Navier-Stokes  equation  for  incompressible  fluids, 

-vp  +  }i-^-^^   +  pF  =  p  du/dt, 
by  recognizing  that 

u^  =  0,  ug  =  0,  u^  =  u(r), 
and  integrating.  The  density  p  is  implicit  in  the  dp/dz 
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term  which,  in  gravitational  fields,  is  sinply 

dp/dz  =  -pg^ 
where  g  is  the  gravitational  constant. 

For  a  pipe  which  is  not  circular  in  cross-section, 
the  ass'junption 

Ui   =   0,   U2   =   0,   U3   =  •u(xi,X2),        (2,1) 
yields 

Pjjv2u  =  dp/dz.  (2.2) 

Exact  solutions  for  this  equation  have  "been  obtained 
for  cross-sections  in  the  shape  of  concentric  circles, 
ellipses,  and  equilateral  triangles.   Lamb  (194-5),  in 
reviev;ing  these  solutions,  points  out  that  the  analysis 
of  laminar  flow  in  a  cylindrical  conduit  is  identical  in 
mathematical  form  to  the  analysis  of  torsion  in  a  uniform 
cylindrical  bar  and  of  fluid  motion  in  a  rotating  cylin- 
drical case,  the  cylinders  in  each  case  having  the  same 
cross-saction,  Tiedt  (1969)  adds  the  reminder  that  the 
analog  is  valid  without  modification  only  if  the  boundary 
of  the  cylinder  is  simply  connected. 

Davies  and  vJhite  (1928)  obtained  the  relationship 
between  the  pressure  gradient  dp/dz  along  a  rectangular 
duct  and  the  Reynolds  number 

Re  =  ^ f.      t   ^ 

where  61  and  62  are  the  half-width  and  half-depth  of  the 
duct. 
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In  the  same  year  Cornish  (1923)  published  the  solution 
to  equation  (2.2)  in  a  rectangular  conduit  in  the  form  of 
a  Fourier  series: 

u  =  --r^^  ((5i2  .  y:2)    +  (2.3) 

32o-i^^(-l)"^    ^Q3  rSTTx1c05h(STr//2oO    \ 
t:3     s"^..     s3  [2oi  Jcosh(sn62/26i)J 

Cornish  successfully  related  the  corresponding  volumetric 
flow  rate  to  e>:pariinental  pressure  gradients.  Further 
data  v/hich  support  this  relationship  are  those  of  Nikuradsa 
(1930)  and  Lea  and  Tadros  (1931),  v/hich  also  shox^  the 
predicted  dependence  of  pressure  gradient  upon  average 
flow  rate.  The  acc-oracy  of  the  velocity  profile  must  be 
inferred  from  Eckert  and  Irvine  (195°)  who  obtained 
excellent  agreement  between  local  velocitj^  measurements 
and  the  Fourier  series  solution  to  equation  (2*2)  for 
triangular  cross- sections.  The  Cornish  solution  is  clearly 
irxaccurate  \jhen   violations-  of  equation  (2,1)  occur  due  to 
secondary  flows  arising  from  convective  effects.  Such 
flows  are  commcr.  even  in  circular  pipes  as  demonstrated 
most  recently  oy   Joimson  (197^) » 
::_on--J'\ewtonian  Flv.lds  in  Rectansjjlar  Conduits 

Of  tlie  various  forms  of  the  Theological  equation 
which  have  been  or  will  ba  suggested  for  Milling  Yellow 
in  this  dissertation,  only  one  has  been  investigated  in 
pipes.  Christiansen,  Ryan,  and  Stevens  (1955)  related 
pressure  gradients  to  average  flov;  rates  for  a  Powell- 
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Eyring  fluid,  but  their  analysis  was  limited  to  circular 
pipes. 

It  will  be  shovm  in  Chapter  Five  that  at  low  flow 
rates  two  of  the  empirical  expressions  for  Milling  Yellow 
reduce  to  the  form 

T  =  K  g^^O 
which  is  the  one- dimensional  form  of  a  power- law  fluid. 
Povfer-lav;  fluids  are  defined  by 

v;here  t^.  ^  is  the  stress  tensor,  gj_j  is  the  rate-of- 
deformation  tensor,  H  is  the  second  invariant  of  g^j? 
and  1{_  and  m  are  constants.  For  the  limited  range  over 
v;hich  the  power  lax/  applies  to  Milling  Yellow,  m  =  1/3 • 
Pox'/er-law  fluids  have  been  investigated  in  rectangular 
pipes  by  several  authors. 

Schechter  (196I)  used  variational  methods  to  obtain 
the  pressure  drop  along  rectangular  pipes  for  power- law 
fluids  having  m  =  0.5)  0.75?  and  1.0  when  the  pipe  aspect 
ratio  &1/62  was  0.25,  0.5,  0.75?  and  1.  He  tabulates  the 
coefficients  to  be  used  in  the  series  solution 

u  =  uSAg  sin(NsTDc/2&i )  sinClTs' TTy/2&2)     (2.4) 
to  obtain  local  velocity  values. 

^"flieeler  and  Wissler  (1965)  elected  to  solve  the  same 
problem  by  finite  difference  methods.  Taking  the  solution 
for  a  Newtonian  fluid  (m  =  1)  as  a  starting  point, 
successive  approximations  for  the  velocity  distribution 
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were  obtained  for  values  of  m  between  0,4  and   1.  Both  the 
stability  of  the  solution  and  the  rate  of  convergence 
decreased  with  m.  Below  m  =  0,5  stability  was  a  serious 
problem,  and  several  hundred  iterations  were  required  for 
convergence  at  the  lowest  value  of  m.  For  square  pipes 
it  was  found  that: 

-dp/ds  =  7.^942  (1.7330  m""'  +  (2.5) 

5.8606f  K  if'/26i+m  . 

For  0.4- <  a  <  1.0,  the  constants  in  this  relationship  were 
accurate  to  four  significant  figures.  Velocity  profiles 
obtained  by  this  method  were  not  published,  but  Wheeler 
and  Wisslsr  state  that  the  profiles  obtained  hy   their 
method  could  be  differentiated  numerically  several  times. 
In  contrast,  differentiation  of  Schechter's  profiles  led 
to  erratic  results.  The  empirical  expression  given  above 
was  verified  experimentally  using  power-law  constants 
obtained  for  their  medium  (sodium  carboxymethylcellulose 
solutions  of  various  concentrations)  by  averaging  the 
measurements  made  in  a  circulaT  pipe  and  a  Couette  visco- 
meter. For  Rejmolds  numbers  less  than  2000  there  was 
ei'cellent  agreement  between  the  predicted  pressuTe  drop 
and  the  corresponding  Reynolds  number. 

Aral  and  Toyoda  (1968)  considered  short  rectangular 
conduits  with  pov;er-law  fluids  having  values  of  m  from 
0,3  to  1.  They  provide  average  wall  shear  rates  in  terms 
of  an  effective  radius: 

^eauiv.  "  26i 62/(5i+&2). 
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The  velocity  distribution  obtained  for  m  =  0,4,  61/62  ~   2, 
is  also  provided  together  with  the  corresponding  shear- 
rate  distribution. 

A  return  to  variational  methods  v;as  provided  by 
Rothene3'-er  (1970)  who  obtained  a  series  of  non-linear 
equations  for  the  coefficients  bg^  of  the  poljoionial 


E   -b 


,2S  „2t 


^   =  JlT^^st^^^  y  "^  (2.6) 

S  =  0  x  =  o 

by  substitution  into  the  non-linear  partial  differential 
equation  governing  pov;er-law  substances  in  pipe  flow: 

dz     21k/  l4i  \[{dxl  \dyl   ^  L^^-^  "^  ^y^J 

Llbxi  ^  i5y/  J   l\5x/  5x2 

+  2  ~  ~  b^u    /5u'  2  d ^u"l  1 
ox  by   oxdy  "^  \djj      by^jj   ' 

where  m'  =  (l-n)/2m.  The  substitution  was  made  at  each 

point  (xjy)  of  a  lattice  distributed  across  one  quadrant 

of  the  cross- section.  This  set  of  equations  was  linearized 

by  substituting  into  the  non-linear  terms  the  values  of 

bg^  obtained  in  the  previous  iteration.  For  the  first 

iteration  the  ITavier-Stokes  solution  obtained  by  Cornish 

(1928)  was  used.  The  boundary  condition  was  met  by  setting 

u  =  0  in  equation  (2.6)  for  lattice  points  along  the  wall 

and  adding  the  resultant  set  of  linear  equations  to  the 

linearized  set  obtained  by  substitution.  The  decision 

not  to  write  equation  (2.6)  in  a  form  which  satisfied  the 

boundary  conditions,  as  was  done  by  Schechter  (I96I),  was 
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dictated  by  Rotheneyer ' s  intent  to  provide  a  nothod  which 
was  appropriate  for  cylinders  of  arbitrary  cross-section. 
In  general  the  number  of  lattice  points  (",y)  v.'as  greater 
than  the  number  of  bg-^  so  that  the  system  of  linearized 
equations  was  overdetermined.  The  e::tra  d:jgrees  of  freedom 
were  used  to  minimize  the  error  due  to  the  b^^  estimate,  a 
least  squares  fit  being  employed.  The  iterative  process 
ended  when  the  computed  flow  rate  through  the  cross-section 
differed  by  less  than  1  percent  from  the  previous  iteration. 

Substances  other  than  Hevrtonian  fluids  and  power- law 
substances  have  received  little  attention  in  flows  through 
rectangular  conduits.   Sokolovskii  (1966)  considered  a 
dilatant  material  with  the  response: 

For  rectangular  pipes  he  considered  orJ.y  freely  dilantant 
movement  (t^  ->  0)  for  which  the  lines  of  constant  velocity 
form  a  set  of  rectangles,  one  inside  the  next. 

Greenberg,  Born,  and  VJetherell  (I96O)  solved  by  finite 
difference  methods  the  torsion  probleri  for  a  square  cylinder 
composed  of  a  material  obeying  the  Ram/oerg- Osgood  stress- 
strain  lav;.  The  fluid  analog  of  a  Ramberg-Osgcod  solid  is 
a  DeHaven  fluid,  defined  by  the  relationship: 

^     =     y^o   g/(l  +  Izz") . 
The  values  of  n  for  which  DeHaven  (1959a,  19?9b)  employed 
this  relationship  were  much  smaller  than  the  values 
preferred  by  Greenberg,  Dorn,  and  '.Tetherell  (i960). 
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Hanzawa  and  Ishikawa  (1970)  investigated  the  problem  of 
Greenberg  et  al.  after  greatly  simplifying  the  boundary 
conditions  by  replacing  the  straight  walls  by  concave 
siirfaces.  Comparable  results  were  obtained  where  the 
stresses  were  highest. 

Litvinov  (1968)  used  variation  methods  upon  empir- 
ical rheological  data  for  polypropylene.  After  expressing 
T  and  g  in  the  form 

Jdi  g  +  d2  g^  +  d3  g3       g  <g^ 
W  g  +  d5  g  ^  gg 

where  T(g)  was  determined  in  one- dimensional  flow,  he 
assumed  that 

and  minimized  the  error  introduced  by  the  coefficients 
over  the  cross- sect ion  of  his  flow. 

In  a  follow-up  of  their  1965  study,  Wheeler  and 
Wissler  (1966)  measured  the  velocity  distribution  of  a 
0.9  percent  solution  of  sodium  carboxymethylcellulose 
flowing  in  a  square  pipe.  By  observing  the  movement  of 
suspended  particles  at  12  locations  in  and  surrounding 
one  quadrant  of  the  cross- section,  they  found  deviations 
of  up  to  7  percent  from  the  velocity  profile  obtained  by 
Wheeler  and  Wissler  (1965).  The  direction  of  the  vari- 
ations was  consistent  with  the  hypothesis  that  there  was 
secondary  flow  within  the  cross-section.  To  test  this 
hypothesis,  the  authors  chose  to  model  the  liquid  as  a 
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Stokasian  fluid*  with  a  constant,  but  non-zero,  normal 
stress  function.  The  velocity  distributions  thus 
obtained  gave  qualitative  support  to  the  hypothesis 
that  secondary  flov/s  v/ere  present.  The  method  by  which 
the  distribution  was  calculated  is  stated  in  general 
terms  and  the  constants  which  were  obtained  for  the 
Stokes  equation  were  not  published. 


*Stokesian  fluids  are  discussed  at  considerable  length 
in  Chapter  Five. 


CHAPTER  THREE 
OPTICAL  ANALYSIS 

As  stated  in  the  Introduction,  when  a  birefringent 
liquid  flows  b8tv;een  two  polarizers,  a  pattern  of  fringes 
is  observed.  The  purpose  of  this  chapter  is  to  derive  a 
method  for  determining  the  amount  of  light  which  emerges 
from  a  given  location  on  the  second  polarizer  (hereafter 
called  the  analyzer).  The  form  of  the  resultant  rela- 
tionships will  determine  the  parameters  and  functions 
which  are  necessary  to  compute  the  location  of  fringes 
for  a  given  flow  field* 

The  discussion  is  in  three  parts.  The  first  considers 
flows  in  which  the  velocity  may  be  regarded  as  constant 
along  any  given  light  path.  Such  flo'v^/s  will  be  designated 
as  two-dimensional  and  have  a  direct  parallel  in  the  two- 
dimensional  models  analyzed  by  the  traditiorial  methods  of 
photoelasticity.  The  second  part  of  the  discussion  will 
consider  steady  flov/s  in  vrhich  the  velocity  varies  along 
the  light  paths.   It  vrill  bs  shovm  that  the  results  of  this 
throe-dimensional  analysis  reduce  to  these  of  the  two- 
dimensional  case  vjhen  the  limiting  case  of  negligible 
varia.tion  along  each  light  path  is  considered.  The  final 
part  of  this  chapter  lists  the  successive  steps  to  be 
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carried  out  in  turn  to  obtain  the  fringe  pattern  from  a 
hypothetical  velocity  distribution  in  the  cross-section. 

Tyro-Dimensional  Flov: 
Ihe   analysis  of  tvro-diciensional  floirs,   of  bjjrefringent 
fluids  occurs  in  many  places.  An  example  is  Thur stone  and 
Schrag  (1962).  Consider  Figure  2  in  which  a  txTO-dimensional 
flo\7  field  in  the  yz-plane  is  observed  by  polarized  light 
moving  through  the  flow  field  in  the  positive  x-direction. 
Neglecting  attenuation,  the  electric  field*  has  the  form 

E  =^  Eo  cos(2it:.:/A)  (-sin  y  5  +  cos  y  k) 
where  Eq  is  the  amplitude  of  the  wave,  A  is  the  wavelength, 

A 

and  Y  is  the  angle  which  the  polarizer  makes  with  k. 

The  propagation  of  the  vector  E  is  a  function  of  its 
direction  in  the  birefringent  medium.   Specifying  the  E 
directions  for  maximum  and  minimum  propagation  speeds  by  ni 
and  ti2  respectively  and  assuming  that  n^  .cij  =  0,  the  vector 
E  can  be  resolved  intc  components  parallel  with  rii  and  no: 
E  =  Eo  [  cos(2TTntx/A)  cos(Y-y)  ni  - 

cos(2TTn2x/A)  sin(Y-^y)  in-2  ]  • 
Here  -^^   is  the  orientation  angle  sho\m  in  Figure  2,  and  n-j 
and  n2  are  the  refractive  indices  of  the  m.edium  v;hen  E  is 
parallel  with  ni  and  na  respectively. 

IJhen   the  light  emerges  from  the  flow  field  at  x  =  b 
the  only  light  which  will  pass  through  the  analyzer  is 


*In  this  dissertation  vectors  are  designated  by  boldface  (E) 
and  unit  vectors  are  denoted  by  a  circumflex  (ni  ) . 
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that  component  v/hich  is  parallel  v;ith  n^   where 

ni  •  n^  =  cos  (y'  -  ^)  > 

n2  •  na  =  sin  (y'  -  i]r). 
In  that  case  the  emergent  amplitude  is 

Ee  =  E  .  £ia 
or 

Eg     =     Eq   [cos(2TTni  6/A)    cosCy-^]/)    cosCy'-ilr)      + 
cos(2TTn2^/A)   sin(Y->I/)   sinCy'-t)]    • 
Setting 

n  =  (ni  +  n2)/2 
An  =  n2  -  ni 
and  expanding  the  cosines  yields 

Ee     =     Eo   {  [cos(2TT6rx/A)    cos(7r6£in/A)   + 

sin(2TT6n/A)   sin(n-6/\n/A)  ]  cos(y-iA)cos(y' -v) 
+    [co3(2TT6n/A)    cos(TT'6An/A)   - 
sin(2Tr6"n,/A)   sin(TT6An/A)  ]   sin(y-iy)sin(y'-Tj,')  }  • 
This  simplifies  to 

Ee     ~     Eo  [  cos(2TT6n/A)   cos(TrFAn/A)   cos(y-y' )     + 
sin(2TT6n/A)    sin(iT6i:>n/A)    cosC'f^y'-2^|<-)]    . 
For  a  dark  polarizing  field,   y  -  y'    =  ~/2  v/ith  the 
familiar  result  upon  substitution, 

Eg     =     Eo   sin(2rr6n/A)    sin(TT6An/A)    sin  2(y->Tr) , 
which  is  used  in  two-dimensional  optical  stress  analysis 
of  solids    (e.g.  Dally  and  Riley,   1965,   P^  171). 

For  a  fringe   to  be  observed  it  is  necessary  and 
sufficient  that  an3'-  of  the  sine  factors   in  this  equation 
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equals  zero.  Each  of  these  conditions  will  be  discussed 
individually. 

The  first  factor  reflects  the  periodic  variation  in 
E  due  to  the  v/ave  nature  of  light.   The  associated 
frequency,  about  lO''''"  Hz,  is  too  rapid  for  eye  or  camera. 
For  this  reason  the  constant  amplitude  is  usually  replaced 
by  Eo  5  V7here 

Eo  =  Eo  sJJi(2Tr&n/A). 

The  middle  factor  is  responsible  for  the  colored 
fringes  Imovm  as  isociiromatics  which  are  seen  v;hen  the 
flow  field  is  illmainated  with  white  light.   In  mono- 
chiTomatic  light  the  same  namie  is  retained  for  those 
fringes  resulting  from  the  condition 

/^rib/x      -     IT 
where  N  is  an  integer  knovm  as  the  fringe  order. 

The  final  factor  is  responsible  for  the  black  fringes 
Iciovra  as  isoclinics  along  which  the  principal  optical 
a:>:es  are  parallel  with  the  polarizers.  This  condition 
is  knovrn  as  extinction  and  the  so-called  e-tlnction  angle 
is  defined  by 
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In  monochromatic  light  it  is  difficult  to  distinguish 
isochuromatics  from  isoclinics.  This  difficulty  may  be 
avoided  by  using  circularly  polarized  light  for  v/hich 
Eg  =  Eo  sin(TTAno/A) . 
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The  derivation  of  this  expression  is  straight.forv;ard,  but 
tedious,  and  it  may  be  found  in  many  references  including 
Dally  and  Riley  (1965,  £2..  174-179). 

Three-Diinensional  Flov/ 

The  eicperimental  stress  analysis  of  three-dimensional 
solids  is  accomplished  by  locking  the  deformation  in  place 
and  then  cutting  the  model  into  slices  for  tv;o-dimensional 
analysis.  The  availability  of  this  technique,  which 
provides  local  stress  distributions  along  any  desired 
path,  has  inhibited  interest  in  the  study  of  three- 
dimensional  fringe  patterns  as  such.  Some  studies  have 
been  carried  out  using  scattered  light  polariscopes  which 
have  the  effect  of  placing  a  temporary  analyser  or  polar- 
izer at  a  selected  plane  within  the  throe- dimensional 
model*  This  technique  is  described  by  Van  Daele-Dossche 
and  Van  Geen  (1969). 

In  liquids  the  direct  three-dimensional  analysis  of 
birefringent  patterns  to  obtain  velocity  fields  is  not 
feasible  due  to  the  variety  of  conditions  which,  in 
principle,  could  lead  to  the  same  fringe  configuration. 
On  the  other  hand,  there  seems  no  theoretical  objection 
to  the  inverse  method:  assuming  a  flow  distribution  and 
determining  the  resultant  fringe  pattern.  The  corres- 
ponding analysis  follows. 
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In  three-dimensional  flox.^s,  even  if  the   streamlines 
are  parallel,  the  principal  optic  axes  will  be  oriented  in 
throe- dimensional  space.  As  a  result,  the  directions  n-i 
and  ^2  must  he  recognized  as  lying  parallel  to  the  principal 
axes  of  the  ellipse  formed  by  the  intersection  of  a  tiiree- 
dimensional  ellipsoid  with  the  plane  orthogonal  to  the 
light  path  at  the  point  of  interest  (Sonmerfeld,  1964-,  pp. 
139-147).  Care  must  be  taken  not  to  confuse  the  directions 
rii  and  nz   with  the  pro;iected  axes  of  the  ellipsoid.  The 
form.er  are  orthogop-al;  the  latter,  in  general,  are  not.   In 
th-ree-djinensional  flows  the  directions  n^  and  n2  v/ill  vary 
along  the  light  paths  as  will  ni  and  n2,  the  magnitudes  of 
their  respective  refractive  indices. 
Assujnptions 

In  the  analysis  which  follows  three  assumptions  are 
mad.e  regarding  the  optical  ellipsoid.  These  are  discussed 
separately. 

Assmi-ption  1.  —  The  properties  of  the  optical  ellipsoid 
are  completely  defined  by  three  cha-racteristies:  the  diff- 
erence between  the  length  of  the  longest  axis  and  the  tv:o 
shorter  aires  (assumed  equal),  the  magnitude  of  the  incli- 
nation of  the  longest  azcis  to  the  principal  flow  direction, 
and  the  direction  of  that  inclination.  These  variables  are 
sho\'jn  in  Figure  3  and  are  denoted  by  An,  tV,  and  ©  respec- 
tively. The  first  two  are  easily  identified  with  the  bire- 
fringence and  the  orientation  angle  of  the  previous  section. 
The  third  parameter,  ®,  ^rill  be  identified  as  the  rotation 
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FIGURE  3»  Optical  ellipsoid,  showing  parameters  9y  ']/•, 
and  An, 
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angle  and  is  necessary  to  describe  variation  of  the  optical 
properties  as  the  direction  of  the  flov;  changes  along  the 
light  path.  The  magnitudes  of  the  ellipsoidal  a:^es  vary 
so  little,  An  «  n,  th:it  it  can  be  assimed,  as  nsu.al,  that 
the  mean  value  of  the  refractive  index  is  constant.  The 
previous  assumption  that  the  two  shorter  optical  axes  are 
equal  in  length  has  been  made  principally  in  the  interest 
of  economy  since  it  reduces  the  number  of  optical  para- 
meters to  a  manageable  number. 

Assumption  2.  —  The  optical  properties  are  uniquely 
defined  by  the  local  shear  rate,  iluj.ierous  authors,  among 
them  Truesdell  and  Noll  (1965)}  have  pointed  out  the 
critical  importance  of  history  in  the  description  of  the 
properties  of  a  material.   In  the  present  case  the  flow  is 
steady.   If  a  fading  memory  is  assumed  for  the  material, 
then  after  a  short  time  the  history  of  the  flov;  may  be 
neglected.  The  neglect  of  strain  in  the  definition  of  the 
optical  properties  i'ollovis   from  experiments,!  xrork  v/ith  the 
medium  which  indicates  that  the  elastic  properties  are 
negligible  except  at  very   low  rates  of  shear.  Work  in 
this  area,  principally  by  Thurstone  and  his  associates 
^Sj3jl   'thurstone  and  Sclirag,  1962),  iir.s   discussed  in  the 
previous  chapter.  A  unique  dependence  of  the  optical 
properties  upon  the  shear  rate  is  consistent  v;ith  those 
theories  vrhich  explain  birefringence  in  term.s  of  the 
continuous  rotation  of  microelem.ents  within  the  fluid 
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as  described  by,  say,  Boeder  (1932)  and  Kuhn  and  Kuhn 
(194-3).  Acceptance  of  these  theories  is  not  necessary 
for  acceptance  of  Assiinption  2  vjhich  is  not  contradicted 
by  any  experimental  evidence. 

Assiunption  3»  ~  The  inclination  of  the  longest  axis 
of  the  ellipsoid  occurs  in  the  direction  of  the  local 
veloci-^3;-  gradient.  That  is,  the  rotation  angle  is  given  by 


©  =  tan' 


-1 


5g75x 


where  g  is  the  magnitude  of  the  shear  rate.  The  coincidence 
of  principal  directions  for  certain  optical  and  Theological 
properties  has  been  suggested  by  Lodge  (1956)  and  has 
strong  heuristic  appeal. 
Definition  of  Effective  Optical  Protjerties 

The  properties  of  the  optical  ellipse  at  any  point  on 
a  light  path  vrithin  the  tliree-dimensional  flow  field  follow 
ijnmediately  from  the  assiornptions.  As  shovrn  in  Appendi;^  A, 
the  effective  birefringence  v/ill  be 

All  =   [(Ai2+A2^-^A3-)co3  2<I^-  2(AiA4+A3Aj)sJJi'^co3^ 
+  (A.2.-A52)sin^^]-^ 
-  [(Ai^+A2^+A3^)sin2-«{^'i-  2(Ai  A^+A3A5)sin  tcos"^ 
+  (A4.2+A;2)eos24-]  -V2 
where  "^  is  the  effective  orientation  angle, 


^    =    -   tan' 


-1 


2  (AiA^+A^A^) 

Ai2+A2  2+A32_A42_A52 


L^i 


and  the  variables  A^  through  A;;  are  defined  in  terms  of 
the  optical  parameters  An,  iV,  and  ®.  For  example, 
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Ai  -     cos  i},'-  sin  9  /  (ii  -  ~j-)  • 
Values  for  A2  tiiron.gh  Aj  are  given  in  Appendix  A. 

It  is  appropriate  to  consider  conditions  for  light 
paths  parallel  with  and  close  to  the  side  v;alls  where  ©  has 
a  limiting  value  of  tt/2.   Substitution  of  this  result  leads 
after  considerable  manipulation  to  the  result 

t^]^j     =  An  sin^2^y. 
At  lo\i   flow  rates,  when  it   approaches  tt/4-, 
A?!.,  =  An. 

VV 

iins.  lysis 

To  analyze  flows  in  vrhich  /^N  and  "^  vary  along  each 
light  path,  it  is  conv^enient  to  employ  the  Poincare  (l339) 
model  in  vrhich  any  change  in  the  condition  of  a  polarized 
light  beam  may  be  represented  'by   the  movement  of  a  point  on 
the  surface  of  a  unit  sphere.  Procedures  for  use  of  the 
Poincare  sphere  are  found  in  many  texts™  The  present  sign 
convention  follows  Hartshorne  and  Stuart  (1970)  and  is 
illustrated  in  Figure  4.   Each  point  on  the  sphere  is 
expressed  in  terms  of  angular  coordinates  4'  and  il  which 
correspond  respectively  to  latitude  and  longitude  in  terres- 
trial navigation.  '"V/est"  and  "Korth^'  being  the  positive 
directions. 

Two  points  are  designated  on  tlie  Poincare  sphere; 
P,  \rhich  represents  the  condition  of  the  polarized  light, 
and  H,  which  represents  tlie  principa.l  optic  axis  of  the 
medium.  These  points  arc  determined  as  follows. 
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FIGURE  4.   Poincare  sphere  showing  P,  the  condition  of 
the  polarised  light  beain,  R,  the  principal  (fast)  optic 
axis  of  the  mediiin,  and  AP,  an  arc  on  the  surface  repre- 
senting the  change  in  polarization  which  occurs. 
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Poijit  P.  —  In  the  general  case  polarized  light  is 
elliptically  polarized  having  the  general  form  (for  propa- 
gation along  the  x-axis) : 

E  =  As  cos  (cot  +  f  )  i:i  +  Bg  cos  (cot  +  -f  )  ^2  ? 
v/here  0',  a",  Ag,  and  Bg  are  constants,  only  two  of  which 
are  independent,  and  J.^    and  K2  are  orthogonal  cartesian 
unit  vectors  in  the  yz-plane.   In  the  present  case  the 
form  selected  is 

E  =  ~  [ii  cos  (cot  -  f )  +  ia  cos  (wt  +  f)]  , 
wiiich  yields  on  rotation  through  the  angle  (e  -  it/4)  : 


E  =  ^  [-[sin  (6  -  -J)  cos  (6ot  -  f )  - 


-  £0 

cos  (€  -  f )  cos  (cot  +  f)  ]  J 

+  [cos  (e  -  ■^)  cos  (cot  -  ■2)  "^ 

sin  (G  -  f)  sin  (cot  +  f)]  k}  . 
Further  trigonometric  manipulation  yields: 

E  =  Eo  [-(sin  e  cos  -f  cos  cot  +  cos  e  sin  ^   sin  f-it)  j 
+  (cos  e  cos  -f  cos  cot  -  sin  e  sin  -^  sin  (^t)  k]. 
This  form  of  representing  E  has  two  useful  properties 
(Figure  5).  The  variable  e  can  be  identified  with  the 
angle  v/hich  the  principal  axis  of  the  elliptically  polar- 

A 

ized  light  makes  v^ith  the  principal  flow  direction  [k..     The 
variable  a   reflects  the  eccentricity  of  the  ellipse  for 
v;hich  the  major  and  minor  axes  have  magnitudes: 


Further,  the  variables  e  and  cr  are  simply  related  to 
the  coordinates  of  the  Poincare  sphere.  Specifically, 
P(n,<5)   =  P(26,cr). 

The  dependence  of  the  electric  field  vector  upon 
time  is  usually  neglected  in  studies  utilizing  polarized 
light, and  texts  such  as  Hartshorne  and  Stuart  (1970) 
utilize  the  parameters  e   and  a   directly  without  stating 
E(t)  explicitly. 

Point  R.  —  Point  R  on  the  Poincare  sphere  represents 
the  optical  properties  of  the  medium.   If,  as  has  been 
assumed  in  all  previous  studies  using  Milling  Yellow  as 
the  birefringent  medium,  the  solution  acts  as  a  linear 
wave  plate,  then  the  point  R  is  located  on  the  "equator" 
of  the  Poincare  sphere  at 

R(il,$)   =  R(2?-,0) 
where  "^I^  is  the  effective  orientation  angle.  In  this  case 
point  R  represents  the  location  of  the  "fast"  optical 
axis  parallel  to  ni .   If  there  were  evidence  of  isomerism 
or  helicity  in  the  molecular  structure  of  Killing  Yellow, 
or  if  a  degree  of  circular  polarization  v;ere  observed 
under  conditions  of  zero  shear,  then  these  effects  would 
require  redesignation  of  point  R  at  that  value  of  $ 
corresponding  to  the  eccentricity  of  the  elliptical  wave 
plate  (see,  for  example,  Shurcliff  and  Ballard,  1964). 
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FIGURE  5.  Elliptically  polarized  light  in  the  yz-plane, 


r+Air 


FIGURE  6.   Projection  of  OP 
on  CR   and  definition  of  ir. 


FIGURE  7.   Definition 
of  Ar.  Angle  is 
2ttANax/X. 
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Once  points  ?  and  R  have  been  designated  the  change  in 
polarization  follows  directly.  The  Poincare  sphere  is 
constructed  such  that  the  change  in  P  due  to  a  ir.ediun  with 
properties  represented  by  Pi  follovrs  a  countor-cloci-r^.fise 
circular  path  about  the  axis  0R  with  an  included  angle  of 
2TT^.'yAx/A  relative  to  this  axis. 

To  obtain  the  associated  vector  equations,  designate 
the  cartesian  coordinates  of  the  sphere  by  Xi ,  X25  and  X3 
with  corresponding  unit  vectors  £1  ,  ^25  and  2^3.  Locate 
P  and  H  b;^  the  unit  vectors  IP  and  (R  whore 

R  =  ^  =  2Ci  sin  2^  +  2:2  cos  2<ir 
and 

P  =  (^  =  ;^i  cos  a   sin  2e  +  kz   cos  cr  cos  2e  + 

z:3  sin  a,  (3»1) 

Referring  to  Figure  6,  the  projection  of  OP  on  CR 
is  (!P»R)R  and  the  radius  of  the  path  is 

r  =  P  -  (P-R)  R. 
The  change  in  F  is  AP,  a  circular  arc  vrith  chord  /^r,  where 

^r  =  P  -  Pq. 
Since  the  chord  subtends  an  angle  of  2TTaxi/\x/A .  it  follows 
from  Figure  7  that  the  length  of  t>x:   is 

Ar  -     2r  sin(-A::£kx/A) 
and  its  direction  is  determined  by  the  conditions: 

R-AC   =   0; 

(r«Ajr     =     -r  Ar   sinCtrAII^c/A)  ; 
(f?xr)^iMr     =      (P-R)   r  Ar   cosCttAUax/A)  . 


54 


V/hen  the  indicated  vector  operations  are  carried  out, 
three  simultaneous  equations  result  from  which  the 
components  of  aip  are  found  to  be :  (3»2) 

cTi      =2  sin(T:l\KAx/A)  cos  2^ [sin  a   cos(ttANax/A) 

-  cos  cT  sin(TTANAx/A)  sin  2ie-'^)]   ; 

AT 2     =     -2  sinCiT Max/A)  sin  2^[sin  a   cosCttANax/A) 

-  cos  a   sin(TTANAx/A)  sin  2(e-->l0]  ; 

Ar3  =  -2  sin(T:ANAx/A)  [sin  a   sinCTiANAx/A) 
+  cos  a   cos(TrANAx/A)  sin  2(e--5r)], 
Consider  a  homogeneous  flow  of  thickness  ax  =  6 
having  optical  properties  AN  =  An  and  4''  =  ^f  which  are 
constant  along  any  given  light  path.  Upon  this  two- 
dimensional  flow  let  light  fall  which  is  plane- polarized 
(c  =  0)  by  passage  through  a  polarizer  oriented  such  that 
Y  =  e.   If 

An6/A  =  IT 
where  N  is  an  integer,  or 

r  =  y  -  "2 

it  can  be  seen  that 

-^ri   =  Ar2  =  Ar3  =  0 
and  the  polarization  of  the  light  is  the  same  ^vhen  it 
leaves  the  flo\/  as  when  it  enters.   In  a  "dark  polarizing 
field"  the  analyzer  is  oriented  at  right  angles  to  the 
polarizer,  and  for  such  a  field,  if  Air  =  0,  a  fringe  will 
result.  Necessarily  the  conditions  for  which  Atr  =  0  are 
identical  with  those  for  isochromatics  and  isoclinics  in 
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the  previous  section  on  two-dimensional  flow. 

For  more  general  two-dimensional  conditions  it  must 
be  recognized  that  there  will  occur  two  types  of  fringes: 
those  due  to  total  extinction  of  plane- polarized  light 
and  those  due  to  partial  extinction  of  elliptically 
polarized  light.  Besides  the  cases  already  considered, 
there  exists  only  one  other  condition  for  which  the 
light  is  plane-polarized  and  completely  extinguished  as 
it  leaves  the  flow.  This  occurs  when 

An5/A  =  II  +  ^ 
v;here  II  is  an  integer,  and  simultaneously 

t  =     i(AY  ±  ^) 
where  N  is  again  an  integer  and  hy   is  the  angle  between 
the  principal  directions  of  the  polarizer  and  analyzer. 
Note  that  both  conditions  must  be  satisfied  for  a  fringe 
to  be  observed. 

In  general,  light  emerging  from  a  flo\/  field  X'jIII  be 
elliptically  polarized.   In  the  discussion  of  how  point  P 
is  determined  on  the  Poincare  sphere,  the  electric  field 
vector  was  expressed  in  the  form: 

E  =  Ec  l_-(sin  e  cos  -f  cos  .->t  +  cos  e  sin  —  sin  Mt)oi 
+  (cos  c  cos  ■§•  cos  ^it  -  sin  e  sin  -j   sin  v.t)  kj  , 
Upon  passage  through  an  analyzer,  the  emergent  amplitude 
of  such  light  would  be  Eg  =  S.n^,  where 

cig_  =  -j;  sin  y'  "^  2£  cos  y'  ? 
or 
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Eg  =     Eo  [cos(€-y')  cos  tot  cos  -f  - 

sin(e  -y'  )  sin  cot  sin  -^  J  . 

The  intensity  of  the  light  varies  \;ith  the  square  oi 

the  electric  field  (see,  for  example,  Feyninan.  Leighton, 

and  Sands,  1963,  p.  31-10): 

I  =  kEgS 

and  the  average  value  of  the  intensity  for  periods  of 

time  much  larger  than  the  period  of  the  light  waves  will 

be: 

■=  _  liia  k  f*'   ^. 

1   —   t'-»<»  "TT  \   J-  Cl."C 

=   ^[C0S2(€-Y')  C0S2  f   + 

sin2(e-Y'  )  sin^  -^j   . 
The  conditions  for  a  fringe  when  the  light  is  plane- 
polarized  (cr  -   0)  have  already  been  discussed.   If  the 
light  is  circularly  polarized  (a  =  ±  ■^) ,   then 

T  =  kEo2/2, 
a  result  which,  predictably,  does  not  depend  upon  the 
characteristic  direction  m^  of  the  analyzer.  For  any 
ether  condition  of  light  leaving  the  fluid  there  will 
be  sone  analyzer  angle  y'  ''•^^hich  will  rnininiize  I.  This 
condition  may  be  obtained  formally  b:\'-  differentiating 
I(y')  with  respect  to  y'  ^'■-'^''-■-  setting  the  result  equal 
to  zero.  The  result. 
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is  consistent  v/ith  v/hat  would  be  predicted  from  looking 
at  Figure  5  and  can  be  used  to  determine  the  necessary 
condition  for  a  fringe  to  occur  when  the  emergent  light 
is  elliptically  polarized. 

For  tv;o-dimensional  flows,  since  the  entrance 
conditions  are  cr  =  0,  e  =  y?  the  initial  polarization 
is  given  by  substitution  into  equation  (3.1)  to  obtain 

A 

Po  =  X-,  Sin  2y   +  K2   cos  Sy. 

A 

The  change  in  P  is  Ar,  for  which  ^srhen  cr  =  0,  &  =  y? 
M   =  An,  "^  =  T}r,  and  ax  =  6: 

Ari   =  -2  siii2(TT^n6/A)  cos  2^|/  sin  aCy-^)? 

APa     =     2  sin^(TrAn5/A)    sin  2\Ir  sin  2(y-\|r), 

Ar3      =     -2  sin(iTAn6/A)    cos(TTAn6/A)    sin.  2(y-^|r). 
On  emergence  the  condition  for  a  fringe  requires  that 

P(e,cr)     =     P(-|  +  y',cr). 
Setting 

P(|  +  V  ,a)     =     P(y,0)   +  Air 
and  solving  for   cr  yields 

a     =     sin"''    [-sin(2TTAn6/A)    sin  2(y->]/)]  , 
which  defines   the  ellipticity  of  the  emergent  light,   and 


tan^CTTAna/A)      =       -sin  2(T'-y) 

sin  2(v'-H^-2\I/ 


y-2^) 

which  is  the  condition  for  a  minimum  to  occur  when  the 
light  is  elliptically  polarized. 

It  has  been  shown  in  the  previous  paragraph  that 
the  Poincare  sphere  can  be  used  to  obtain  a  complete 
description  of  the  polarization  of  light  passing  through 


po 


a  two-dimensional  flou.  These  relationships  will  now  be 
used  in  an  iterative  scherfie  to  determine  the  polarization 
of  a  light  beam  passing  along  a  path  for  which  the 
effective  optical  properties  are  knoim  esplicitly  but  are 
no  longer  constant. 

Integration  to  Obtain  Fringe  Patterns 
The  polarization  of  a  light  beam  moving  through  a 
birefringent  mediiom  has  been  obtained  in  terms  of  a 
position  vector  IP  which  designates  a  point  on  the  surface 

A 

of  the  Poincare  sphere.  The  incremental  change  in  IP  may 
be  expressed  in  terms  of  the  chord  Ar : 
[P(x+Ax;e+Ae,cr+ACT;£kri,'9')   = 

A 

P  (x ;  €  5  0 ; AIT  ,^)  -^  act  ( Ax ;  e ,  cr  jAN  ,-<^) . 
By   choosing  values  of  ax  small  enough  so  that  AN  and  '^ 
may  be  regarded  as  constant  for  the  increment,  the  change 
in  IP  along  the  entire  light  path  may  be  obtained  by 
summation: 

IP(x=6i  )  =  P(x=-6i)  +  i^r^. 
The  effective  parameters  AIT  and  "^  are  functions  of  An, 
il-,  and  ©  which  vary  in  turn  with  the  local  shear  rate 
g(x,y),  itself  a  function  of  the  spatial  coordinates  x 
and  y.  Figure  8  provides  a  schematic  representation  of 
the  means  by  vrhich  the  summation  is  to  be  performed.  The 
corresponding  steps  are  tabulated  below: 

1.  —  Choose  the  y-coordinate  of  the  licht  beam.  Set 
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Figure  8.   Schewatic  representation  of  successive  deter- 
minations of  the  variables.  The  subscript  s  has  been 
omitted  frora  the  variable  names.   The  notation  (0)  indicates 
that  all  variables  except  y  are  reset  at  their  initial 
(x  =  -61 )  values.   Question  marks  indicate  decisions. 
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2.  —  Determine  the  polarization  P  as  the  light  beam 
enters  the  flov;  hj   setting  6  =  y  a^^  c  =  0  in  equation 

(3.1). 

3.  —  Obtain  the  shear  rate  gg^Xgjy)  from  a  distribution 
calculated  from  those  in  Chapter  Six  or  elsewhere.   The 
subscript  indicates  that  this  is  the  s   iteration. 

4.  —  It  has  been  assumed  that  the  rotation  angle  0g 
is  defined  by 


®g  =  tan" ^ 


5gs/5y 


Obtain  ©g(gg) . 

5.  —  Obtain  AngCgg)  from  the  optical  relationships 
of  Chapter  Four  or  elsei'/here.  Obtain  VsCgs)  i^  a  similar 
manner. 

6.  -  ""J^gCAngj^^gjOs)  and  /ii^s(Ans,tg,®g)  are  defined  in 
Appendix  A.   Obtaiji  "^g  and  Allg. 

2.  ~"  Choose  a  trial  va.rae  for  AXg* 

8.  —  Repeat  steps  3  through  7  with  Xs-m  =  Xs  "^  ^^s 
to  obtain  the  corresponding  terms  with  subscripts  s+i . 

9.  -  Sat  ^=  (%+i  +  ^s)/2. 

10.  -  Set  AIT  =  (AHs+i  +  AHs)/2. 

11.  —  Determine  the  fractional  variations  j  C'^-^c)/'^ 

and  (AII-AH5)/AiI,  of  the  optical  coefficients  for  the  interval 
AXg.   If  either  variation  exceeds  a  prescribed  level,  say 
f  =  0.01,  reduce  the  value  of  AXg  and  repeat  steps  8  tlirough 
11.   Omit  this  stcT)  if  AN  or  "^  =  0. 
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12.  —  Obtain  Ar^  from  equation  (3«2). 

13.  -  Set  Ps-t-1  =  Ps  "^  ^-^^s' 

14.  —  From  the  definition  of  P  the  2C3-component  is 

P3     =     sin  a. 
Hence , 

cTs+1      =     sin"^  [(^3)3-5-1]  • 
Obtain  crs+i  . 

15.  -  From  the  definition  of  P  it  is  also  clear  that 
£5+1      =    \  tan"UPi/P2)s+i    . 

Obtain  £3+1 . 

16.  —  If  AXg  +  Xs4-i  <  5i ,  repeat  steps  8  through  1?. 
If  AXg  +  Xs-H  >  5i  J  set  AXg  =  61  -  xs-fi  and  repeat  steps 
8  through  15.  If  ax^  =0,  go  on  to  step  17. 

17.  -  Calculate  the  relative  intensity  of  the  light 
emerging  from  the  flow  field  and  analyzer  at  coordinate  y. 
Recall  that: 

^       =  cos2(t-.Y')Gos2(cT/2)  +  sin2(t-Y')sin2(0/2). 

18.  —  Choose  a  new  y-coordinate  and  repeat  steps  2 
thj?oiigh  17. 

19.  —  When  the  final  y-coordinate  has  been  chosen 

(y  =  62)  and  the  final  relative  intensity  has  been  calcu- 
lated, plot  I/Io  as  a  function  of  y  to  obtain  the  fringe 
pattern  predicted  by  the  relationships  chosen  in  steps  3 
and  5. 


CHAPTER  FOUR 
DSTER14INATI0N  OF  OPTICAL  PROPERTIES 

This  chapter  discusses  the  determination  of  tiro 
optical  properties  of  Milling  Yellow:  the  birefringence, 
An,  and  the  orientation  angle  ^^r.  Both  of  these  quantities 
have  been  defined  in  Chapter  Three.  They  will  be  discussed 
separately,  the  orientation  angle  only  briefly  at  the  end. 

Birefringence 

The  birefringence  is  that  property  of  an  optically 
active  fluid  which  is  responsible  for  that  lart  of  the 
fringe  pattern  knovm  as  isociiromatics.  The  necessary 
condition  for  these  fringes  is  that 

AnsA  =  K 
where  An  is  the  birefringence,  6  is  the  thickness  of  the 
flow  field  through  v/hich  the  light  passes,  A  is  the 
wavelength  of  the  incident  light,  and  N  is  an  integer 
known  as  the  fringe  order.   It  is  clear  that,  except  for 
a  constant  of  proportionality,  the  birefringence  and  the 
fringe  order  niay  be  used  interchangeably  for  a  given 
experimental  arrangements 

The  data  of  Peebles,  Prados,  and  Honeycutt  (1965) 
indicate  that  at  a  given  concentration  and  temperature 
Milling  Yellow  has  a  birefringence  (fringe  order)  which 
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shoi.;s  a  progressive  non-linea.r  increase  with  shear  rate. 
This  increase  suggested  that  the  data  be  replotted  with 
the  square  of  the  fringe  order  as  the  dependent  variable. 
Figure  9  shows  the  result  when  this  is  done  for  the  four 
most  concentrated  solutions  reported  by  Peebles,  Prados, 
and  Honeycutt.  The  hypothetical  relationship 

K2  =  l:ig  +  k2 
appears  adequate  for  the  range  of  shear  rates  shown. 

An  experimental  verification  of  this  form  was  carried 
out  with  the  Milling  Yellow  used  in  the  present  disser- 
tation.  In  contrast  to  the  data  of  Peebles,  Prados,  and 
Honeycutt,  which  were  measured  in  a  concentric  cylinder 
polariscope  at  nearly  constant  shear  rates,  the  data  for 
the  present  study  vrere  obtained  near  the  wall  of  a  nearly 
square  rectangular  conduit  which  is  described  in  the  next 
section.  In  the  succeeding  sections  will  be  found  an 
analysis  of  Milling  Yellow's  optical  response  when  flov;ing 
slowly  thJ70ugh  such  a  channel,  some  experimental  data,  and 
a  discussion,  thereof. 
Apparatus 

The  Milling  Yellow  solution  v/as  prepared  as  described 
in  Appendix  B.  Ihe   nearly  square  conduits  used  in  these 
experim.ents  have  been  described  previously  by  Lindgren 
(1962.  1963)  v/ho  constructed  them  to  observe  the  trans  it  ion 
betvreen  lam.inar  ajid  turbulent  flows  in  bentonite  sus- 
pensions. A  schematic  of  the  apparatus  is  shown  in 
Figure  10. 
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Overhead  tank 


Rectangular 
conduit 


Light 
source 


3  Punp  and  motor 


FIGUEIE  10.   Sch'^^natic  of  experimental  at)r)3Tatus .  The 
polarizer  and  analyzer  are  represented  "by  vertical  lines 
to  the  right  and  left  of  the  rectangular  conduit  between 
the  light  source  and  the  camera. 
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The  Milling  Yellow  was  normally  stored  in  a  lovrer  tanl: 
which  is  lined  v;ith  epo:-;y- bound  fiberglass.  From  the  taiilc 
there  is  a  gravity  feed  into  a  Mo;/no  EM- 2304  special  appli- 
cation pump  which  provides  positive  displacement  with 
minimum  shearing  of  the  fluid. 

At  intervals,  between  experim.ental  runs,  the  stopcock 
above  the  pump  was  opened  and  the  covered,  polystyrene, 
overhead  tank:  was  filled.  When  pumping  was  complete, 
normally  a  matter  of  seconds,  the  pump  was  turned  off  and 
the  stopcock  closed  to  prevent  siphoning  thiTough  the  puj^ip 
or  the  introduction  of  air  into  the  pipe. 

The  fluid  level  in  the  overhead  tanlc  can  be  maintained 
at  a  constant  head  by  pumping  continuously  and  permitting 
the  excess  fluid  to  return  to  the  lower  tani:  through  an 
overflow  pipe;  hov/ever,  this  procedure  results  in  undesir- 
a.ble  temperature  rises  and  to  the  entrainment  of  air  in 
the  Milling  Yellow.  VJhen  it  v;as  found  that  the  fluid  level 
in  the  overhead  tanli  had  a  negligible  effect  upon  flow  rate 
through  the  rectangular  conduit  (see  Appendix  C),  use  of  the 
overflox-;  pipe  was  restricted  to  providing  protection  against 
accidental  overfilling  of  the  upper  tank. 

A  gravity  feed  from  the  overhead  tanic  leads  to  a  set 
of  parallel,  vertical  conduits,  only  one  of  which  was  used 
at  a  given  time.   The  construction  of  these  conduits  is 
shovm  in  a  cut-away  view  in  Figure  11. 
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FIGURE  11.  Cut-av;ay  view  of  rectangular  corxduit 
showing  roughened  walls,  gaskets,  and  spacing  wires, 
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The  nearly  square  conduits  are  constructed  by  sand- 
wiching two  square,  12-iiijn  plexiglass  rods  between  two 
plexiglass  strips  having  widths  of  36  mm.  Plastic  gaskets 
in  V- joints  seal  the  flow  channel,  and  the  assembly  is 
bolted  together  along  the  length  of  the  conduit.  Thin 
wires  run  parallel  with  but  outside  the  gaskets  to 
maintain  constant  internal  dimensions.  On  certain  of 
the  channels  two  facing  surfaces  were  covered  with 
grinding  cloth  to  provide  a  kno\m  roughness.  The  conduits 
are  4-. 87  meters  in  length  with  the  internal  dimensions 
sho\m  in  Table  I. 

TABLE  I 
INTERNAL  DIMENSIONS  OF  CHANNELS  (Lindgren,  I963) 

Channel    Height  of  Roughness    Distance,  mm,  between: 
Number         Elements         Strips,  25i  Rods,  2b 2 

1  Polished  plexiglass      13.^3      13.26 

2  0.03b'  to  0.044  mm       13.^3      12.40 
6        0.59  to  0.70  mjn        13.^2     11.92 

Because  of  the  roughness  elements,  the  flow  can  be 
viewed  through  the  side  walls  only  in  the  clear  charjiel, 
and  even  in  this  case  the  x'iew  is  unss.tisfactory  due  to 
the  gaskets  which  prevent  a  viei/  of  flov;  along  the  walls. 
Tiirough  the  front  x.'alls  there  is  an  unimpeded  viev;  of  the 
flow  in  all  of  the  channels.  Nikuradse  (1933)  and  Moody 
(1944)  found  that  surface  roughness  plays  no  significant 
role  in  the  lajninar  flow  region,  and  Lindgren  (1963)> 
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using  the  present  apparatus,  reported  no  significant 
difference  in  flov/s  through  tubes  with  two  walls  roughened 
and  flows  through  tubes  v;ith  four  v;alls  roughened  v;hen 
the  transition  region  v;as  studied.  Hence,  the  principal 
effect  of  the  rougliness  elements  in  the  present  dissertation 
is  the  reduction  in  cross-sectional  area  which  results. 
Returning  to  Figure  10,  note  that  each  chamiel  was 
provided  with  nanoneter  taps  on  the  smooth-faced  sides  of 
the  channel.  The  nanometer  fluid  \-i3.s   carbon  tetrachloride 
(specific  gravit;^'-:  1.584). 

The  f lov;  x^ras  normally  illuminated  by  a  sodium  vapor 
lamp  from  which  light  passed  through  an  array  of  polari^iers, 
through  the  transparent  channel  and  flow  field,  and  t.hrough 
an  array  of  analyzers.  The  orientation  of  the  sets  of 
elements  in  the  thjree  array  configurations  used  during  the 
experiments  is  given  in  Figure  12  and  Table  II.  The  method 
of  aligriment  is  given  in  Appendix  D.   On  four  occasions 
(runs  45,  418,  419,  and  424)  the  arrays  v;ere  replaced  by 
circular  polarizers  differing  in  phase  by  n/2   radians. 

The  fringes  were  photographed  with  an  Exacta  Varex  I la 
35  am  camera  with  a  Jena  58   mm  lens.  A  bayonet  extension 
was  manufactm-ed  for  the  lens  permitting  an  object  23  cm 
from  the  film  plane  to  be  focused  upon  and  photographed. 
Kodak  Tri-X  film  \:as   used  and  commercially  doveloiDed  a.t  an 
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TABLE   II 


VALUES   OF  Y  Al^iC  y'    FO^   SUCCESS ITE  ELEMSIITS 
OF  POL/>RIZIKG  ARRAYS 


Conf igiiratiorx  for 
runs   123    &    2l8 


Y 

I 

Y 

0 

Tr/2 

-tt/2 

0 

0 

0 

0 

tt/2 

-tt/2 

0 

0 

0 

0 

tt/2 

-u/2 
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Conf ipiiration  for 
runs  130  &  131 


I 

0 

I 

n/2 

Y 

-7T/2 

0 

-n/4 

tt/4 

0 

C 

0 

0 

tt/4 

-n/4 

0 

tt/2 

-tt/2 

0 

tt/4 

-tt/4 

0 

0 

0 

0 

-:t/4 

Tr/4 

0 

Tr/2 

-tt/2 

Configuration  for 
all  other  runs  with 
plane- polarized 
arrays 


Y 
0 

-  tr/4 

tt/2 

tt/4 

0 

0 

0 

0 

0 


tt/2 

Tr/4 

0 

-tt/4 

Tr/2 

Tr/4 

0 

-tt/4 

Tr/2 


Y 
-tt/2 

Tr/2 

tt/2 


-n/2 


-tt/2 

-tt/4 

0 

tt/4 

-tt/2 


*Circular  -noiarization  used  during  rujis  45?    4l8,    419,    and 
424. 


72 


equivalent  exposure  index  of  ASA  I6OO,  A  typical  exposure 
was  a  lens  setting  of  f3.5  vrith  a  shutter  speed  of  1/1 50 
second. 

The  choice  of  flov;  channels  and  the  rate  of  flow 
were  controlled  by  a  stopcock  at  the  base  of  the  flow 
channel.  Flow  rates  were  measured  by  collecting  the 
fluid  at  the  outlet  above  the  lower  storage  taiik   and 
timing  the  efflux  to  the  nearest  0,1  second  v;ith  a 
Junghans  timer.   In  a  few  instances  the  efflux  volume 
was  measured  to  the  nearest  milliliter  with  a  250-ml 
graduate,  but  in  general  the  efflux  was  caught  in  a  clean, 
dry  vessel  and  its  mass  measured  on  a  Harvard  Trip 
Balance  to  the  nearest  0,1  gm.  Prior  to  each  Xveighing 
the  temperature  of  the  efflux  was  measured  to  0^1  C°   with 
an  ordinary  laboratory  thermometer. 
Procedure 

The  fluid  was  circulated  in  the  system  for  about 
twenty  minutes  to  insure  uniformity  of  temiperature  and 
concentration  bet^veen  the  upper  and  lower  tanlcs,  VJhen 
this  v/as  accomplished,  as  indicated  by  constant  temper- 
ature readings  in  the  lower  container,  the  camera  was 
cocked  so  that  the  shutter  would  release  12  seconds  after 
the  preset  mechanism  was  actuated.  The  lower  stopcock 
was  adjusted  until  the  desired  flow  was  obtained  as  indi- 
cated by  the  birefringent  fringe  order  at  the  wall  of  the 
channel.  The  efflux:  container  was  placed  under  the  fluid 


73 


outlet  and  simultaneously  the  timer  was  started.  After 
an  appropriate  interval  the  preset  mechanism  on  the 
camera  was  tripped  so  as  to  release  the  shutter  at  the 
midpoint  of  the  efflux  interval.  The  time  at  which  the 
shutter  was  heard  was  noted,  and  at  approximately  twice 
this  time  the  container  w^as  removed  and  the  timer  stopped. 
The  temperature  of  the  efflux  was  measured,  and  the  flow 
rate  checked  by  confirming  that  the  fringe  order  at  the 
wall  remained  unchanged.  This  done,  the  flow  was  stopped 
and  the  mass  of  the  efflux  measured.  Finally,  the  fluid 
container  v;as  cleaned  and  dried,  and  the  camera  and  timer 
readied  for  the  next  run.  For  the  last  t^^ro  runs  manometer 
readings  were  talcen  just  prior  to  cutting  off  the  flow. 
As  necessary,  the  overhead  tanl:  '.•ra.s  refilled. 
Prediction  of  Results 

Consider  the  steadj^  flov;  of  a  fluid  in  a  rectangular 
conduit  having  dimensions  of  26i  and  262*  A  force  balance 
yields 

(26-!)(262)  dp  =  \-(^6t  +  ^^2)   <iL 
where  dp/dL  is  the  pressure  gradient  along  the  conduit 
and  %j   is  the  average  shearing  stress  along  the  perimeter. 
For  a  IIev;tonian  fluid  in  a  rectangular  conduit,  Cornish 
(1928)  showed  that 

dp/dL  =     26iU>ij/62Sc 
v/here  u  is  the  mean  flovj   velocity,  yiy^   is  the  viscosity, 
and  So  is  a  geometric  constant  having  the  dimensions  of 
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areaJ 

2 


>o  - 


25i52   12862   nT  „-5 

3  T75      g^,_3 


>   s~^  tanii(sTT6i/262)  • 


The  dimensions  &i  and  62  i^ay  be  interchanged  in  the  last 
two  equations  without  affecting  anything  but  the  rate  of 
convergence  of  the  series.  Nunerical  values  for  So  are 
tabulated  along  with  some  other  constants  in  Appendix  G. 

Define  the  apparent  viscosity  of  a  non-Newtonian 
fluid  by 

y^^     =  (62S0/2&1U)  dp/dL. 
Substitution  into  the  force  balance  yields 

^^,  =  26i2)i^u/So(6i+62). 
Replacing  the  average  shear  stress  at  the  xvall  by  the 
product  of  the  average  shear  rate  and  the  average  viscosity 
there, 

g^^  =  26i2)aaU/So(6i+62)JI^^. 
The  viscosity  is  lowest  at  the  wall  where  the  shear 
rate  is  highest.  Replacing  )x^   by  JL,,^  -  ajx,   obtain 

g^  =   [26i2-a/So(6i+62)](l  -  a;i}I^,'''). 
3y  hypothesis, 

N2   =  kig  +  k2, 
so  at  the  wall: 


w 


2ki  61  ^u/So( 61+62 )1(1-A;aja^,^~  )  +  k2  .   (^.1) 


If,  a.s  Peebles  and  Liu  (1965)  have  sho>m,  the  viscosity  is 
nearly  constant  at  low  flow  rates,  aji  may  be  neglected  and 
H^^2  =  Z^.i  61  %/So(  61+62)  +  k2. 
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The  mass  flow  rate  is 

G  =  4p6i62^> 
where  p  is  the  density,  so  at  low  flow  rates: 

N.^^   =  ki6iG/2So62p(6i+62)  +  ^2.  (4.2) 

As  G  rises  the  assumption  that  a}i  is  negligible  becomes 
invalid  and  the  slope,  dTI^^VdG,  should  decrease. 

In  the  vicinity  of  the  wall  the  flow  approaches  two- 
dimensionality,  for  which  isochromatic  fringes  v;ill  be 
observed  when 

V;         vl 

Substitution  and  rearrangement  lead  to  the  result 

86162S0  (6t  +  52)      4-61'^ 

vrfiere  k-j  and  k2  are  to  be  determined  experimentally  at 
the  low  flow  rates  for  v/hich  the  relationship  is  valid. 
Experimental  Data 

Tests  of  the  equation  just  given  were  conducted  with 
the  results  shovm  in  Tables  III  to  J.VI  and  Figures  13  to 
16,   It  will  be  noted  that  the  temperature  varied  during 
these  tests. 

The  run  numbers  reflect  the  date  in  1973  upon  which 
the  data  were  talven.  Run  420A,  for  example,  was  the 
first  run  conducted  on  April  20,  Run  420D  was  the  fourth 
rim  on  the  same  day.  Although  data  vrare  obtained  for  two 
other  runs,  4203  and  420C,  these  data  are  not  presented 
because  the  flov;  rates  were  too  high  for  the  fringe  order 
at  the  v;all  to  be  accurately  determined. 
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TABLE  III 
FRINGE  ORDER  AT  WALL  -  RUN  123 


Fringe 
Order"^ 

Efflux, 
gm 

Time, 
sec 

Flov/  rate, 
gm/sec 

Tenp, 

2+ 

1590.8 

15.6 

102.0 

22.0 

2 

1505.2 

15.3 

98.4 

21.9 

1-V2+ 

1289^9 

18.0 

71.7 

22,0 

1-1/2 

14-53.2 

20.8 

69.9 

22.3 

2-1/2 

1076.7 

9.6 

112.2 

21.8 

1+ 

1321.2 

29.1 

45.4 

22,2 

1/2+ 

1436.3 

59.9 

23.98 

22.2 

1 

1391.5 

46.3 

30.05 

22«1 

V2 

1322.1 

151.1 
TABLE   iV 

8.75 

22.2 

FRINGE   ORDER  AT  WALL 

-  RLII  130 

Fringe 
Order* 

Efflurt, 
gn 

Time, 
sec 

136.2 

Flow  Rate, 
gm/sec 

8.19 

Temp, 

1 

1115.7 

20.0 

2+ 

937.0 

36.7 

25.53 

19.7 

2 

1251.1 

45.4 

27.56 

20.3 

3 

1032.8 

15.8 

65.37 

20.1 

^Halved  fringe  orders  were  so  identified  xdien  a  fringe 
occurred  at  the  x^^all  in  a  light  polarizing  field 
(analyzer  parallel  with  polarizer).  A  plus  (+)  indicates 
distinct  separation  of  the  fringe  from  the  wall. 
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TABT,^,  V 

i<'hii:ge  op 

'.DER  AT  WALL 

-  RUIT  131 

Fringe 
Order* 

Efflux, 

gm 

Time, 
sec 

39.1 

Flov;  rate, 
gm/sec 

26.32 

TenT3, 

1 

1029.0 

22.8 

2 

1555.5 

21.6 

72.02 

22.4 

2 

765.0 

25.2 

30.36 

21.2 

2 

826.4 

27.5 

30.05 

21.1 

2 

678.2 

22.6 

30.01 

21.0 

1 

695.7 

57.3 

12.14 

— 

2+ 

718.7 

18.2 

39.49 

21.0 

3 

1493.5 

22.5 

66.38 

21.0 

3-1/2 

17^9.2 

14.2 

123.18 

21.0 

3-1/2 

1434.6 

10.8 

132. S3 

20.8 

(+) 

395.8 

190.5 

2.088 

20.6 

1+ 

885.C 

69.8 

12.68 

20.6 

TABLE  VI 
ORDER   0?  PJED  FRINC-S  AT  V/ALL  -  RUII  2l8 


Fringe 
Order 

Efflu>:j 

gm 

Time, 
sec 

35.8 

Flow  rate, 
gn/sec 

Temp 

1 

199*>i< 

5.56 

23.8 

2 

218** 

19.8 

11.01 

23.9 

3 

212** 

11.7 

18.12 

22.8 

4 

1001.2 

20.5 

48.84 

22.3 

4(?) 

846.2 

18.3 

46.24 

22.2 

*See  note,  Tables  III- IV. 

**Killiliter  measurement  converted  to  grams 
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TABLE  VII 
HIINGE  ORDER  AT  WALL  -  RIBI  42 


Fringe 
Order 

EffltD:, 
gm 

Time, 
sec 

28.1 

Flow  rate, 
gm/sec 

0.736 

Temp, 

1 

20.7 

24.3 

io5 

29.8 

22.3 

1.337 

24.6 

2.3 

42.1 

17.4 

2.416 

24.7 

2.45 

108.8 

33.7 

3.23 

24.8 

3.5 

122.6 

30.5 

4.02 

24.8 

4 

131.0 

25.5 

5.14 

25.2 

5 

145.1 

23.8 

6.10 

24.7 

TABLE  VIII 
FRINGE  ORDER  AT  WALL  -  RUIT  44 


Fringe 
Order 

Efflux, 
gm 

Time, 
sec 

33.5 

Flow  rate, 
gm/sec 

2.04 

Tent), 

2.4 

68.3 

23.8 

3 

80.8 

33.2 

2.43 

24.2 

2.9 

80.0 

34.1 

2.35 

24.1 

3.2 

99.7 

33.0 

3.02 

24.9 

1.5 

85.3 

59.8 

1.43 

24.7 

4 

113.3 

29.9 

3.79 

24.7 

5 

163.35 

30.4 

5.37 

24.3 

5.35 

130.4 

20.0 

6.52 

24,3 

6.4 

219.1* 

24.0 

9.13 

24.2 

^Milliliter  measurement  converted  to  grans 
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TABLE  IX 
FRINGE  ORDER  AT  WALL  -  RIJIT  44A 


Fringe 
Order 

Effl-ox, 
gm 

Time, 
sec 

34.9 

Flow  rate, 
gni/sec 

1.194 

Temp, 

1.35 

41.7 

24.7 

l.S 

57.2 

34.8 

1.643 

24.9 

2.2 

68.2 

29.3 

2.288 

25.1 

2.35 

105.4 

36.3 

2.903 

25.2 

3.35 

66.3 

14.8 

4.51 

25.2 

5 

132.2 

20.4 

6.48 

25.1 

5.5 

175.4 

20.7 

8.47 

24.3 

6.5 

221.1=i< 

17«0 

13.01 

25.1 

TABLE  X 
FRINGE  ORDER  AT  WALL  -  RUIT  45 


Fringe 
Order** 

Effl-LDC, 

gra 

Time, 
sec 

100.9 

Flov;  rate, 
gm/sec 

0.191 

Temp , 
°C 

1 

19.3 

*** 

2 

16.6 

55.2 

0.301 

*** 

3 

52.9 

71.1 

0.738 

20.2 

3 

56.5 

50.3 

1.123 

20.6 

4 

66.6 

44.7 

1.490 

20.7 

6 

90.01 

39.4 

2.23 

20.7 

6 

109.5 

40.2 

2.72 

20.6 

7 

110.7 

30.1 

3.68 

20.7 

*Milliliter  measurement  converted  to  grams 

'^'-^First  run  v;ith  circularly  polarized  light.  Fringe  orders 

controlled  more  carefully  than  usual. 

***Efflu:<:  volume  too  small  for  tcm-osrature  measurement 
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TABLE  X  (Continued) 


Fringe 
Order 

EfflllK, 

gm 

TlTTie, 
sec 

29.4 

Flow  rate, 
gm/sec 

Temp, 

7 

153.9 

5.23 

20.8 

10 

211.1* 

26.4 

7.99 

20.8 

11 

272.4* 

18.2 
TA3T.F,  XI 

14.96 

20.4 

FRINGE  ORDER  AT  WAT.T, 

-  RUII  418 

Fringe 
Order 

Efflicc, 
gm 

Time, 
sec 

44.0 

Flow  rate, 
gm/sec 

Temp, 

1.27 

65.7 

1.49 

24.7 

2 

92.6 

36.1 

3.90 

24.8 

2.25 

110.4 

25.5 

4.33 

24.8 

3 

134.5 

22.3 

5.90 

24.9 

3.32 

158.0 

20.2 

7.32 

24.8 

3.46 

156.0 

17.4 

8.97 

24,7 

4.3 

166.1 

15.2 

10.93 

24.6 

5.0 

l3l,0 

13.7 

13.21 

24.7 

5.6 

188.7 

11.3 

16.70 

24.8 

6 

222.4 

10.3 

21.59 

24.8 

♦Milliliter  measurement  converted  to  grams 
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TABLE  XII 
FRiriGS   ORDER  AT  WALL  -  RUN  419 


Fringe 

Efflux, 

Time, 

Flov;  rate. 

Temp, 

Order 

gm 

sec 
60.2 

gm/sec 
0.78 

°C 

1 

47.1 

23.7 

1.7 

128.9 

60.0 

2.15 

23.8 

2.0 

124.3 

39.7 

3.13 

24.1 

2.35 

94.2 

21.9 

4.30 

24.2 

3.28 

124.5 

20.0 

6.23 

24.2 

3.50 

135.9 

19.7 

6.90 

24.2 

4.30 

157.2 

18.1 

8.69 

24.1 

5.0 

168.6 

15.9 

10.60 

24.1 

6.0 

149.4 

11.6 

12.38 

24.1 

TABLE  XIII 
FRINGE   ORDER  AT  WALL  -  RUIT  420A 


Fringe 
Order 

Effl-ox, 

gE 

Time, 
sec 

60.0 

FlovT  rate, 
gm/sec 

I.S3 

Temp, 

1.45 

109.7 

23.2 

2 

55.7 

20.4 

2.73 

23.4 

2.36 

73.1 

20.0 

3-66 

23.8 

3.1 

100.7 

19.7 

5.11 

23.8 

3.4 

128.0 

19.7 

6.50 

23.3 

4.2 

118.3 

16.1 

7.35 

23.7 

4.3 

175.7 

21.4 

8.21 

23.7 

5.1 

159.7 

16.1 

9.92 

23. 3 

5.6 

192.5 

15.3 

12.18 

23.3 

6.3 

167.3 

11.2 

14.94 

23c3 
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TABLE  XIV 
FniNGS   ORDER  AT  Iv'ALL  -  RUII  420D 


Fringe 
Order 

Effl^o^c, 

Time, 
sec 

35.8 

Flow  rate, 
gm/sec 

2.34 

Tenp, 
°C 

1.28 

83.6 

25.5 

2.22 

127.7 

22.5 

5.67 

25.5 

3.3 

154.7 

15.9 

9.73 

25.5 

4.5 

191.2 

12.1 

15.80 

25.5 

5.4 

165.75 

8.3 

19.97 

25.5 

6.6 

185.22 

6.4 

28^94 

25.5 

7.8 

189.4 

4.8 

39.46 

25.3 

TABLE  XV 
FRINGE   ORDER  AT  V/ALL  -  RlHI  424 


Fringe 
Order 

Efflir:c, 
gm 

Time, 
sec 

40.2 

Flov;  rate, 
gn/sec 

Tenp, 

2.1 

100.7 

2.50 

24.9 

3.3 

163.7 

29.6 

5.53 

25.3 

4.3 

128.1 

20.2 

6.34 

24.4 

5.4 

134.6 

20.3 

9.09 

24.3 

6.0 

188.5 

14.7 

12.82 

24.3 

7 

211.9 

10.2 

20.73 

24.5 

8 

199.5 

7.8 

25.57 

24.2 
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Fringe 
Order 

1.58 

2.26 

2.43 

3.6 

4.4 

6 

7 


TABLE  XVI 
FRINGE  ORDER  AT  WALL  -  RUII  425 


Efflirc, 
gm 

Time, 
sec 

50.1 

Flow  rate, 
gn/sec 

1.99 

Temp 
°C 

99.8 

24.2 

141.0 

40.3 

3.50 

24.2 

146.0 

40.5 

3.60 

23.8 

141.6 

24.8 

5.71 

23.8 

153.3 

19.7 

7.78 

23.8 

175.2 

14.2 

12.34 

23.8 

185.6 

11.4 

16.00 

23.8 

84 


N. 


w 


:  CC 


Flow  rate,  grysec 


FIGURE  13.   Fringe  order  at  ■'.■Jail  bs   function  of  mass  flov/ 
rate.  Runs  :-0 -123 ;  ~© -131; -0-130 :-0-:a3. 
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Disci-isslon  of  Results 

This  discussion  is  in  three  parts.  The  first 
describes  the  preliminary  testing  used  to  determine  the 
feasibility  of  the  method.  The  second  part  describes  the 
manner  in  v;hich  the  data  \<fQre   interpreted.  The  last  part 
gives  the  numerical  values  \ihich  were  obtained  for  the 
optical  coefficients  used  to  calculate  the  birefringence. 
Preliminary  tests 

The  objective  of  these  tests  was  to  check  the  working 
order  of  the  apparatus,  refine  the  measurement  technique, 
and  confirm  the  validitj''  of  the  form  of  equation  (4.2). 
These  objectives  were  largely  fulfilled. 

Run  123  was  conducted  in  ciiannel  6.  The  data  for 
this  run  may  be  represented  by  a  straight  line  when 
plotted  as  shoxm  in  Figure  13 «  There  is  an  upward 
departure  suggested  by  the  last  datum  which  is  at  odds 
v/ith  the  falling  off  predicted  by  equation  (4-.1),  The 
upward  departure  may  be  due  to  temperature. 

Run  130  v/as  also  conducted  in  channel  6,  Ko  firm 
conclusions  are  possible  from  these  te\-i   data  except  that 
the  birefringence  is  higher  than  in  the  previous  test, 
a  result  consistent  with  the  lov/er  mean  temperature  and  a 
concentration  increase  due  to  evaporation. 

Run  131  I'ra.s  conducted  in  channel  2,  The  data  are 
well  represented  by  a  straight  line  at  low  flow  rates  and 
fall  off,  as  predicted,  as  the  flov;  rate  increases. 
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Run  218  was  conducted  in  channel  6  in  v/hite  light. 
This  v.'as  the  only  run  in  either  the  preliminary  or  the 
final  test  series  in  •'.-/hich  the  sodium  vapor  laiap  was  not 
used  for  illujnination.  Tvra  sets  of  isoch-Tomatic  fringes 
were  observed,  one  red  and  the  other  green.  Other  wave- 
lengths were  hard  to  detect  either  because  they  were 
strongly  absorbed  by  the  czediun  or  because  of  their  simi- 
larity to  the  deep  orange  of  I-illing  Yellow.  The  data 
plotted  in  Figure  13  are  those  from  the  red  set  of  fringes. 
The  data  are  too  few  to  serve  as  a  test  of  equation  (4.2). 

In  general  the  prelLminary  tests  support  the  analysis 
which  predicted  that  the  square  of  the  fringe  order  would 
vary  linearly  with  the  mass  flo;;  rate. 
Interpretation  of  data 

If  a  dark  polarizing  field  is  defined  as  one  in  which 
polarizer  and  arj^^lyzer  are  at  right  angles  (Ay  -  -5-  in  Table 
11)5  then  it  is  evident  that  two  different  dark  polarizing 
fields  were  used  in  the  final  testing.  Associate  those 
sets  in  which  either  the  polarizer  or  the  analyzer  is 
parallel  v;ith  the  principa.1  flow  direction  (y  =  0  or  y'  ~  '-') 
with  "porallel"  fields  and  associate  those  in  which  both 
analyzer  and  polarizer  are  oriented  across  the  flow  lines 
(^y  -  ±  ^^   y'  =  ^  j)   with  "crossed"  fields.   In  the  prelim- 
inary tests  the  dark  f  j.elds  were  of  the  parallel  type  cn].y. 

The  equation  describing  the  emergent  electric  field 
in  ti.-o-dimensional  flows  when  Ay  =  -j  v/as  given  in  Chapter 
Three  as : 
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Eq  =  Eo  sin(Tr'6An/A)  sin  2(y~<\r) , 
This  relationship  predicts  that  the  periodic  fringe 
patterns  will  be  the  same  for  all  dark  polarizing  fields, 
but  that  the  magnitude  of  the  emergent  field  will  decrease 
progressively  as  the  polarizer  and  analyzer  become  more 
and  more  nearly  parallel  with  the  optic  axes.  The  latter 
prediction  is  \7ell  supported  by  photographs  of  the  fringe 
patterns.  At  low  flow  rates,  when  all  theories  and 
measurements  predict  orientation  angles  near  tt/4  radians 
(see  Jerrard,  1959?  for  examples),  the  crossed  field  is 
dark  xirith  diffuse  fringes,  while  the  parallel  field  is 
light  with  broad,  lower-order  fringes  standing  out  in 
sharp  contrast  to  their  background.  At  higher  rates  of 
flow  a  lov/ered  orientation  angle  results  in  a  darkened 
parallel  field,  while  the  fringes  of  the  crossed  field  are 
sharp  and  clear. 

The  prediction  of  the  two-dimensional  theory  that  the 
phase  of  the  fringe  pattern  is  independent  of  the  orien- 
tation of  the  dark  field  is  not  supported  by  the  photo- 
graphic evidence.  The  fringe  pattern  in  the  crossed  field 
differs  by  half  a  fringe  order  from  its  parallel 
counterpart.  This  is  not  surprising,  since  the  orientation 
of  the  polarizing  field  cannot  be  neglected  in  the  three- 
dimensional  flow  which  is  actually  present. 

This  raises  tvro  questions,  one  theoretical  and  one 
practical.  The  theoretical  question  is:  "Can  equation 
(4,2),  v/hich  postulates  two-dimensionality  at  the  v/all, 


be  used  to  describe  a  three-dircerxSional  flow  field?"  To 
this  question  the  ansv;er,  based  upon  the  data,  is:  "I'es. 
provided  that  an  error  of  ±  1/2  fringe  is  permissible, 
since  that  seens  to  be  the  maximum  deviation  for  the  t-.-;o 
extreme  orientations."  The  practical  question  is:  "Given 
two  sets  of  "dark  field"  patterns,  which  should  be  used 
to  calculate  the  birefringence?"  The  ansv/er  is:  "VJhich- 
ever  can  be  read  with  the  greater  precision." 

With  these  considerations  in  mind,  an  error  of  1/2 
fringe  has  been  accepted  in  the  calculation  of  the  bire- 
fringence, and  the  crossed  field  has  been  chosen  to  provide 
the  data  since  the  only  fringes  which  are  not  clear  in  that 
field  are  those  of  the  first  and  second  order  which  are  too 
broad  for  accurate  dtermination  in  any  field.   In  circu- 
larly polarized  light,  of  course,  only  one  set  of  fringes 
is  observed  and  a  choice  of  patterns  is  avoided. 

To  the  1/2-fringe  bias  introduced  by  the  choice  of  a 
polarizing  field  must  be  added  the  experimental  deviation 
of  the  measurements  them.selves.  From  equation  (4.2): 

ki   =   [252(6i+52)Sop/5i]ci(I\/)/dG. 
Since  Sq  is  an  explicit  function  of  6i  and  62?  the  inde- 
pendent variables  are  01,  62?  P?  ana.  the  experimentally 
determined  slope  d{l\/)/dQ.    Since  the  four  variables  are 
independent,  it  follows  that  (e.g.  Beming,  1964,  p.  40) 
the  variance  in  k-j  is 
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where  cr^CX^)  designates  the  variance  of  the  s"*^^  independent 
variable  X^.   It  is  evident  that  the  relative  variance 
a^(X2)As  of  the  slope  defined  by 

ki'   =  dCN^/)/dG 
is  much  larger  than  the  relative  variance  of  the  other 
variables.  That  is, 

kl  '  Xg 

If  this  assumption  is  valid,  then 

r  oi^i 


(J2(ki)   = 


CT^Cki') 


.d(ki'). 

to  a  good  approximation. 

If  ki '  is  obtained  from  the  data  by  the  method  of 
least  squares,  then  according  to  linear  regression  theory 
(see,  for  example,  Volk,  1969)  the  estimated  variance  of 
ki '  is: 

s2(ki')  =  s^(V)/2'G% 

=   1^2' (x\2)  -  ki'S'GII^/]/(TTp-2)S'G2 
vfhere  s^ClI-;^)  is  the  estimated  variance  of  the  value  of 
I'k^j^   at  G  predicted  by   the  straight-line  approximation 
and: 

S'G^  =  EG^  -  (>:G)^/>lp, 

S'GN^.^2  ^    z;(GN^^2)  _  (SG)(>:N^/)/Np, 

Summations  are  for  the  Np  pairs  of  data  (N^.r^j  G-)  upon  which 
the  least  squares  analysis  was  performed. 
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The  estimated  variance  s^(ki ' )  is  based  upon  the 
dispersion  in  the  data  about  the  straight  lines  in  Figures 
13  to  16.  The  half- fringe  variation  introduced  by  the 
arbitrary  choice  of  a  polarizing  field  may  be  presumed  to 
have  a  systematic  bias  and  must  be  added  separately.  This 
may  be  done  by  increasing  the  variance  of  the  estimate 
s2(N  2)  by  (1/2)2.  There  results: 


v; 

s^(ki>)  = 


2n,  .^   -  Z'CNyr^)  -  ki'E'GTy^  4-  (l/2)2(Np^2) 


(N  -2)E'G2 


Neglecting  the  relative  variances  of  the  other  variables, 

s2(ki)  =  (ki/ki')  s2(ki'), 
as  assumed  earlier. 

The  estimated  variance  of  k2  can  be  obtained  imme- 
diately.  Since  ka  is  the  predicted  value  of  11^^^  v;hen  the 
flow  rate  is  zero,  the  estimated  variance  of  ka  is  the 
same  as  that  of  N.,^^  increased  by  the  systematic  bias  of 
(1/2)2: 

s2(k2)   =  ^H\^)   +  (1/2)- 

i:'(K,.,2)  ,  ICi'S'GIIy;^  +  (1/2)^(Nt3-2) 

Recalling  that 

1^2  =  kis  +  k2 
by  i::^pothesis,  and 

IT  =  An6/A , 
where,  in  this  case,  6  =  26i ,  the  variance  of  the  bire- 
fringence \-jill  be  proportional  to  the  variances  of  the 
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respective  coefficients  in  the  defining  equation: 
An2  =   (A2ki/46i2)  g  +   (A2k2/4-&i2) 

if  variations  in  A  and  6t  are  neglected. 
Computation  of  optical  coefficients 

The  optical  coefficients  ki ,  ka,  (A2ki/46i2), 
and  (A^kaA&i^)  have  been  defined  in  the  preceding 
paragraphs.  Values  for  these  constants  and  their  devia- 
tions are  provided  in  Tables  XVII  and  XVIII,  and  the 
corresponding  birefringences  are  plotted  in  Figures 
17  and  18,  Figure  1?  also  includes  an  extrapolation  of 
some  of  the  data  of  Peebles,  Prados,  and  Honeycutt  (1965)» 

As  indicated  in  Appendix  B,  two  stock  solutions  of 
Killing  Yellow  were  prepared.  The  data  in  Table  XVII  and 
Figure  17  are  for  the  original  stock  solution  which  was 
prepared  in  April  1972.  The  concentrations  cited  in  the 
plot  were  obtained  by  drying  samples  and  vreighing  the 
residue  as  described  in  Appendix  L.  The  effect  of  temper- 
ature upon  the  data  is  particularly  evident  in  Figure  17 
where  testing  occurred  over  a  short  enough  time  period  so 
that  the  concentration  of  the  solution  did  not  change 
greatly.  It  is  presumed  that  the  higher  birefringence  of 
the  less  concentrated  solution  of  Peebles,  Prados,  and 
Honeycutt  is  due  to  a  different  proportion  of  inorganic 
salts  in  the  commercial  dye. 

The  data  in  Table  XVIII  and  Figure  I8  are  for  the 
fresh  stock  solution  prepared  in  March  1973.  The  tests 
were  conducted  soon  after  the  solution  was  prepared,  and 
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day  by  day  changes  in  the  properties  of  the  medium  can 
still  be  observed.  The  birefringence  plotted  in  Figure 
18  shows  simultaneous  increases  from  run  to  run  due  to 
drops  in  temperature  and  to  changes  in  the  composition 
of  the  medium.  For  a  further  discussion  of  the  stability 
of  the  solution  and  changes  ^^rith  time,  see  Appendix  H, 
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Orientation  Angle 
In  this  dissertation  the  isoclinic  fringe  pattern 
in  two-dimensional  flow  fields  has  been  defined  in  terms 
of  the  orientation  angle  >];,  which  is  the  angle  between 
the  principal  flow  so^ls   and  ni  ,  the  "fast"  optical  axis. 
The  orientation  angle,  rather  than  the  more  familiar 
extinction  angle  X?  was  chosen  because  the  latter  variable 
is  defined  only  in  t'.v'o-dimensional  fields  when  the  polar- 
izer and  analyser  are  at  right  angles.  The  relationship 
between  the  orientation  angle  and  the  extinction  angle 
has  been  given  previously  at  the  end  of  the  section  on 
t^'/o-dimensional  flows  in  Chapter  Three,  Recall  that  when 
X  is  defined  it  has  the  value: 

r  I  Y  -  ^i^  I  ;y  -  ^y|  ^  f 

L  T  -  I  Y  -  ^1  |Y  -  Yi  >  4 

VJhen  -^   is  known,  X  can  be  determined  immediately  by  this 
relationship. 
Assumption  of  Form 

The  orientation  angle  (ertinction  angle)  of  Hilling 
Yellow  at  various  sliear  rates  was  not  measured  in  Lha 
present  dissertation,  Peebles,  Prados,  and  Honeycutt 
(1965)  found  that  the  extinction  angle  decreased  from 
45°  at  zero  shear  rates  to  an  asymptotic  value  in  the 
vicinity  of  20°  at  higher  shear  rates.  Although  it  is 
probably  coincidental,  Taylor  (1969)  obtained  19.5''^  as 
the  maximiuii  inclination  of  a  long  rod  falling  in  a 
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viscous  fluid.  There  nay  exist  sorne  model  ;;hich  relates 
these  two  numbers  since  the  hydrodynamics  of  flows  around 
rods  and  ellipsoids  is  a  common  basis  for  theories  of  bire- 
fringence, but  no  such  claim  is  m.ade  here.  Rather,  the 
value  of  19*5°  is  adopted  in  the  absence  of  a  better 
claim  by  some  other  value  in  the  sajne  experimental  range. 

Although  Peebles,  Prados,  and  Honey cutt  did  not 
choose  an  empirical  form  for  their  data,  it  appears  that 
the  equations 

tan  2%     =  tan  2/oo  +  ^Tg"^ 
or 

tan  2^  =  tan  2;]/co  +  l^^S"^ 
describe  then  reasonably  well  as  shoim  in  Figure  19.  This 
form  was  derived  by  Tsvetkov  (1951)  for  the  extinction 
aP-.gle  of  deformable,  anisotropic,  molecular  ellipsoids 
which  entrap  their  solvent o  Acceptance  of  the  Tsvetkov 
form  is  based  upon  its  convenience  and  not  the  model: 
hov/ever,  the  Tsvetkov  model  is  consistent  \;ith  the  possi- 
bility mentioned  in  Chapter  Tv/o  that  Milling  Yellox;  has  a 
lyo tropic  mesoDhaso. 
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FIGURE  19,  Extinction  angles  measured  by  Peebles,  Prados, 
and  Honeycntt  (1965)  replotted  to  obtain  straight  lines. 


CIL^PTSR  FIVE 
DETERl-ilKATIOrT  OF  RKSOLOGICAL  PROPERTIES 

This  chapter  discusses  the  method  by  which  the  rheo- 
logical  properties  of  Milling  YelloX'/  v;ere  measured.  The 
chapter  is  divided  into  five  parts.  The  first  discusses 
empiricism  in  rheological  formulations  and  concludes  v/ith 
an  equation  using  tvjo  constants  \vhich  sa.tisfies  data 
reported  by  Pebbles,  Prados,  and  Ploneycutt  (1965)  ija  thu?ee- 
constant  form.  The  second  part  of  the  chapter  analyzes 
the  Hoppler  Rheo-Viscometer,  an  instrument  generally 
considered  unsuited  for  the  measurement  of  basic  rheo- 
logical constants.  The  third  part  tabulates  the  data 
obtained  v;ith  the  Rheo-Viscometer.   The  fourth  part  uti- 
lizes these  data  to  construct  a  functional  relationship 
between  shear  stress  and  shear  rate  in  Killing  Yellow. 
The   fifth  and  final  part  considers  the  range  of  shear 
stresses  for  which  the  constructed  relationship  is  valid. 
Empiricism  in  Rheological  Formiilations 

The  late  Professor  Isadore  Aindur  once  enlivened  a 
particularly  abstruse  derivation  of  a  thermodynamic 
relationshiT)  with  the  comr.ent:* 


*I.  Amdur,  Physical  Chemistry  Lectures,  Massachusetts 
institute  of  Teclmology,  CambrD.dge,  Massachusetts,  SDring  j 
1955 
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Gentlemen,  you  expect  science  to  be  simple. 
You  are  mistaken.  It  is  not  science,  but  scientists, 
who  are  simple. 

The  Theological  behavior  of  Milling  Yellow  is  not 
simple.  As  the  review  in  C'fapter  Two  indicates,  neither 
continuum  mechanics  nor  the  methods  of  analytic  model 
construction  are  able  to  provide  useful  constitutive 
relationships  which  can  be  applied  to  the  solution  of 
practical  problems.   In  their  absence  aji  empirical 
approach  becomes  necessary,  but  such  an  approach  must  be 
provided  with  sufficient  justification  not  to  fly  in  the 
face  of  established  principles.  Perhaps  a  reasonable 
compromise  is  represented  by  the  so-called  Stokesian 
fluid*  The  discussion  which  follows  is  taken  in  large 
part  from  Frederickson  (1964-,  2E«  68-89). 

Consider  an  incompressible  Stokesian  fluid  as  defined 
by  the  constitutive  relationship 

where  the  stress  tensor  may  be  represented  as  a  "poly- 
nomial" in  the  rate-of-deformation  tensor  g^-n.  If  the 
fluid  is  isotropic,  the  polynomial  expansion  reduces  to 

^±i     =     -P  ^ij  ^  cli  gj.  +  d2  gi:-gki 
where  p  is  the  pressure,  6^j  is  the  Kronecker  delta,  and 
di  and  d2  are  polynomials  of  the  three  principal  invar- 
iants, I,  H,  and  m,  of  the  rate-of-deformation  tensor. 
For  flows  which  are  steady  and  axial, 
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Ux      =  0,  Ua  =   0,  U3   =  ■U3(xi,:c2) 
the  tensor  invariants  become  simply 

I  =  m  =  0, 

n  =  -.[(ou3/dxi)2  +  (dU3/dx2)^]A. 
Thus 

which  is  reniniscent  of  Sricksen' s  (1960a): 

'^ij  "  -P^j  -  ^;gij  -  i^Qn-^P)Si_kgkj  -^  2^n^Si/St 
where  /o^  is  the  intrinsic  viscosit;^'"  and  0^^  and  p  are  the 
primary  and  secondary  normal  stress  functions.  The 
stress  functions  have  been  calculated  by,  among  others, 
Giesekus  (I962)  and  measured  by  Tanner  (1970) •  The 
magnitude  of  9^^  is  snail  in  most  fluids  and  S  is  smaller 
leading  Bird*  to  conment  humorously  that  integrated  over 
the  last  decade  of  investigation,  p  is  zero.  The  presence 
of  the  Jaumann  derivative 

/"^cr  .  .        Act  .  .  Ac-  .  . 

reflects  local  anisotropy  in  the  medium.  The  tensor  Wj_j 
is  the  vorticity  or  spin  tensor. 

In  the  present  stiidy  the  effect  of  normal  stress 
fionctions  upon  the  rheological  behavior  of  Killing  Yellow 
is  neglected  leading  to  the  basic  relationship  for  a 
"generalized  Ilewtonian  fluid"  with  the  simple  form: 


*R.  Byron  Bird,  Rheology  Lectures,  University  of  Florida, 
Gainesville,  AxDril  1973 
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The  sacrifice  in  dropping  the  terns  with  6^  and  p  is 
less  than  it  appears.  Ericksen  (1956)  has  shovm  th-at  if 
the  surfaces  of  constant  velocity  for  a  flow  are  neither 
planes  nor  circular  cylinders,  the  only  Stokesian  fluid 
vrfiich  can  be  descrihed  by  the  conditions  for  steady  axial 
flow, 

Ut    =   0,   U2   =   0,   U3   =   U(X1,X2), 

is  one  in  which  9^  and  p  are  either  directly  proportional 
to  nj   or  zero. 

For  steady  axial  flow,  the  components  of  t-j_,-  for  a 
generalized  Uewtonian  fluid  are 

T-,  1    =  122       =       1^3  3   ~   -Pj 
T^l  2   =   '^21    =   Oj 

Ti  3  =  T31   =  -/^(5u/dxi  )/2, 

T23   =       1^32   -   -/«(&u/9X2)/2. 

Squaring  the  expressions  for  0:13  and  1:23  and  adding: 
T^^^   +  T23^   =  /)72[(d-u/9xi)^  -  (5u/5X2)^]/4. 
The  scalars  t  and  g  are  defined  by  comparing  the  terins  in 
this  expression  with  the  corresponding  ternis  in 

Recognizing  that  the  intrinsic  viscosity  was  assumed 
to  be  a  function  of  H  and  that  gV4-  =  -ZE,  the  final 
general  form  for  the  empirical  relationship  is  just 

or  with  equal  generality, 
g  =  gC'^). 
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In  thair  analysis  of  viscosity  rneasurements  made  in 
a  capillary  viscometer,  Peebles,  Prados,  and  Honeycutt 
(1965)  found  that  Killing  Yellow  is  well-represented, 
though  not  uniquely,  by  the  Powell-Eyring  equation: 

^  =  >ioog  +  Ap(;io->ioo)  sinh""Ug/Ap). 
The  flexibility  of  empirical  representations  is  indicated 
by  an  alternate  formulation  of  the  same  data  for  shear 
rates  belo\j  4-00  sec"''.  Peebles,  Prados,  and  Koneycutt 
provide  a  plot  of  apparent  viscosity  versus  wall  shear 
stress.  If  these  same  data  are  replotted  as  )i^  versus 
1/t,^^2j  as  in  Figure  20,  a  straight  line  is  obtained; 

For  flov;  in  capillary  tubes  the  classical  approach 
of  Rabinowitsch  (1929)  and  llooney  (1931)?  as  modified 
by  Krieger  and  Maron  (1952),  leads  to  the  relationship 
betv7een  shear  rate  and  shear  stress  at  the  capillary 
wall: 


where 


and 


gy-r,,  =   a/ji^)  +  ^  d(i/ji3)/dT,, 


c,,^  =  H'.p/2L^ 


yi^     =      ttaprVSL^Q. 
The  variables  in  these  expressions  are  R,  the  capillary 
radius,  Ap/L^,  the  pressure  drop  along  capillary  length  L^, 
and  Q,  the  volumetric  flow  rate  through  the  capillary-. 

Ifiien  the  equation  for  the  straight  line  in  Figure  20 
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(cms/dj/ne) 


FIGUEE  20,   Data  of  Peebles,  Prados,  and  Honeyc-utt  (19^5) 
replotted  to  obtain  a  linear  relationship. 
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is  substituted  into  Krieger  and  Karon' s  equation,  the 
result  is^ 

}ico-^,;^  +  (3k/2) 
g(^)   =   ^3  na,^^  ^   (3k/2)  . 

This  relationship  differs  qualitatively  from  the 
Povell-Eyring  equation  in  being  of  the  form  gCx)  rather 
than  T(g),  Further,  it  is  expressed  in  terms  of  tiro 
constants,  yLao    and  k,  rather  than  the  three  constants 
used  in  the  Powell-Ejrring  equation.  Inspection  of  the 
straight-line  fit  in  Figure  20  also  suggests  that  the 
expression  for  gi'z)   given  above  may  not  be  valid  for 
values  of  t  above  about  500  dyne/cm^.  V/ere  it  not  for 
this  restriction,  at  high  values  of  t  the  above  equation. 
v7ould  reduce  to: 

gCi:)  =  t/jlioo  (t  large). 

At  very  small  values  of  t  it  is  evident  that 

g(T)  =  3T3/2k  (t  small), 

but  this  represents  an  extrapolation  of  the  experimental 
data  of  Peebles,  Prados,  and  Honeycutt  much  belo^y  the 
lowest  shear  stress  reported. 

Besides  these  relationships,  a  largo  nimber  of 
alternate  empirical  forms  are  available  for  consideration. 
Skelland  (I967,  pjo.  ^-14)  lists  four  equations  for  fluids 
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\rlth   yield  stresses,*  seven  for  fluids  without  yield 
stresses,  and  four  for  fluids  flovjirig  in  tubes  where  the 
stress  at  the  wall  is  related  to  the  volujnetric  flow  rate. 
Although  the  form  of  nany  of  these  equations  is  influenced 
by  nolecular  theory  or  other  ar^alysis,  the  spirit  is 
empirical,  and  the  constants  are  determined  by  experiment. 
Analysis  of  Sliding  Ball  Viscometer 
Since  standardized  samples  of  Killing  Yellow  cannot 
be  prepared  vrith  predetermined  properties,  and  because 
the  properties  of  the  solution  change  v/ith  time,  the 
constants  which  appear  in  the  rheological  relationships 
describing  Killing  Yellow  must  be  determined  hj   exper- 
iment.  In  previous  tests  (e.g.  Peebles,  Prados,  and 
Honeycutt,  19^5)  the  instrument  used  has  been  a  capillary 
viscometer.  In  the  present  instance  the  instruiaent  used 
was  the  sliding  ball  viscometer,  sometimes  called  a 
consistometer  (Blair,  194-9)  p.  146).  The  sliding  ball 
viscometer  has  not  been  considered  an  appropriate 
instriunent  for  the  determination  of  basic  rheological 
properties.  For  this  reason  an  analysis  h.as  been  carried 
out  by  vrfiich  the  rheological  constants  may  be  educed  from 
measurements  made  v/ith  the  sliding  ball  viscometer.  The 
pages  which  follow  describe  this  analysis. 


*in  fluids  with  yield  stresses,  t(0)  t^   0.  Although  Milling 
Yellov/  has  a  yield  stress  (Thurstcne  and  Schrag,  1962),  it 
is  negligible  in  most  applications. 


Ill 


Related  Investigations 

The  sliding  ball  viscometer  consists  of  a  cylindrical 
container  inclined  at  a  slight  angle  from  the  vertical 
within  which  a  closely  fitted  sphere  is  driven  dovmv;ard 
through  the  liquid  by  a  known  force  transmitted  along  a 
rigid  shaft.  The  present  design  of  the  Hoppler  Rheo- 
Viscometer  v;hich  is  shovm  in  Figure  21  differs  little 
from  the  original  (Hoppler,  194-la).  The  variables  of 
interest  are  the  force  applied  to  the  ball  and  the  time 
required  for  the  ball  to  move  downward  through  a  set 
distance.  Advantages  of  the  instrument  include  small 
sample  sizes  (less  than  20  ml),  easy  control  of  temper- 
ature and  evaporation,  speed  and  simplicity  of  operation, 
excellent  reproducibility,  and  the  traditional  advantage 
of  such  instruTients:  the  duration  of  shear  is  short.  The 
principal  disadvantage  is  also  the  traditional  one:  neither 
the  shear  stress  nor  the  shear  rate  is  constant  within  the 
instrument.  Blair  (1949)  has  discussed  these  advantages 
and  disadvantages  without  reaching  a  rigid  conclusion* 

The  sliding  ball  viscometer  v;as  developed  in  Germany 
at  the  outset  of  Wox^ld  War  II.  Hoppler 's  (1941a,  194rD, 
1942)  descriptions  were  either  qualitative,  or  the  shear 
stress  in  the  instrument  was  treated  as  constant  and  equal 
to  the  force  transmitted  by  the  shaft  divided  by  the  cross- 
sectional  area  of  the  ball,.  Follov;ing  the  war,  analyses  of 
the  sliding  ball  instru-ment  are  rare  (e.g.  Lindgren,  1957) 
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FIGURE  21.  Hoppler  Rheo-Viscometer.  ITumbered  ite'ns: 

(1)  scale  beam  transmitting  load  to  siiaft  and  ball; 

(2)  shaft  and  sliding  ball;  (3)  counterweight  to  maintain 
offset  position  of  ball;  (4)  fluid  sample;  (5)  calibrated 
mass  determining  load  on  ball;  (6)  imused  masses;  (7)  dial 
indicating  position  of  ball  in  cylinder  to  nearest  0,1  mm; 
(8)  slider  to  compensate  for  buoyancy  effects  arising  from 
samples  of  varying  density;  (9)  eccentric  >:n.ob  which  is 
used  to  raise  and  release  beam;  (10)  water  from  constant 
tem.pera-ture  circulator;  (11)  precision  thermometer; 

(12)  receptacle  for  next  sample. 
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and  to  the  present  day  the  instriiment  manual  furnished 
by  the  manufacturer  speaks  only  of  the  "shear  stress" 
defined  by  Koppler's  original  papers. 

If  the  sliding  ball  viscometer  is  absent  from  the 
literature,  the  same  cannot  be  said  for  its  cousins: 
the  rolling  ball  viscometer,  the  ball  and  tube  flovrmeter, 
and  the  falling  cylinder  viscometer.  Ewing  (1925)  set 
some  sort  of  record  for  ex  post  facto  reminiscences  by 
describing  his  experiments  with  the  ball  and  tube  flow 
meter  in  I876.  He  points  out  that  the  ball's  position 
in  the  tube  is  not  stable  when  the  tube  is  vertical  and 
provides  photographs  of  colored  dye  flo^-zing  up  the  center 
of  the  tube,  around  the  ball,  and  up  the  wall  while  the 
inclination  of  the  tube  maintains  the  ball  off-center  in 
an  eccentric  position.  Happier ^s  (1933)  measurements  with 
a  rolling  ball  viscometer  demonstrated  the  dependence  of 
fall  times  upon  the  ratio  of  ball  radius  to  cylinder  radius 
and  showed  that  accurate  viscosity  measurements  could  be 
made  with  a  variety  of  materials  following  calibration 
v;ith  water.  Eoppler's  photograph  of  dye  flowing  around 
the  ball  is,  unlike  Ewing' Sj  symmetric  above  and  below  the 
ball. 

Hubbard  and  Brovm  (194-3a)  derived  a  general  equation 
for  the  rolling  call  viscometer.  A  calibration  constant 
was  necessary  for  this  equation  and  was  found  to  vary  by 
about  10  percent  when  used  with  fluids  ranging  from  air  to 
ben.sene  (1943b).  An  analytic  form  for  this  constant  was 
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provided  by  Lewis  (1953) •  The  corresponding  analyses  of 
the  falling  cylinder  viscometer  are  those  of  Lawaczeck 
(1919)  and  Lohrenz,  Smith,  and  Kurata  (i960).  Christoph- 
er son.  and  Dowson  (1959)  have  considered  the  falling, 
rotating  ball.  Bauer  and  DuPuis  (I967)  discuss  drag  on  a 
sphere  rolling  in  a  closed  tube.  Floberg  (1968)  is  the 
author  of  a  treatise  on  the  ball  and  tube  flo\^meter  and 
the  rolling  ball  viscometer.  All  of  these  studies  have 
been  concerned  with  Newtonian  fluids. 

The  extension  of  this  work  to  non-Newtonian  fluids 
requires  an  accurate  description  of  the  flow  in  the  visco- 
meter annulus.  For  a  close  fitting  sphere  or  cylinder 
the  usual  first  assumption  h^s  been  to  assume  that  the 
flow  pattern  is  essentially  that  of  two-dimensional  flow 
along  a  channel  having  as  its  boundaries  two  circular 
surfaces. 

Except  for  the  elementary  case  of  Newtonian  flow  in. 
a  concentric  annulus  (e.g.  Lea  and  Tadros,  1931)?  the  first 
solution  of  this  problem  appears  to  be  that  of  Stevenson 
(1948),  v;ho  solved  the  analogous  problem  of  torsion  in 
hollow  tubes.  The  analogj^  between  steady  axial  flows  and 
torsion  in  cylindrical  bars  h^as  been  discussed  in  Chapter 
Two.  For  unit  eccentricity,  v;hen  the  ball  or  cylinder 
slides  along  the  cylindrical  wall,  the  bounding  surface  is 
simply  connected  satisfying  the  necessary  condition  for  the 
mathematical  solutions  of  the  two  problems  to  be  the  same. 
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Ashare,  Bird,  and  Lescarboura  (1965)  considered  the 
flow  of  non-Newtonian  fluids  in  a  concentric  annulus  and 
found  that  non-Newtonian  flow  could  be  approximated  by- 
slit  flow  if  one  is  interested  only  in  the  average  fluid 
velocity  in  the  channel.  Chen,  Lescarbouraj  and  Swift 
C1968)  studied  the  effect  of  eccentricity  upon  the  flow, 
obtaining  an  exact  expression  for  tlie  velocity  of  a 
Newtonian  fluid  in  an  eccentric  annulus.  An  extension  of 
their  method  to  non-Nev;tonian  fluids  was  also  indicated. 
Both  Ashare^s  group  and  Chen's  applied  theii*  results  to 
the  falling  cylinder  viscometer.  They  assumed  that  the 
cylinder  was  long  enough  so  that  entrance  and  exit  effects 
could  be  neglected. 

The  sliding  ball  viscometer,  despite  a  superficial 
resemblance  to  the  rolling  ball  instriment,  resembles  the 
falling  cylinder  device  when  the  flow  is  analyzed.  The 
similarity  lies  in  the  fact  that  neither  the  sliding  ball 
nor  the  cylinder  rotate  as  they  move  dox^mward.  The 
difference  in  the  two  falling  bodies  is  more  apparent  than 
real.  At  the  close  tolerances  for  which  the  sliding  ball 
viscometer  is  used  in  the  present  dissertation,  it  may  be 
sho^iTti  that  the  resistance  of  the  ball  occurs  over  an 
"effective  length"  which  is  less  than  10  percent  of  the 
ball's  diameter.  Over  this  effective  length,  the  radius 
of  the  "cylinder"  which  the  ball  approximates  changes  by 
less  than  1  percent. 
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Analys  is 

Paralleling  the  derivation  for  constitutive  rela- 
tionships in  a  capillary  as  carried  o-at   by  Rabinowitsch 
(1929)5  Mooney  (1931),  and  Krieger  and  Maron  (1952),  it 
is  first  necessary  to  define  the  volumetric  flow  of  the 
liquid  upward  in  the  cylinder  as  the  ball  descends.  The 
ball  moves  so  slowly 
compared  with  the  mean 
speed  of  the  liquid 
(less  than  1  percent) 
that  the  speed  of  the 
ball  can  be  neglected 
and  the  flo'v;  becomes 
that  of  a  non-Nev;tonian 


FIGURE  22. 
Geometry  in  annulus. 


liquid  in  an  eccentric  annulus.  Using  the  parameters 
sho\m  in  Figure  22,  one  obtains: 

Q  =   j'^u  dA  =  2  f  [  j'  u(a+C)dC]  dG 

where  Q  is  the  volumetric  flov;  rate,  u  is  the  axial  speed 
at  coordinates  (CjG),  and  6  is  the  width  of  the  annulus • 
The  m.agnitude  of  6  is 

b     =     &^(l-cos  e)/2  ±  5jjj2(l-cos  2e)/l6a  +  Oib^J/a^) 
where  the  sign  of  the  second  and  higher  order  terms 
depends  upon  whether  5  is  measured  along  an  outrv'ard  normal 
from  the  ball  or  an  inward  normal  from  the  cylinder.  Since 
6rn/s.  =  0,0088  for  the  Hoppler  Rheo-Viscometer  used  in  the 
present  study,  no  serious  error  is  introduced  if  one 
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neglects  the  terms  in  b^,^   and  higher  to  obtain 
6  =  6j^(l-cos  e)/2. 
Consider  the  term 
I  =  5>(a+C)dJ: 
in  the  expression  for  Q  just  given.   Set  C  =  '7o  +  7 

and  obtain 

I     =     fuCa  +  I  +  />7o)   d/770  +    5^^(a  +  I  +  ^0)   d/yo- 
In  a  narrow  annulus  ii(-/>/o)    =  ■u(/77o)5    so 

I     =     2fS(a  +1+^0)   d/r^o. 
Integrating  by  parts, 


2u 


(a  H-  |),>^o  +  ■I'Yo^' 


u  =  0,  />/o   -  6/2 


-  2  il 


2'"( 


du. 


Recognizing  that  the  first  term  is  zero  at  both  limits, 
and  replacing  'T^o  by  C  -  2?  obtain 


u. 


I  =  2  1  (a  +  •|)(f  -  O   du. 


r    <  ^ 


Substituting  into  the  original  integral  for  Q  yields 
Q  =  Aj   (a  4-  |)(-|  -  ^)  du  d9. 

o  o 

The  unl-movm  constitutive  relationship  vrfiich  is  sougW 
for  in  the  viscocetric  measurer.ents  is 

du/dC  =     gCi:). 
Substituting, 

Q  =  4  f  f(a  +  ^)(|  .  C)  g(T)  dC  de. 
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At  this  point  it  is  necessary  to  drav;  a  force  balance 
on  an  annular  element  of  the  fluid  subtending  an  arc  d9, 
having  a  thickness  dC  and  a  length  z.   It  is  assumed  that 
the  change  in  the  shearing  stress  in  the  angular  direction 
dG  is  negligible  compared  with  the  change  in  the  direction 
dCj  and  that  the  flov;  is  steady  and  axial  so  that  the 
shearing  stress  is  unchanged  throughout  the  length  z. 
Under  these  assumptions  the  force  balance  can  be  ^inritten 
as 

p(a  +  C)  de  dC  =  -  dT  (a  +  C)  dG  z 
which  reduces  to 

-dr/dC  =  Ap/z. 
The  change  in  pressure  from  one  side  of  the  ball  to  the 
other  is  due,  in  the  Hoppler  Rheo-Viscometer ^  to  the  force 
applied  to  the  ball.  Hence,  identify  the  pressure  drop  Ap 
along  the  annular  element  vrith  the  force  applied  to  the 
shaft  divided  by  the  cross-sectional  area  of  the  ball. 
This  is  the  quantity  which  the  instrument  manual  designates 
as  P.  the  "average  shear  stress"  in  the  fluid.*  There  is 
no  analytical  expression  for  the  "effective  length"  L 
which  replaces  z: 

Ap/z  =  P/L. 
The  effective  length  \j111   be  obtained  experimentally. 

Since  P  and  L  are  constants,  the  force  baHance  can- 
be  integrated  from  C  =  0  to  C  =  &/2  to  obtain 


*The  units  of  P.  according  to  the  manufactuj-er ,  are  gn/cm.2 
A  conversion  to  dyne/cm^   is  required  before  the  numerical 
values  given  for  P  can  be  used  consistently. 
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T   -   P(-|  -  C)/L. 
At  the  wall, 

T^,^  =  P6/2L  =  P6jj,(l-cos  8)/4L. 

Substituting  the  value  for  the  shear  stress  into  the 
integral  for  Q  yields 

Q  =  ^   (a  +  I)   Tg(T)  dT  dG  . 

P'=^   '^o  Jo 

It  is  a  simple  matter  to  rex^rrite  Q  in  terms  of  the  "fall 
time"  t  for  the  ball.  Applying  the  continuity  requirement, 

Q  =   ^^^annulus  =  (UA)^^!!  =  |  ^raS 
where  ,f  is  the  distance  (normally  3  cm)  through  which  the 
ball  is  driven  in  time  t.  Substitution  and  rearrangement 
leads  to: 

In  the  analyses  v;hich  follow  it  will  be  noted  that  -|  has 
usually  been  neglected  compared  with  a  in  order  to 
reduce  the  complexity  of  the  integration  procedures. 
Although  justified  by  the  relative  magnitudes  of  the 
quantities  (-  =  0.0044a),  this  assum.ption  has  been  made 
for  convenience  rather  than  necessity.  Typically  a 
rigorous  integration  has  been  possible  even  when  it  was 
not  performed. 

For  a  Ilev/tonian  fluid  the  response  of  the  viscometer 
can  be  obtained  directly.  Substituting  gC-r)  =  -z/^jl,^   into 
equation  (5.1)  J^ields  (when  the  6/2  term  is  neglected): 
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1/t  =  — ~ —  \  1  t2  dT  de? 

Substituting  the  value  for  the  stress  at  the  \-Jall   given 
on  the  previous  page, 

1/t  =  3-^^:^4-g^  (1-cos  e)3  de, 

=  5P6ja3/96;ijja^L. 
Hence,  for  a  Newtonian  fluid, 
?l^     =   55ia3pt/96aJ^L. 

Calibration  of  the  viscometer  yields  a  numerical 
value  for  the  ratio  5b^^/96a.JlL   frOK  vrhich  the  effective 
length  may  be  calculated.  When  this  is  done  for  the 
Hoppler  instrument  the  effective  length  is  fcujid  to  be 
0,1365  cm,  compared  with  a  ball  diameter  of  1,596  cm. 
Details  of  this  calibration  may  be  found  in  Appendix  E, 
Since  the  credibility  of  the  present  analysis  depends 
upon  the  concept  of  a  fixed  effective  length,  it  appears 
desirable  to  give  some  attention  to  the  experimental  value 
for  L  and  the  inferences  to  be  drax^m  from  its  magnitude. 

In  a.ssuming  that  the  flow  past  a  sliding  ball  Gs,a  be 
modeled  on  the  flovr  in  an  eccentric  aimulus  of  the  same 
cross- section  it  is  implicit  that  the  ball  be  replaced  by 
a  cylinder  of  length  L.^  Three  questions  naturally  arise: 


*VJhen  this  a:ialysis  is  applied  to  a  falling  cylinder  visco- 
meter the  length  of  the  actual  cylinder  should  approxiiiiate  L. 
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1.  —  Why  use  a  ball  instead  of  a  cylinder? 

2.  —  Hov;  valid  is  the  assumption  that  the  ball  can 
be  modeled  by  a  cylinder? 

3.  —  Hov;  constant  is  the  effective  length  from  fluid 
to  fluid  and  from  one  force  on  the  ball  to  another? 
These  questions  will  be  answered  individually. 

1.  Vfny  use  a  ball?  -  All  of  the  analysis  which  is 
provided  in  this  chapter  models  flow  past  a  ball  at  the 
point  of  minimum  clearance  on  flow  past  a  cylinder  of 
the  same  dimensions;  hence,  the  analysis  should  be  valid 
without  significant  change  for  a  sliding  cylinder  visco- 
meter, TTaturally  the  primary  ii^petus  for  the  analysis 
was  the  need  to  analyze  the  instrument  which  v;a3  actually 
used.   It  is  valid  to  inquire,  nevertheless,  v/hether  an 
instrument  with  a  fitted  cylinder  v/ould  have  been  prefer- 
able to  the  present  instrument.   It  appears  lilcely  that 
the  sliding  ball  possesses  the  advantage  that  a  smaller 
force  is  required  to  overcome  the  shearing  stresses  in  the 
eccentric  annulus  than  would  be  the  case  with  a  cylinder 
of  sufficient  length  that  entrance  and  exit  effects  could 
be  neglected.  On  the  other  hand,  the  effective  length  of 
the  sliding  cylinder,  even  if  a  correction  for  end  effects 
is  necessary,  corresponds  more  closely  to  an  actual 
physical  dimension  than  the  effective  length  of  the  slidrcig 
ball.  VJhich  instrument  would  be  preferable  in  a  general 
case  would  depend  upon  the  rela.tive  importance  of  these 
considerations . 


122 


2.  How  v/ell  does  the  ball  model  a  cylinder?  —  In 
terns  of  changes  in  the  radius  of  the  cylinder  which  the 
ball  models,  the  assumption  that  the  ball  is  a  cylinder 
over  the  effective  length  L  =  O.1365"  cm  is  a  good  one.      If 
it  is  assumed  that  the  position  of  the  effective  length  is 
such  that  it  extends  equally  above  and  below  the  point  of 
minimurn  clearance,  then  the  change  in  the  radius  of  the 
ball  is  about  0,72  percent.  Despite  the  instability  of 
flow  in  divergent  chan-nels,  as  calculated  by  liillsaps  and 
Pohlhausen  (1953)  extending  the  work  of  Jeffery  (1915)5 
Hamel  (19l6)  and  others,  it  appears  tlmt  at  the  Reynolds 
nusibers  which  are  present  in  the  chamiel  (see  Appendix  J") 
a  bad:  flow  does  not  occur  within  the  computed  effective 
length.  The  width  of  the  eccentric  annulus  does  change 
significantly  over  the  effective  length.  Above  and  below 
the  ball,  at  points  corresponding  to  the  limits  of  the 
effective  length,  the  channel  has  a  width  twice  that  at  the 
point  of  minimum,  clearance.  The  angle  which  the  tangent  to 
the  ball  nal'es  with  the  cylinder  wall  at  these  points  is 
about  5  degrees. 

An  alternate  miodel  has  been  considered  by  Sestak  and 
Ambros  (1973)  '•'ho  assum.ed  tliat  the  flow  remains  parallel 
v;ith  the  cylinder  walls,  but  that  the  flov;  rate  decreases 
in  proportion  to  the  increase  in  cross-sectiona.1  carea.. 
Integrating  along  the  boundary  of  the  sphere  in  the 
direction  of  flow,  they  obtained  viscometer  responses  for 
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llev/tonian  and  povrer-law  substances.  Their  results  for 
the  latter  extended  a  similar  analysis  by  Bird  and 
Turian  (1964).  V/hile  the  model  of  Sestai^  and  Ambros 
overcomes  one  assimption  of  the  present  model  —  that 
the  flov;  rate  is  unchanged  along  the  effective  length 
of  the  annulus  —  it  retains  another:  that  the  radial 
components  of  the  velocity  are  negligible.  VJhile  this 
assumption  appears  reasonable  for  integrations  orthog- 
onal to  the  flow  path  at  the  point  of  minimum  clearance, 
as  in  the  present  dissertation,  it  appears  less  reason- 
able \'fhQii   the  integration  occurs,  as  it  does  in  the 
Sestalc  and  Ambros  analysis,  in  the  direction  of  flov/. 
Since  the  assumption  made  in  the  present  dissertation 
permits  the  determination  of  the  rheological  properties 
vrhen  the  form  of  the  constitutive  relationship  is  not 
Imovm,  a  determination  which  is  not  possible  when  the 
model  of  Sestak  and  Ambros  is  employed,  it  would  appear 
that  the  advantages  of  modeling  a  cylinder  by   means  of  a 
ball  overweigh  the  corresponding  objections* 

3-     How  constant  is  tlie  effective  length?  —  The 
constancy  of  the  effective  length  as  va.rious  fluids  are 
placed  in  the  viscometer  is  a  m.atter  of  conject'ore.   It 
is  clear  from  the  calibration  in  Appendix:  E  tiiat  the 
change  is  negligible  for  xievrtonian  fluids  through  a  wide 
range  of  shea.r  rates  and  between  ile^rtonian  fluids  when  the 
viscosity  varies  by  more  than  a  factor  of  ten*  The 


assiomption  made  in  this  dissertation  is  that  the  change 
in  the  effective  length  is  also  negligible  for  non- 
Newtonian  substances. 
Integration  of  the  viscometer  equation 

xhe  viscometer  equation,  equation  (5.1)?  has  already 
been  integrated  for  a  liev/tcnian  fluid.  Another  important 
class  of  fluids  for  which  the  equation  may  be  integrated 
directly  are  the  power- law  fluids  vrhich  may  be  represented 

in  the  one- dimensional  form 

m 
T  =  K  g 

where  K  and  m  are  constants.  Substitution  into  the  visco- 
meter equation  yields 

1/t  =   (4-LV7Ta2Xp2K^/in)r  (a  +  |)[\^  dT  dG. 
Neglecting  6/2  compared  \ilth   a, 

1/t   =   (^LVrraJlP^K''^)]  It"'  dT  dG. 

o  o 

Integrating  once  and  setting  (l+2m)/rr;  -   M, 


1/t  =   (4LVTTa«p2>FJ'^)J  (P6^/4-L)^'(l-cos  9)^'  dG, 


\  2^^     sin2i-x?'  (2  dcp') 


where  29'  =0.   Integration  yields 

1/t  . ^Jl^^ pVm 

2^^'~  ^yJTaJL^"-  2l-^>r(ll+l) 

where  r(M+V2)  and  F(M+1)  are  gamma  functions.  Thus  the 
response  of  the  viscometer  to  a  power-law  substance  shows 
the  sajne  exponential  characteristics  as  in  one-dimensional 
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shear  flov;.  The  model  of  Sestak  and  Ambros  (1973)  yields 
a  siinilar  result. 

The  treatment  of  materials  with  yield  stresses  (such 
as  a  Bingham  plastic)  or  bilinear  responses  is  similar, 
but  a  separate  consideration  of  the  limits  of  integration 
is  required.  An  example  of  a  bilinear  material  having 
the  response 

t/)1q  t  <  Tc 


g  = 


^-^  -^^-^c 


is  treated  in  Appendix  F, 

If  the  form  of  the  constitutive  equation  is  unlcno^'/n 
or  resists  integration  in  explicit  form  it  is  possible  to 
obtain  a  plot  of  g(T)  versus  t  from  experimental  values 
of  P  and  t. 

Replace  t^^  in  equation  (5«l)  "by  P6/2L  to  obtain  (a>>5jjj) 

PVt  =   (^LV^aOrr*'"-^  g(T:)  dT  de.  (5.2) 

Differentiating  with  respect  to  P  yields 

d(PVt)/dP  =  (P/7Ta^)  i  52  g(5P/2L)  de. 
Expand  g(T)  as  a  power  series  in  t: 

g(T)   =   OCo  +  CCi  T  +  CC2  T  /2!  +  ,. 

At  the  wall,  where  t^.  =  oP/2L  =  P5rr,(l-cos  9)/4L, 

Si^j)      =  g(5P/2L)  =  £aj,6ciV(l-cos  e)''/4.'^''n! 

Substituted  into  the  expression  for  d(p2/t)/dP  with 
6jjj/4-L  set  equal  to  po}  there  results: 
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^2  ,,  X         .   2 


P   dP        4naA  Jo  ^o  ^^"^  ' 

which  is  integrated  to  yield: 

1  d(PVt)  ^   &m  ^  anPo^P^2^"^^  ^  1.3.5- » (2n+3) 
P   dP       4-TTaiJ^o     n!     2  2.4.6- .(2n+4-) 

The  quantity  on  the  left  hand  side  of  this  equation  is 
obtained  experimentally.  By  least  squares  analysis  or 
other  suitable  technique,  fit  a  polynomial  of  the  form 

1  d(P^/t)        J-   D  J.   -02  ,      •□n 
p  — ^p — ^  =  ao  +  arP  +  aaP''  +  ..  a^^P 

to  the  experimental  data.  Alternately  one  may  find  an 
empirical  fit  in  the  form 

1/t  =  bo  +  biP  +  b2P2  +  ..  b^P^, 

in  which  case 

2-n  ~  (n+2)t>n. 
Once  the  coefficients  have  been  obtained  for  the  exper- 
imental data,  a  term  by  term  comparison  yields: 

aii^  ^a-;^  2-4-6--(2n+3) 
^^    6ni^l3o^2n  1.3-5--(2n+4)  * 

Substituted  into  the  series  expansion  for  gC"^): 

gCo;)  =  -alf  _aTi_  1.3.5"(;2n-f3l  ^n        (5,3) 
6rn2^c2^Po^  2-4-6"(2n+4) 

This  is  the  desired  constitutive  relationship. 

As  a  simple  example  of  this  technique,  consider  again 
the  lTev;tonian  fluid.  Previously  it  was  found  that 

1/t  =  5?6^^3/96pjjaJ^L, 


1P.7 


so  that 

bi   =  55j^V96p,yaiL, 

and 

"bo  =  b?.  =  b3  =  b4  =   ..   =0, 

Hence , 

ai      =     3bi      =      5&ni^/32>ijjaA'L. 

Recalling  that  Po  =  ^n/^^  ^^^  substituting  into 
equation  (5»3)5 

g(T)  =  -z/jx^, 
as,  of  course,  it  should. 

If  the  analytical  form  of  the  experimental  da.ta  can 
be  induced,  then  equation  (5.2)  nay  be  solved  by  quadrature 
without  the  necessity  for  a  poxN^er  law  expansion.  Recall 
that  differentiation  of  the  viscometer  equation,  equation 
(5.1),  yielded 

d(PVt)/dt  =      (P/TTai)£52  g(5P/2L)  de. 
Setting 

-^m  =  P5m/2L' 

and  recalling  that 

T-^,;  =  P6j^(l-cos  e)/4L, 

substitute 

2t,/Tjj.  =  (1-cos  e) 
and 

^::L =  ae 

to  obtain: 
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P       dp  Tiai'TmS  4     "  -Z^^^/ 

The  left  side  of  this  equation  may  also  be  expressed 
in  terms  of  z-^i 

1  d(PVt)      ^     1  d(TmVtni) 

P       dP  \       dTin  ' 

v/here  the  subscript  on  the  symbol  tj^  is  a  reminder 
that  in  place  of  t(P)  one  must  iTrite  t{2'Li:^/h-^) . 

The  result  of  this  manipulation  is  an  equation  of 
the  form 

^    dT^     -  na^  -io  -^w   g^S/^  y^p^^ 

which  is  recognizable  as  Abel's  equation  for  which  the 
solution,  as  given  by  Tricomi  (1957?  pp.  39-40)  among 
others,  is: 


Xr,, 


dT'T 


v%  '^o  V^m  n-; 


where 


«V  =  ^i<^  (5.., 

°m     '-^  ^m 


and 


T'(s.)  =  '^^'^"^^ 


d'^m 


'^-='^.; 


"m  "n'; 
The  corresponding  constitutive  equation  is  sinply 

^^^^       t2   "^  1:3,2  Jo^  ^  ^•■^''  ^T-^;  •  ^^•^'' 

To  illustrate  the  use  of  this  equation,  consider 
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for  the  third  time  a  Newtonian  fluid.  Again, 

1/t     =      5?6^V96jijja^L, 
so 

in  which  case, 

T(0)/t2    =    lim  (Tjj,^0)      5'i^ffiVl6;i^.jT2  (-^  < -^ci^* 

=     0 

T'C-t,,)     =     l5'^,Vl6;aj,, 
and 

as,   again,    it  must. 


*Recall  that  t^  =  P6^/2L  is  the  mai:imiim  value  of  the 

^U  ill 

shear  stress  anyi-ihere   in  the  annulus. 


130 


Experimental  Results 

Besides  the  measurements  made  to  determine  the  rheo- 
logical  constants,  several  other  investigations  v;ere 
carried  out.  These  included  the  change  in  apparent 
viscosity  of  Milling  Yellox/  with  time  (Appendix  H)  and 
with  concentration  and  temperature  (Appendix  I).   It 
was  found  that  the  range  of  the  Hoppler  Rheo-Viscometer 
could  be  extended  to  viscosities  less  than  the  minimum 
of  4  centipoise  recommended  by  the  manufacturer  by 
observing  the  restriction  (Appendix  J): 

P/}ia^  <  15.?  gm/(cm-cP)2. 
Confirmation  of  the  calibration  of  the  instrument  using 
distilled  water  (Appendix  E)  is  an  example  of  this  usage. 

During  the  time  that  the  medium  was  being  used  to 
measure  the  birefringence,  as  described  in  Chapter  Four, 
nine  samples  were  taken.  Each  of  these  v/as  tested  in  the 
Hoppler  Rheo-Viscometer  at  20°  and  2^^   C.  These  data 
are  tabulated  in  Tables  XIX  to  XXVII.  The   fall  times 
listed  are  for  a  displacement  of  3 "00  cm.  The  temper- 
atures are  those  of  the  precision  thermometer  located 
adjacent  to  the  viscometer  cylinder  in  the  water  bath  of 
the  instrument. 
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TA3T,R 

XIX 

VISCOMST 

ER  KEASUREKENTS  -  RUIT 

42* 

P,   gin/cm2 

t,    sec 

T,    °C 

P,    gm/cin2 

t,    sec 

T,    °G 

160 

3.2 

20.00 

80 

4.7 

25.00 

3.4 

20.02 

4.6 

25.00 

3.4 

20.02 

4.7 

25.00 

3.4 

20.02 

3.8 

20.02 

80 

6.3 

20.02 

60 

6.7 

25.00 

6.0 

20.02 

6.8 

25.00 

6.1 

20.00 

6.8 

25.00 

6.0 

20.02 

6.5 

20.02 

40 

16.8 

20.02 

40 

12.1 

25.00 

16.4 

20.02 

12.1 

25.00 

16.6 

20.00 

12.1 

24.98 

16.8 

20.02 

16.5 

20.02 

30 

26.9 

20.00 

30 

18.6 

25.00 

26.^ 

20.00 

18.6 

25.00 

26.3 

20.00 

18.6 

25.00 

26.6 

20.00 

27.0 

20.00 

20 

54.8 

20.00 

20 

35.5 

25.00 

54.9 

20.00 

35.6 

25.00 

55.6 

20.00 

35.7 

25.00 

54.5 

20.00 

54.4 

20.00 

15 

58.6 
53.5 

25.00 
25.00 

15 

99.4 
96.8 

20.00 
20.00 

58.5 

25.00 

99.3 

20.00 

10 

125.3 

25.00 

95.1 

20.00 

128.1 

25.00 

95.6 

20.00 

128.3 

25.00 

10 

264.7 
249.1 

20.00 
20.02 

1.5 

237.2 

25.00 

240.  "^r 

20.00 

5 

603.1 

25.00 

241.2 

20,00 

604.4 

25.00 

^Sample  42  contained  sedinentation.  This  caused  varia- 
bility in  the  run  at  20°  C.  The  rim  at  25^'  C  was  conducted 
with  fluid  decanted  fron  the  toD  of  the  sample  bottle. 
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TABT.F; 

XX 

VISCC 

)L1ETER  MEASUR] 

SI4EII 

p.    gn/cn2 

t,    sec 

T, 

°C 

Li 

160 

3. 

,A 

20. 

.00 

3. 

.5 

20, 

.00 

3. 

.4 

20. 

,00 

80 

6, 

,1 

20. 

.00 

6. 

,1 

20. 

,00 

6. 

.2 

20. 

,00 

40 

17. 

.2 

20. 

,00 

17. 

.2 

20. 

.00 

17. 

,2 

20. 

.00 

30 

27. 

.4 

20, 

.00 

27. 

.4 

20. 

,00 

27. 

.3 

20, 

.00 

20 

51^ 

,0 

20, 

.00 

57. 

.7 

20, 

.02 

58. 

.6 

20, 

.00 

58. 

.0 

20, 

.00 

51. 

.1 

20, 

.00 

15 

103. 

.3* 

20, 

.00 

108. 

,ZL* 

20, 

,00 

100. 

.6 

20, 

.00 

100. 

.2 

20, 

.00 

100. 

.2 

20, 

.00 

gm/cc^  t,  sec  T,  °G 

80  5.1  25. CO 

5.1  25.02 

5.0  25.00 

60         6.9  25.00 

7.1  25.00 
7.0  25.02 

40        12.3  25.00 

12.4  25.00 

12.4  25.02 

12.3  25.02 

30         19.1  25.00 

19.1  25.00 

19o3  25.00 


20 


15 


36.8 

25.00 

36.9 

25.00 

36.9 

25.00 

60.5** 

25.00 

60.8** 

25.00 

61.2** 

25.00 

61.5^* 

25.00 

62.9** 

25.00 

60.5 

25.00 

60.7 

25.00 

60.5 

25.00 

10  254.6        20.04  10  136.6  25.00 

254.1        20.04  136.0  25.00 

254,3        20.00  136.8  25.00 

7.5  255.7  25.0c 

5  1853.9      20.00  5  672.4  25.00 

666.6  25.0c 


*Re.jected  as  spurious 

**Rejected  as  spurious  after  a  pause  to  clean  viscometer 

ball  produced  a  reproducible  result. 
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TABLE  XXI 
VISCOI^ETER  ISASimEIlEriTS  -  RUK  4l8 


P,  gm/cns 

tj^  sec 

hu 

°C 

P,  gra/cm2 

t,  sec 

T, 

°C 

160 

3.6 

20 

.00 

80 

4 

.9 

25. 

.00 

3.6 

20 

.00 

4, 

.9 

25. 

.00 

3.5 

20 

.00 

4, 

.9 

25. 

.00 

80 

5.9 

20 

.00 

60 

6, 

.7 

25. 

.00 

?-9 

20 

.00 

6 

.8 

25. 

.00 

6.0 

20 

.00 

6 

.7 

25. 

.00 

6.0 

20 

.00 

40 

16.6 

20, 

.00 

40 

11, 

.9 

25. 

,00 

16.6 

20, 

.00 

11, 

.9 

25. 

.00 

16.6 

20, 

.00 

11, 

a 
» / 

25. 

.00 

30 

26.4 

20, 

.00 

30 

18. 

.3 

25. 

.00 

26.5 

20. 

.00 

18. 

.2 

25. 

,00 

26.5 

20, 

,00 

18. 

.3 

25. 

,00 

20 

56.3 

20. 

.00 

20 

34. 

.6 

25. 

,00 

54.0 

20. 

.00 

34. 

.8 

25. 

,00 

53.8 

20. 

.00 

34. 

c 

25. 

,00 

54.0 

20. 

.00 

53.8 

20. 

.00 

15 

92.8 

19. 

.96 

15 

56, 

.3 

25. 

,00 

92.7 

20. 

.00 

56. 

.3 

25. 

^QQ 

93.4* 

20, 

.00 

56. 

0 

*  ^ 

25. 

00 

92.8 

20. 

.00 

92.9 

20, 

.00 

10 

220.8 

20. 

,00 

10 

118. 

.3 

25. 

00 

221.0 

20, 

,00 

119. 

,0 

25. 

00 

220.7 

20, 

,00 

119. 

,1 

25. 

00 

7.5 

445.0 

20. 

,02 

7.5 

210. 

,0 

25.00 

5 

1253.1 

20. 

,00 

5 

485. 

,8 

25. 

00 

1240.7 

20. 

,00 

482. 

1 

25. 

00 

*Rej acted  as  spurious 
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TA3I£  XXII 
VISCOMETER  I-IEASIREI-SKTS  -  RUIJ  419 


P,    gn/cin2 

t,    sec 

T,    °C 

P,    gffi/cn2 

t,    sec 

T,    °C 

160 

3.5 

3.4 

3.5 

20.00 
20.00 
20.00 

80 

4.9 
4.7 
4.7 

25.00 
25.00 
25.00 

80 

5.9 
6.0 
6.0 

20.00 
20.02 
20.02 

60 

6.7 
6.9 

6.8 

25.00 
25.00 
25.00 

40 

16.4 
16.3 
16.3 

20.00 
20.00 
20.00 

40 

12.0 
12.0 
12.0 

25.00 
25.00 
25.00 

30 

25.9 
25.8 
26.0 

20.00 
20.00 
20.00 

30 

18.4 

18.3 
18.2 

25.00 
25.02 
25.00 

20 

52.0 
52.1 
52.1 

20.00 
20.00 
20.00 

20 

34.6 
34.6 
34.7 

25-00 

25.00 
25.00 

15 

90.3 
89.6 

89.7 
89.8 
89.8 
89.8 

20.00 
20.00 
19.98 
20.00 
19.98 
20.00 

1? 

55.7 
55.7 

25.00 
25.00 
25.00 

10 

211.4* 
213.3 
214.3 
214.5 

20.04 
20.00 
20.00 
20.00 

10 

116.4 
116.9 
117.1 

25.00 
25.02 
25.00 

1>5 

419.6 

20.04 

7.5 

205.6 

25.02 

5 

1150.9 
1157.1 

20.00 
20.00 

5 

474.9 
477.2 

25.00** 
25.00 

*Rej acted  as  spurious 

**Difficulty  with  circulator.  Temperature  varied  during 

Tvn   from  24.94°  to  25.04°  C. 
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TAET.K 

XXIII 

VISC012]TER  MEASUREI-SITTS  -  RM 

420 

r,    gm/cn^ 

t,    sec 

T,    °C 

P,    gm/cin2 

t,   sec 

T.    °C 

160 

3.7 

3.6 

3.5 

20.00 
20.00 
20.00 

80 

5.1 
5.1 
5.1 

25.02 
25.02 
25.02 

80 

6.0 
6.1 

6.1 

20.00 

20.00 
20.00 

60 

6.8 
6.8 
6.9 

25.00 
25.00 
25.00 

40 

17.0 
16.9 
16.9 

20.00 
20.00 
20.00 

40 

12.1 
12.1 
12.2 

25.00 
25.00 
25.00 

30 

26.9 
27.0 

27.1 

20.00 
20.00 
20.00 

30 

13.6 
18.6 
18.8 

25.00 
25.00 
25.00 

20 

54.8 
55.1 
55.1 
55.0 
55.2 

20.00 
20.00 
20.00 
20.00 
20.00 

20 

35.4 
35.4 
35.4 

25.00 
25.00 

10 

227.1 
227.1 
227.4 

20.00 
20.00 
20.02 

10 

122.8 

123.3 
122.9 

25.00 
25.00 
25.00 

1.5 

440.3 

20.00 

1.5 

214.0 

25.00 

5 

1220.3 
1218.2 

20.00 
20.00 

5 

497.6 
497.4 

25.0c 

25.  CO 

'Meas-urenent  omitted 
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TABT.F, 

XJZIV 

VISCOMETJhlK  MEASUREI-lEIITS  -  RUli 

420B 

P,    gm/cEi2 

t,    sec 

T,    °C 

P,    gm/cm^ 

t,    sec 

T,    °C 

160 

3-7 

3.6 
3.7 

20.00 
20.00 
20.00 

80 

5.0 
5.0 
5.0 

25.00 
25.00 
25.00 

80 

6.0 
6.1 
6.0 
6.0 

20.00 
20.00 
20.00 
20.00 

60 

6.8 
6.8 
6.8 

25.00 
25.QO 
25.00 

40 

16.9 
17.0 
16.9 

20.00 
20.00 
20.00 

40 

12.0 
12.0 
12.0 

25.00 
25.00 
25.02 

30 

27.1 
27.1 
26.9 

20.00 
20.00 
20.00 

30 

18.4 
18.4 
18.4 

25.00 
25.00 
25.00 

20 

55.3 
55.3 
55.2 

20.00 
20.00 
20.00 

20 

34.9 
34-. 8 

34.9 

25.CO 
25.00 
25.00 

15 

96.0 
96.4 
96.2 
96.3 
96.3 

19.98 
20.00 
20.00 
20.00 
20.00 

15 

56.4 
56.3 
56.5 

25.00 
25.00 
25.00 

10 

230.0 
230.7 
231.5 
230.3 
230.7 

20.00 
20.00 
20.00 
20.00 
20.00 

10 

119.6 
120.2 
120.2 

25.02 
25.00 
25.02 

1.5 

448.7 

20.00 

7.5 

214.1 

25.00 

5 

1237.5 
1236. A 

20.00 
20.00 

5 

508. 3 
500.0 

25.00 
25.00* 

^Difficulty  with  circulator.  Temperature  varied  during 
run  from  24.98° to  25.02°  C. 
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TABT,^ 

XXV 

VISCOI-STER  ^T;/lSURElEEKTS  -  RUII 

424 

P,    gm/cm^ 

t,    sec 

T,    °C 

P,    gm/cn^ 

t,    sec 

T,    OC 

160 

3.6 

3.5 
3.4 

20.00 
20.00 
20.00 

80 

4.9 
4.8 

4.9 

25.00 
25.00 
25.00 

80 

6.3 
6.3 

6.4 

20.00 
20.00 
20.00 

60 

7.1 
7.1 
7.1 

25.00 

25. CO 

25.00 

40 

17.8 
17.8 
17.8 

20.00 

20.00 
20.00 

40 

12.6 
12.7 
12.6 

25.00 
25.02 
25.00 

30 

28.5 
28.6 
28.4- 

20.00 
20.00 
20.00 

30 

19.6 
19.6 
19.5 

25.00 
25.00 
25.00 

20 

58.7 
58.8 
58.8 

20.00 
20.00 
20.00 

20 

37.5 
37.7 
37.6 

25.00 
25.00 
25.00 

15 

102.0* 
102.3* 
102.5^- 
102.8* 
103.1* 

20.00 
20.02 
20.00 
20.00 
20.00 

15 

61.8 
61.7 
61.6 

25.00 
25.00 
25.00 

10 

254.2 

253.9 
253.4 

19.98 
19.98 
20.00 

10 

134.6 
134.8 
134.9 

25.00 
25.00 
25.00 

7.5 

51c.  8 

20.00 

7.5 

243.1 

25.00 

5 

1480.8 

20.00 

5 

591.3 
591.4 

25.00 
25.01* 

*xhe  systematic  increase  in  these  valu.es  is  unexplained. 
*=^Difficulty  v/ith  circulator.  Temperatures  varied  during 
ruji  from  24.98°  to  25.04°  C. 
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TABT.E 

JJNl 

VISCOMETER  MEASUF.EMEITTS  -  RUIJ 

424B 

P,    gn/cm2 

t,    sec 

T,    °C 

P,    gm/cin2 

t,    sec 

T,    °C 

160 

3.3 
3.4 
3.4 

20.00 
19.98 
19.98 

80 

5.0 
4.9 
5.0 

25.00 
25.00 
25.00 

80 

6.4 
6.4 
6.4 

20.00 
20.00 
20,00 

60 

7.3 

7.3 
7.3 

25.00 
25.00 
25.00 

40 

18.2 
18.2 
18.2 

20.00 
20.00 
20.00 

40 

12.9 
13.0 
13.0 

25.00 
25.00 
25.00 

30 

29.3 

29.3 
29.3 

20.00 
20.00 
20.00 

30 

20.1 
20.1 

20.1 

25.00 
25.02 
25. CO 

20 

60.9 
60.9 
61.1 

20.00 
20.00 
20.00 

20 

38.9 
39.0 
39.0 

24.94 
25.00 
25.00 

15 

107.3 
107.1 
106.9 

20.00 
20.00 
20.00 

15 

64.3 
64.4 

64.3 

25.00 
25.00 
25.00 

10 

267.4 

271.3 
268.0 

20.00 
20.00 
20.00 

10 

141.1 
139.6 
140.7 

25.02 
25.02 
25.00 

1.5 

549.4 

20.00 

7.5 

257.6 

25.00 

5 

1632.6 

20.00 

5 

634.8 
640.6 

25.02 
25.00 
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TABLE  XXVII 
VISCOL^TER   LiEASroSI-rSIITS  -  RUN   425 


P,   gm/c!n2 

t,    soc 

T,    °C 

P,    gm/cm^ 

t,   sec 

T,    oQ 

160 

3.5 
3.5 
^.5 

20.00 
20.00 
20.00 

80 

5.1 
4.9 
5.0 

25.00 
24.98 
24.96 

80 

6.4 
6.4 

6.4 

20.00 
20.00 
20.00 

60 

7.3 
7.4 

7.3 

25.00 
25.00 
24.98 

40 

18.3 
18.2 
18.3 

20.00 
20.00 
20.00 

40 

13.1 
13.2 
13.2 
13.1 

25.00 
25.00 
25.00 
25.00 

30 

29.5 
29.6 
29.6 

20.00 
20.02 
20.00 

30 

20,3 
20,3 

20.4 

25.00 
25.00 
25.00 

20 

61.5 
61.6 
61.8 

20.00 
20.02 
20.00 

20 

39.6 

25.00 
25.00 
25.CO 

15 

109.7 
109.6 
109.6 

20.00 
20.00 
20.00 

15 

65.4 

65.3 
65.4 

25.00 
25.00 
25.CO 

10 

273.8 
274.6 
274.6 

20.00 
20.00 
20.00 

10 

144.2 
144.7 
145.0 

25. CO 

25.00 
25.00 

1*5 

562.2 

20.01* 

7.5 

264.9 

25.00 

5 

1661.7 

20.00 

5 

656.2 
660.8 

25.00 

'Difficulty  with  circulator.     Temperature  varied  durir.^ 
'un  from  20.00°    to   20.02°    C.  ^^--^u  uux  ^„^ 
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Discussion  of  Results 

This  discussion  is  in  two  sections.   The  first 
describes  the  reduction  of  the  data  to  an  analytic  form 
by  the  method  of  least  squares.  The  second  section  is 
concerned  with  the  calculation  of  the  rheological  constants 
for  Milling  Yellov;  from  this  analytic  form. 
Curve  Fitting 

Three  of  the  experimental  runs  are  plotted  in 
Figure  23.   Inspection  of  the  preceding  tables  indicates 
that  the  plotted  data  are  typical. 

In  all  of  v;hat  follows  the  values  of  P  which  appear 
in  Tables  XIX  to  XXVII  have  been  reduced  by  0.28  gm/cm^ 
to  correct  for  instrujnental  friction.  The  need  for  this 
correction  is  provided  in  Appendix:  E, 

The  data  in  Figure  23  appear  to  have  a  slope 
d(ln  P)/d(ln  t)  v/hich  has  a  value  of  -2  at  lox-;  values  of 
P  and  i/hich  approaches  -1  as  P  becomes  large.  A  straight 
line  of  slope  -2  has  been  passed  trirough  txhe  data  for  the 
lower  values  of  P  (run  4-25)  •  The  logarithm  of  the 
difference  between  the  experimental  value  of  t  and  the 
value  obtained  from,  this  straight  line  v;as  plotted  versus 
In  P  and  a  nearly  straight  line  was  obtained.  This 
suggested  a  test  of  the  furLctional  form 

1/t  =  P2/(ki+k2?) 
or 


p 


Pt  =  -~  +  k2  . 


1  000  _ 


Fall 

time,    100 

sec 


~ 


1  0 


'9 


\S 


to 


1  00 


P,  gram/era^ 
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■^■^'^^'2^23.  Response  of  viscometer,  runs:  ©  42;  CD  420- 


O  425. 
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The  scale  of  the  log-log  plot  (Figure  23)  is  too 
small  to  provide  a  serious  test  of  the  hypothetical 
functional  form  just  given,  so  the  data  xrere  replotted 
on  a  cartesian  scale  as  Pt  versus  1/P.  The  result  sho-^'m 
in  Figure  24-  is  for  the  same  experimental  runs  plotted 
previously.  The  predicted  linearity  is  reasonable  only 
at  higher  values  of  P, 

Again  the  departure  from  the  predicted  functional 
form  served  as  the  new  variable  for  improved  curve- 
fitting.  Following  the  example  of  Lindgren  (1957)?  who 
vras  interested  in  the  viscosity  at  infinite  shear  rates, 
the  value  of  the  product  Pt  as  1/P  approaches  zero  was 
taken  from  Figure  24,  and  a  new  correlation  was  attempted 
betvreen  [pt  -  (Pt)oql  and  P,  where  (Pt)co  is  the  value 
for  Pt  as  P  increases  without  limit.  As  sho\.;-n  in  Figure 
25,  a  linear  result  is  obtained  vrhen  P[Pt  -  (Pt)oo]  is 
plotted  versus  1/P,  The  corresponding  functional  form  is 

Pt  =  (Pt)oo+  ki/P  +  k2/P2 
or 

1/t  =  p3  [(P+Ci)"^  -  (P+C2)~VCo. 
Values  of  the  constants  Cq,  Ci  ,  an.d  Co   are  given  in 
Table  XICVIII  for  each  of  the  experimental  runs.  Also 
ta.bulated  are  the  viscosity  values  corresponding  to  (Pt)oo  , 

>.ico  =   0.00534(Pt)„=   0.0053^  (C2-Ci)/Go, 

and  an  estimate  of  the  goodness  of  fit  of  the  functional 
form. 
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FIGURE   24-.      Test  of  functional  for:n:    Pi 
0    ^2;    0     420;    O     425. 
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TABLE  xr/III 
RHEOLOGICAL  CO?ISTAIITS 


Run 

Temp, 

Co 
7733 

Ci 

1.91 

C2 
P2.96 

>l005 

cP 
1.96 

Percent 

error 
in  1/t* 

42 

20 

0.95 

25 

4013 

4.68 

18.45 

1.56 

1.75 

AA 

20 

18892 

-16.52 

20.16 

2.75 

13.7 

25 

2173 

7.36 

14.54 

1.62 

1.19 

4-18 

20 

3301 

8.38 

19.06 

1.90 

3.31 

25 

7262 

1.06 

28.14 

1.43 

1.00 

419 

20 

3669 

9.49 

20.18 

1.83 

1.27 

25 

6660 

1.29 

24.86 

1.51 

0.75 

420 

20 

5072 

8.54 

23.18 

1.85 

2.27 

25 

7479 

1.06 

28.65 

1.45 

1.13 

420B 

20 

6401 

7.66 

26.95 

1.77 

1.39 

25 

6555 

1.65 

25.19 

1.49 

1.13 

424 

20 

7Q^5 

2.51 

22.77 

2.07 

4.04 

25 

6931 

2.27 

26.86 

1.51 

1.17 

424B 

20 

10052 

-0.79 

24.26 

4.28 

5.17 

25 

6567 

2.91 

25.39 

1.56 

0.66 

425 

20 

9042 

1.35 

24.68 

2.07 

4.25 

25 

6846 

2.93 

26.63 

1.54 

1.12 

*I-iaxiEiura  percentage  deviation  \rhen   measured  values  of  1/t 
are  compared  with  values  calculated  using  the  estimate 

1/t  =  P^[(P+Ci)""'-(P+C2)-']/Co. 
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Determination  of  gCx)   from  Experimental  l-leasiirementi 
Previous  analysis  yielded  equation  (^.5) 


-Z^  TV2    Jo  '        */'^-'^V7 


where 

a^Tm  dCTa^/tm) 


T(T^)        = 


&ra^  d'^m 


Using  the  functional  form  for  1/t  obtained  in  the  previous 
section  and  recalling  that 

t(P)   =  t,3(2LTj^/6j3), 
it  follovrs  that 


vrfiere 

do  =  Co6y2L, 

cii  =  Ci6y2L, 

^2  =  C26jjj/2L. 


iience , 


m  6jn3  \   do     LT,+di        T^^-dsJ 

'^m^  r      1  1       1 


do   L(T^+di)~         (^n^cl2)^J 


Y(0)      =      0 


J 


and  g(T)    is   given  in  integral  form  by: 


1^7 


gi^)      = 


2aJfL 


doS^^-r^/z  - 


-  1  \  25t. 


w 


^n 


'^w'^^1   'H;"^d2 


1 


1 


2t, 


w 


J l\    dT,.; 


(V^dTr  '  ('t^^+ci.2)^ji  ^A^ 


The  integration,  though  tedious,  is  not  difficult  and 
yields  the  desired  relationship  betv;een  shear  stress  and 
shear  rate  for  Milling  Yellow: 

3.JIL     r   3di5       3d2^ 


g(T)   = 


do6^3  1  T;(T:+(i^)2     T(T+d2)' 


(5.6) 


■'m 


6628di3       6628d2  3 

105  T         105  -c 


14368di2(T+di)       14368d2^(^+d2) 


10  5t 


10  5t 


374-4di(T+di)2    .    3744d2(T--^d2)^     ^ 
35^  "  35^ 

1024(T+di)^        1024(T;-d2)^      , 

■  — — —  ■     ■     ■     ■     >  ^      — — --    -■-■  -  -■  - ■  T 

35-^ 


371^ 


3d/ 


UA3(T-^di')?     y-rBCT+dO^ 


'  tanh"^  ^/(-u+di  )   - 


[  3d; 


2d: 


1  . •^ 

!  tanh"^  -A/(T+d2)  > 


[yT3(T+d2y?       yT3(T+d2)3j 

At  very  high  values  of  t.  this  becomes  simply 
96a^L(d2-di)'^ 


gCi^oo)  = 


5do6^3 


14- 


p 


At  very  low  values  of  t, 


5 


g(,,)   .  _aiL  2560  (d,-d,)  ^3       (^^0) 
or  in  terms  of  the  constarxts  tabulated  in  Table  XX7III, 


^^  ^     Co&m?   63     G1C2 

In  the  first  part  of  this  chapter  it  was  found  that 
the  data  of  Peebles,  Prados,  and  Honeycutts  (1965)  could 
be  represented  by  the  two-constant  empirical  relationship 

f    >     JtiQoT^  +  (31v^2) 
g(T)  = —  t3  {5»7) 

\-7hich  had  the  limiting  value  as  t  becomes  small: 

g(T)   =   3T3/2k.  (t-0) 

Comparison  of  this  equation  \-jlth   the  one  just  above 
provides  the  constant  k  to  be  used  in  the  tvro- constant 
relationship  based  upon  the  present  data: 

l^     =       21  Co5k?  ^1^2  (^  8) 

5120  aa3  (C2-C1) 

A  similar  comparison  for  very  high  rates  of  shear  yields 
the  value  of  ;ico  previously  cited  and  tabulated: 

>ico  =  ^-^  (Pt)a,  =     0.00534  %:§^  .     (5.9) 
9d  a/(L  C2O1 

Several  typical  plots  of  the  Theological  response 
of  Milling  Yellow  are  provided  in  Figure  26,  based  upon 
equation  (5.6).   Figure  27  provides  a  comparison  between 
equations  (5.6)  and  {5»7)»      Clearly  k  and  jicd  are  poorly 
predicted  by  equations  (5.8)  and  (5.9). 
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It  should  be  noted  that  whatever  form  for  gC-r)  is 
chosen  for  a  given  application,  there  is  a  requirement 
that  do,  di ,  and  da  (Co,  Ci ,  and  C2)  be  positive.  This 
requirement  fails  in  runs  44  and  424B  at  20°  C.  There  is 
no  explanation  for  the  anomalous  measurements  of  run  44 
which,  due  to  the  poor  fit  indicated  by  the  last  coliomn 
of  Table  XXVIII,  was  not  analyzed  further.   In  run  424B 
the  fit  of  the  data  is  only  slightly  poorer  than  that  for 
runs  424  and  42^  suggesting  that  these  data  are  valid. 
Since  Ci  in  this  case  is  small,  it  seems  appropriate  to 
replace  C-,  by  zero  in  the  fit  of  the  data  to  obtain: 

1/t  =  p2(^:~^)C2/Co. 

This  relationship  may  now  be  used  in  equation  (5.5) 
to  obtain  the  corresponding  relationship  between  shear 
stress  and  shear  rate.  This  has  not  been  done  in  the 
present  case. 

Range  of  Application 
The  data  of  previous  investigators,  beginning  with 
Prados  (1957),  have  indicated  that  the  viscosity  remains 
constant  when  the  shearing  stress  is  less  than  about  5  to 
10  dyne/cm^  for  1.3  to  1.7  percent  solutions  at  temper- 
atures around  20°  to  25°  G.  This  point  was  discussed  at 
some  length  in  Chapter  Two.  It  would  be  helpful  if  the 
experimental  measurements  reported  in  this  chapter  were 
able  to  cast  some  light  into  this  area  of  concern.  Unfor- 
tunately, the  Hoppler  Rheo-Viscometer  is  an  inappropriate 
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instrument  for  the  measurement  of  viscosities  at  lower 
shear  stresses.  Even  when  the  force  applied  to  the  ball 
is  reduced  to  its  minimum  value,  so  that  the  ball  falls 
as  slowly  as  0,025  mm/sec,  the  shearing  stress  at  the 
wall  of  the  annulus  will  reach 

T,   =  £&  _  (5x981  dvne/cmM (0.007  cm) 
V     2L  (2) (0.1365  cm) 

=  125  dyne/cm^ 
according  to  the  previously  derived  relationship.  Although 
the  stress  will  become  smaller  at  other  points  on  the  wall 
and  drops  linearly  to  zero  as  one  passes  from  the  v/all  to 
the  center  of  the  annulus,  the  shearing  stress  is  much 
higher  than  the  5  to  10  dyne/cm^  value  at  v/hich  a  constant 
viscosity  might  be  observed. 

Evidence  for  a  constant  viscosity  at  low  shear 
stress  can  be  inferred  from  the  success  of  the  technique 
in  Chapter  Four  which  was  used  to  calculate  the  bire- 
fringence based  upon  the  assumption  that  at  low  flow  rates 
the  viscosity  would  remain  constant.  More  direct  data 
are  available  in  the  form  of  manometer  readings  made  during 
the  last  two  of  the  birefringence  measurements. 

According  to  Cornish  (1928)  the  pressure  gradient  in 
a  rectangular  conduit  flowing  a  Newtonian  fluid  will  be 

dp/dz  =  -  26iU}ijf/62So 
where  6i ,  52 »  and  So  are  channel  parameters  defined  in 
Chapter  Two  with  numerical  values  tabulated  in  Appendix  G, 
TT  is  the  average  flow  rate,  and  )i^   is  the  viscosity.  The 
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manometer  estimate  of  the  pressure  gradient  is  well  known: 

dP/dz  =  (pf  -  p)g(,h/AZ 
where  p^  is  the  density  of  the  manometer  fluid,  p  is  the 
density  of  the  medium,  Zq   is  the  gravitational  constant, 
h  is  the  difference  in  fluid  levels  in  the  manometer,  and 
t.z   is  the  distance  between  the  manometer  taps. 

Values  for  the  manometer  readings,  flow  rates,  and 
corresponding  viscosities  are  tabulated  in  Table  XXIX. 
If  the  first  two  readings  of  run  4-24  are  discarded  due 
to  their  anomalous  temperatures,  the  remaining  data  for 
both  trials  indicate  a  reasonable  likelihood  that  constant 
viscosities  in  the  neighborhood  of  100  and  125  cP  were 
observed  at  the  beginning  of  the  two  runs.  It  must  be 
noted  that  this  conclusion  could  be  invalidated  if  a 
previous  ass-umption  —  that  normal  stress  functions  are 
negligible  for  Milling  Yellow  solutions  —  fails. 

When  the  response  of  a  Newtonian  fluid  having  a 
viscosity  of  100  cP  is  superposed  on  the  previous  plot 
of  Milling  Yellow's  rheological  response  given  in  Figure 
27,  it  will  be  noted  that  the  intersection  occurs  at 
T  =  20  dyne/cm^,  about  twice  the  upper  end  of  the  linear 
range  observed  by  Prados  (1957)  and  others.  In  the  absence 
of  more  convincing  evidence,  it  will  be  assumed  that  this 
intersection  is  valid,  and  the  following  relationship 
for  g(T:)  is  recommended  for  Milling  Yellow  over  the  full 
range  of  shear  stresses : 


TABLE  XXIX 
ESTBIATION  OF  APPARENT  VISCOSITY 
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Run  424: 

Flow  rate, 
gm/sec 

Temp, 
24.9 

Manometer 
head ,  cm* 

8.13 

Apparent 

viscosity, 

cP 

2.50 

157.5 

5.53 

25.3 

8.26 

72.0 

6.34 

24.4 

13.21 

100.5 

9.09 

24.3 

18.80 

99.8 

12.82 

24.3 

23.62 

88.8 

20.78 

24.5 

28.96 

69.8 

25.57 

24.2 

32.77 

60.8 

Run  425: 

Flow  rate, 
gm/sec 

Temp, 
24.2 

Manometer 
head ,  cm* 

5.33 

Apparent 

viscosity, 

cP 

1.99 

129.0 

3.50 

24.2 

9.14 

125.8 

3.60 

23.8 

9.91 

132.5 

5.71 

23.8 

13.72 

115.8 

7.78 

23.8 

17.53 

109.5 

12.64 

23.8 

21.59 

82.2 

16.00 

23.8 

29.46 

88.8 

♦Manometer  fluid  was  reagent  grade  CCl^,  S.G.  =  1,584; 
manometer  taps  were  separated  by  200.0  cm. 
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gCt)  from  equation  (5.d)   '^  ^ '^c 


where 
and 


T   =  20  dyne/cm^. 


When  convenience  of  form  is  an  important  consid- 
eration, gC^)  may  be  obtained  from  equation  (5«7)  over 
the  upper  range  of  shear  stresses;  however,  the  constants 
in  the  equation  must  be  obtained  by  experiment  since 
comparisons  of  equations  (5«6)  and  {5*7)   at  very  low  and 
very  high  shear  rates  do  not  provide  adequate  estimates 
for  the  interrelationship  between  the  constants  in  the  two 
expressions. 


CHAPTER  SIX 

DISTRIBUTION  OF  SHEAR  RATES 
IN  RECTANGULAR  CONDUITS 


Even  if  the  fimctional  dependence  of  the  Theological 
and  optical  properties  upon  shear  rates  is  well  defined 
for  a  birefringent  medium,  the  prediction  of  fringe 
patterns  for  steady  axial  flow  through  a  given  cross- 
section  requires  a  hypothetical  velocity  distribution, 
within  the  channel.  Once  this  is  obtained,  the  distri- 
bution of  shear  rates  follows  by  differentiation. 

The  velocity  distribution  for  Milling  Yellow  flowing 
in  rectangular  conduits  has  not  been  calculated  and  is 
not  determined  in  this  dissertation.  The  distribution  of 
other  fluids  in  rectangular  sections  has  been  discussed 
at  length  in  Chapter  Two.  The  present  chapter  evaluates 
these  distributions  in  terms  of  their  applicability  to 
Milling  Yellow. 

Power- Lav  Fluids 

Of  the  non- Newtonian  materials  for  which  flows  in 
rectangular  pipes  have  been  studied,  only  power-law  fluids 
have  received  significant  attention.  Since  there  is  some 
evidence  that  Milling  Yellow  has  a  behavior  similar  to  a 
power- law  fluid  over  a  certain  range  of  lower  shear  rates. 
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it  appears  appropriate  to  compare  predicteO  pressure  drops 
for  power- law  substances  with  the  pressure  drop  measured 
in  the  rectangular  conduit  with  Milling  Yellow. 

The  relationship  g  =  g{i)   has  been  plotted  in 
Figures  26  and  27  for  run  424B,  At  shear  rates  less  tten 
10  sec"^  these  data  are  well  represented  by  the  power- law 
expression: 

T  =  7.35  g^^. 
Substituting  K  =  7.35  and  m  =  1/3  into  equation  (2,5)? 
obtained  by  Wheeler  and  Wissler  (1965)  for  square  pipes, 
there  results,  for  mass  flow  rate  G, 

-dp/dz  =  38.1  G^/^ 
where  the  sides  of  the  square  in  equation  (2.5)  have  been 
replaced  by  the  geometric  mean  of  the  average  length  of 
the  walls  in  the  nearly  square  conduit.  For  the  manometer 
described  in  the  previous  chapter, 

-dp/dz  =  2.83  h. 
Hence,  if  Milling  Yellow  behaves  as  a  power- law  fluid, 
the  manometer  readings  should  vary  with  the  flow  rate 
according  to  the  simple  relationship: 

h  =  13.5  gV3. 
The  manometer  readings  predicted  by  this  equation 
are  plotted  in  Figure  28  together  with  the  manometer 
readings  given  in  Table  XXIX  for  runs  424  and  425.  Mano- 
meter values  predicted  for  Newtonian  fluids  having 
viscosities  of  100  and  125  cP  are  represented  by  dashed 
lines  in  the  same  figure.  The  g(T)  relationship  for  run 
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FIGURE  28.  Manometer  readings  compared  with  those 
expected  for  Newtonian  fluid  (dashed  lines  having  the 
indicated  viscosities)  or  a  power-law  fluid.  Runs: 
9  4-24;  (D  425. 
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425  (t  =  1.^  g^/3)  differs  so  little  from  that  for  rim 
424  that  the  separate  manometer  prediction  for  run  425 
(h  =  13.7  GV3)  is  not  plotted. 

It  will  be  noted  that  the  gC-r)  relationships  are 
based  upon  viscometric  measurements  made  at  25°  C,  while 
the  average  temperature  for  runs  424B  and  425  in  the 
conduit  were  24.6''  and  23.9"  C  respectively.  Variations 
in  the  apparent  viscosity  of  Milling  Yellow  with  temper- 
ature (Appendix  I)  average  10  to  15  percent  per  degree 
Celsius.  A  correction  of  this  order  would  yield  a  good 
correlation  between  runs  424B  and  425,  but  the  disparity 
between  the  data  and  the  power- law  model  would  persist. 

Since  the  power-law  model  predicts  pressure  gradients 
inaccurately,  it  is  unlikely  that  an  accurate  velocity 
distribution  can  be  derived  from  the  model.  Yet  even  an 
insufficient  solution  for  the  distribution  may  represent 
an  improvement  over  the  Newtonian  assumption  that  the 
viscosity  is  constant. 

Although  Wheeler  and  Wissler  (1965)  have  criticized 
Schechter's  (1961)  variational  coefficients  as  imprecise, 
Schechter's  solution  for  power- law  flow  in  square  pipes 
has  two  advantages:  the  coefficients  are  tabulated  and  the 
velocity  distribution  is  scaled  to  the  measured  flow  rate 
in  the  channel.   In  a  square  pipe  sjamnetry  reduces  the 
number  of  coefficients  in  Schechter's  solution  (equation 
2.4)  from  six  to  four.  Three  of  these  are  plotted  as  a 
function  of  the  power- law  exponent  in  Figure  29.  The 
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foiirth,  which  cannot  be  plotted  on  a  logarithmic  scale, 
has  values  of  0,0289  at  m  =  1,  0.000?  at  m  =  0.75,  and 
-0.0285  at  m  =  0.5.  Extrapolation  gives  A4  =  -0.0472 
at  m  =  1/3 •  Using  this  estimate  and  obtaining  the 
remaining  coefficients  by  extrapolating  in  Figin-e  29  as 
shovm,  equation  (2.4)  becomes: 

u  =  (4G/6162)  {  2.16  sin(TTx/26i)  sin(Try/262)  + 

0.349  [sin(3TTx/26i)  sin(TTy/262)  + 

sin(Trx/2&i)  sin(3ny/262)  ]  - 

0.0472  sin(3'Tx/26i)  sin(3TTy/262)  + 

0.0683  [sin(5"x/2&i)  sin(TTy/262)  + 

sin(nx/2&i)  sin(5TTy/262)]]  • 

The  shear  rate  is  obtained  directly  by  differentiation. 

Some  Alternatives 
Neither  Cornish's  (1928)  relationship  for  a  Newtonian 
fluid,  equation  (2.3)  ,  nor  the  power- law  approximation 
just  cited  may  be  used  with  confidence  to  predict  the 
velocity  distribution  for  Milling  Yellow  in  a  rectangular 
pipe.  Alternatives  require  the  adoption  of  more  realistic 
rheological  relationships  such  as  those  derived  in  Chapter 
Five  and  the  choice  of  a  suitable  analytical  technique 
for  its  utilization.  These  choices  are  not  independent 
since  the  analytical  technique  will  be  influenced  by  the 
complexity  of  the  differential  equation  which  results 
from  substitution  of  the  rheological  relationship  into 
the  equilibrium  equation.  Of  the  analytical  methods 
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siirveyed  in  Chapter  Two,  those  of  Litvinov  (I968)  and 
Rothemeyer  (1970)  seem  to  possess  the  necessary  generality. 
It  also  seems  necessary  to  resolve  the  possibility  raised 
by  Wheeler  and  Wissler  (1966)  that  significant  secondary 
flows  are  present  in  the  channel. 


CHAPTER  SEVEN 
SUMMARY  AND  CONCLUSIONS 

This  dissertation  has  consisted  of  three  nearly- 
separate  studies  which  have  been  discussed  individually 
in  Chapters  Three,  Four,  and  Five.  The  major  findings 
of  these  chapters  may  be  summarized  as  follov;s. 

Optical  Analysis  of  Birefringent  Flov;  Fields 

An  aqueous  solution  of  Milling  Yellow  provides  an 
excellent  nediiun  for  the  visualization  of  flows  in  a 
variety  of  processes.  At  very  low  flow  rates  the  fringe 
patterns  which  result  from  its  birefringence  may  be  used 
to  calculate  directly  the  velocity  and  pressure  distri- 
butions of  two-dimensional  flows, 

A  like  determination  of  three-dimensional  flow  fields 
directly  from  birefringent  patterns  is  not  practical,  but 
the  inverse  process  —  hypothesizing  a  flow  field  and 
checking  it  from  the  fringe  pattern  —  is  indeed  practical, 
even  at  flow  rates  for  which  the  optical  and  rheological 
properties  of  the  medium  do  not  vary  linearly  with  the 
shear  rateo  In  such  cases  it  is  shown  that  the  fringe 
pattern  will  depend  upon  the  effective  optical  parameters 
•^  and  AN  which  are  in  turn  related  explicitly  to  the 
optical  properties  of  the  medium:  An,  \Ja,  and  0. 
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It  is  ass-umed  that  An  and  -^   are  uniquely  determined 
by  the  magnitude  of  the  shear  rate: 
An  =  An(g), 

and  that 

®  =  tan-t 


5g/ay 
dg/5x 


Any  hypothetical  velocity  which  can  be  differentiated 
will  yield  a  distribution  of  shear  rates  from  which  the 
optical  coefficients  nay  be  calculated  point  by  point. 
The  fringe  pattern  may  then  be  computed  by  an  iterative 
process  which,  in  effect,  replaces  the  medium  by  a 
succession  of  linear  wave  plates. 

Optical  Properties  of  the  Medium 

Standard  solutions  of  Milling  Yellow  cannot  be 
prepared  with  reproducible  properties  from  the  commercial 
dyestuff.  For  this  reason  the  birefringence  An  has  been 
measured  in  the  past  either  just  prior  to  a  given  test 
or  just  after, using  a  two-dimensional  flow  regime.  The 
present  study  shows  that  at  low  flow  rates,  in  the  vicin- 
ity of  the  wall,  the  square  of  the  birefringence  varies 
linearly  with  the  mass  flow  rate,  providing  a  simple  means 
of  determining  either  once  the  other  is  known.  Plots  of 
the  fringe  order  at  the  wall,  squared,  versus  the  flow 
rate  confirm  this  finding,  and  birefringences  for  the 
various  test  samples  are  calculated  from  the  data. 

The  orientation  angle  \{f  is  not  measured  in  the 
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present  study,  but  an  empirical  form  is  provided  for 
previously  published  data. 

Rheological  Properties  of  the  Medium 

Previous  authors  have  demonstrated  that  Milling 
Yellow  is  well  represented,  though  not  uniquely,  by  a 
tJriree- constant  empirical  relationship,  the  Powell-Eyring 
equation.  The  same  data  are  shown  in  the  present  study 
to  be  satisfied  by  a  two-constant  relationship  in  which 
the  general  form  g  =  g(T)  is  used  rather  than  the  previous 
form  T  =  T(g),  There  is  evidence  that  the  two-constant 
form  is  invalid  at  shear  stresses  above  about  500  dyne/cm^. 

Data  are  presented,  both  in  the  text  and  in  the 
appendices,  which  show  the  variation  of  Milling  Yellow's 
rheological  properties  with  time,  with  temperature,  and 
with  concentration.  Based  on  measurements  made  with  the 
HSppler  Rheo-Viscometer ,  a  sliding  ball  viscometer,  rheo- 
logical constants  are  provided  from  which  the  relationship 
between  shear  stress  and  shear  rate  for  Milling  Yellow  may 
be  obtained.  The  resultant  expression,  though  complicated, 
requires  only  three  coefficients  which  are  tabulated  for 
the  samples  in  the  present  investigation. 

The  use  of  the  Hoppler  Rheo-Viscometer  to  obtain 
basic  rheological  information  is  supported  by  an  analysis 
of  the  flow  in  the  instrument  and  a  full  discussion  of  the 
assumptions  upon  which  the  analysis  is  based.  A  principal 
advantage  of  the  method  given  is  that  the  constitutive 
form  for  the  material's  rheological  response  to  stresses 
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need  not  be  known  to  employ  the  instrtunent,  and  indeed, 
this  form  may  be  constructed  by  quadrature  from  the  mathe- 
matical form  of  the  experimental  results.  For  the  Hoppler 
Rheo-Viscometer  it  is  necessary  to  assume  that  flow  past 
the  sliding  ball  may  be  modeled  on  flow  through  an 
eccentric  annulus  having  a  fixed  "effective  length"  which 
is  obtained  by  calibration.   In  a  falling  cylinder  visco- 
meter it  is  hypothesized  that  the  effective  length  will  be 
approximated  by  the  actual  length  of  the  cylinder.  The 
effective  length  of  the  Hoppler  instrument  is  shown  to 
vary  less  than  1  percent  for  a  wide  range  of  shear  stresses 
in  the  annulus  and  for  a  change  in  the  viscosity  of  the 
Newtonian  calibration  fluid  of  more  than  an  order  of 
magnitude.  Solutions  are  provided  to  the  viscometer 
equation  for  Newtonian  fluids,  power- law  substances,  bi- 
linear materials,  and,  as  stated  earlier.  Milling  Yellow. 

Other  Findings 
The  prediction  of  fringe  patterns  requires  a  knowledge 
of  the  distribution  of  shear  rates  in  the  conduit  through 
which  the  fluid  flows.  For  the  special  case  of  a  rectang- 
ular conduit  (in  which  the  birefringence  was  measured) 
methods  for  obtaining  such  relationships  were  reviewed, 
but  a  solution  of  the  problem  was  not  attempted. 


APPENDIX  A 

THE  EFFECTIVE  BIREFRINGENCE  AND  ORIENTATION  AI^GLE 
OF  THE  OPTICAL  ELLIPSE 
FORMED  BY  THE  INTERSECTION  OF  THE  OPTICAL  ELLIPSOID 

WITH  THE  PLAI^IE 
ORTHOGONAL  TO  THE  PATH  OF  LIGHT 

Consider  an  ellipsoid  of  circular  cross-section 
(prolate  spheroid)  with  the  surface 

(x'/nz)^  +  (y'/na)^  +  (z'/ni)^  =  1, 
where  ni  and  n2  are  the  major  and  minor  a^ces  respec- 
tively. In  the  initial  position  the  x'y' z' -coordinate 
system  is  coincident  with  the  xyz- coordinate  system. 
If  the  x'y'z*  system  is  rotated  first  through  an  angle 
©about  the  z'-axis  and  then  through  an  angle  -if   about 
the  y'-axis,  then  the  relation  between  the  coordinate 
systems  will  be: 


X' 

y' 


cos  -^     0  sin  -^ 
0    10 
-sin  -^     0     cos  -^ 


cos  ®  sin  ®  0 

-sin  ®  cos  ®  0 

0      0    1 


X 

y 

z 


If  the  center  of  the  coordinate  system  (x',y'»z')  is 
translated  to  the  point  (xo,yo,zo),  then  the  new  relation 
can  be  obtained  by  substituting  (x-xo)  for  x  and  (y-yo) 
for  y  in  the  column  matrix  to  the  right.  Expanding  the 
matrices  and  substituting  into  the  equation  for  the  surface 
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of  the  ellipsoid  yields: 


1 


[(x-xo)cos  f   COS  9  +   (y-yo)cos  ^lr   sin  9  +  z   sin  a]/-]' 


(n-^)^ 


[(x-xo)sin  Q  +   (y-yo)cos  e] 


2 


Qn  -  -^; 

[-(x-xo)sin  >]/•  cos  ®  -    (y-yo)sin  \}f  sin  ®  +  z  cos  AJr]' 

(n  +  -y) 

where  the  major  and  minor  axes  have  the  dimensions  n  i  -j"* 
The  ellipsoid  so  defined  conforms  to  the  description 
of  the  optical  ellipsoid  shown  in  Figure  3.  The  corres- 
ponding optical  ellipse  is  formed  by  the  intersection  of 
this  surface,  F(x,y,z)  =  0,  with  the  plane  x  =  xq: 

[(y-yo)cos  \|f  sin  0  +  z  sin  f]^   _^  [(y-yo)cos  ®]^ 
(.n  -  -2-;  tn  -  -f- ; 


[-(y-yo)sin  -^   sin  @+  z  cos  i]/]^ 


(i5+"W 


1 


Designating 


Ai   =  cos  TJr  sin  ©  /  (n  -  -2~)  j 
Aa  =  cos  0/  (n  -  ^) , 
A3  =  -sin  T|r  sin  ®  /  (n  +  ~), 
A4.  =  sin  \Ir  /  (n  -  ^) , 
A5  =  cos  \lr  /  (n  +  ^) , 
and  substituting  yields  upon  rearrangement: 

(Ai2+A22+A32)(y-.yo)2  +  2(Ai  A4+A3A5)  (y-yo)z 

+  (A/^^+A5^)z^  =  1. 
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To  remove  the  cross-product  term,  set 
(y-yo)  =  -Ni  sin-^  +  N2  cos ''I', 
z  =  Ni  cos  ■^  +  N2  sin^, 
and  substitute  into  the  equation  just  given  to  obtain 

BtNi^  +  B2N1N2  +  B3N2^  =  1 
where 

Bi   =  (Ai2+A2^+A3^)  sin^-  2(Ai  A4.+A3A5)  sin^  cos-J^ 

+  (Aa^+Aj^)  cos^, 
B2   =  2(AiA4.+A3A5)  cos  2*" - 

(A12+A22+A32-A4.2-A52)  sin  2^, 
B3  =  (A12+A22+A32)  cos2^+  2(AtA4+A3A5)  sin^cos^ 
+  (A4^+A52)  siB^'i'. 
If  Ni  and  N2  are  to  be  identified  as  the  principal 
axes  of  the  ellipse,  then  it  is  necessary  that 

B2  =  0 
which  is  satisfied  if 


^.  1  tan-f-^i^^^^;^ 
^     2      1^^^  2+^2^+A32-A4.^-. 


2+A2^+A3^-A4.^-A5^. 
The  birefringence  and  the  orientation  angle  have  been 
defined  respectively  as  the  difference  between  the  lengths 
of  the  major  and  minor  axes  of  the  optical  ellipse  and  the 
angle  between  the  principal  flow  direction  (k)  and  the 
major  optic  axis.  From  these  definitions  it  is  clear  that 
the  effective  orientation  angle  is  identically  ^  as  just 
derived.  The  effective  birefringence  is  simply 


AN  =     ' 


n/E7  ./^ 


or 
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AN     =      [(Ai^+A2^+A3^)sin^  -   2(Ai  A4.+A3A5)sin^   COSTS' 
+   (A42+A52)cos2*]-V2 
-      [(Ai2+A2  2+A3  2)cos2fr  +  2(AiA4+A3A5)sin-<Jr    cos^ 
+   (A4.2+A52)sin2^]-V2    . 


APPENDIX  B 
PREPAHATION  OF  MILLING  YELLOW  SOLUTIONS 

Two  stock  solutions  of  Milling  Yellow  were  prepared, 
the  first  in  April,  1972,  and  the  second  in  March,  1973. 
The  procedures  were  the  same  in  principle  as  those  followed 
by  previous  investigators  such  as  Peebles,  Garber,  and  Jury 
(1953).  Changes  in  the  properties  of  these  solutions  with 
time  are  discussed  in  Appendix  H.  The  manner  in  which  the 
concentration  was  measured  is  discussed  in  Appendix  L, 
Original  Stock  Solution 

Seven  and  one-half  liters  of  tap  water  were  heated 
to  about  60°  C  in  a  galvanized  iron  pail  using  an  immer- 
sion heating  element.  To  the  warm  water  was  added  113 
grams  of  commercial  dye.  Milling  Yellow  NGS,  supplied  by 
Keystone  Aniline  and  Chemical  Company,  Chicago,  Illinois, 
The  mixture  was  stirred  until  all  of  the  dye  particles 
had  sunk  beneath  the  surface,  after  which  an  agitator  of 
the  egg-beater  type  was  introduced  and  heating  continued. 
The  temperature  was  raised  to  and  maintained  at  a  temper- 
ature just  below  the  boiling  point,  typically  93°  C.  No 
attempt  was  made  to  control  the  nominal  1.5  percent  concen- 
tration for  evaporation  losses. 
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Periodically  the  mixture  was  inspected.  The  prepar- 
ation was  considered  complete  when  three  criteria  were  met: 

Clarity.  —  Though  deeply  colored,  the  solution  was 
clear.  Specifically,  the  end  of  a  stirring  rod  could  be 
seen  5  cm  beneath  the  surface  within  the  opaque  container. 

Absence  of  residue.  —  A  sample  of  the  hot  solution, 
taken  from  the  bottom  of  the  pail  with  a  hollow  tube,  could 
be  released  at  the  surface  of  the  fluid  without  observing 
any  undissolved  dye. 

Birefringence.  —  A  small  sample,  rapidly  cooled, 
remained  clear  and  exhibited  birefringence  when  stirred 
or  swirled  between  crossed  polarizers. 

Having  met  these  criteria,  the  solution  was  added 
to  the  polystyrene  container  in  which  the  bulk  stock  was 
stored.  The  earlier  preparations  were  allowed  to  stand 
and  cool  for  a  sufficient  time  that  the  temperature  of 
the  bulk  stock  would  never  rise  above  45°  C  in  the  poly- 
styrene container. 

It  was  noticed  that  any  preparation  which  was  allowed 
to  stand  contained  a  small  amount  of  sedimentation  despite 
the  test  for  residue.  The  warm  solution  has  a  viscosity 
not  appreciably  different  from  water,  and  it  was  easy  to 
decant  the  clear  solution  and  include  the  residue  in  the 
bottom  of  the  pail  at  the  beginning  of  the  next  prepar- 
ation. The  sediment  which  remained  at  the  bottom  of  the 
pail  at  the  end  of  the  preparation  sequence  was  discarded. 


V7i 


Fresh  Stock  Solution 

Preparation  of  the  fresh  stock  solution  followed  the 
procedure  used  for  the  original  stock  solution  with  the 
following  changes: 

The  fresh  dye,  though  provided  by  the  same  supplier, 
was  from  a  different  lot# 

The  galvanized  pail  was  replaced  by  two  3-gallon, 
heavy-weight  aluminum  containers.  Both  were  heated  on 
hot  plates.  One  was  agitated  with  a  magnetic  stirrer. 
The  other  was  not.  No  difference  between  stirred  and 
unstirred  samples  was  evident. 

Distilled  water  was  used  in  place  of  tap  water. 

The  amount  of  Milling  Yellow  used  in  the  preparation 
of  the  separate  batches  was  not  controlled,  two  "heaping 
cupfuls"  being  used  each  time  for  each  container;  however, 
the  total  amount  of  Milling  Yellow  was  weighed  as  a  single 
quantity  before  preparation  of  the  fresh  stock  solution 
began,  and  the  total  mass  of  the  bulk  stock  was  determined 
when  preparation  was  complete.  On  this  basis  the  nominal 
dye  concentration  was : 

a  =    1351.8  gm  of  d7e  ^^^^^     ^     ^^^2  percent. 
83500  gm  of  solution 

Criteria  for  lack  of  residue  and  for  clarity  were 
retained.  The  criterion  for  birefringence  was  eliminated 
as  unnecessao^y. 
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Following  the  separate  preparation  of  the  individual 
batches,  the  entire  bulk  stock  was  heated  to  50°   C  and 
then  cooled  slowly  while  the  stock  was  circulated  between 
two  containers.  The  temperature  difference  between  the 
containers  was  never  greater  than  0.1  C°,  so  that  a 
corresponding  uniformity  of  concentration  was  inferred. 
The  50®  C  temperature  was  determined  by  the  heating  capa- 
city of  the  immersion  heating  elements  which  were  used; 
however,  this  temperature  exceeds  by  a  comfortable  margin 
the  transition  temperature  of  about  30**  C  below  which 
Milling  Yellow  becomes  birefringent. 


APPEITOIX  C 

VARIATION  IN  FLOV/  RATE 
AS  AMOUITT  OF  LIQUID  IN  OVERHEAD  TANK  DECREASES 


Referring  to  Figure  10,  and  recalling  that  the  pump 
is  not  operated  continuously,  it  will  be  evident  that  the 
amount  of  fluid  in  the  overhead  tank  decreases  during  each 
run  with  a  corresponding  reduction  in  the  pressure  head, 
A  test  was  conducted  to  determine  if  this  reduction  had 
a  significant  effect  upon  the  flow  rate  in  the  channel. 

The  flow  rate  used  during  the  test  was  about  six 
times  higher  than  any  used  during  later  testing  so  that 
the  magnitude  of  the  effect  was  increased  proportionately. 
The  procedure  was  as  follows. 

The  test  was  conducted  in  Channel  6,  Thirty  seconds 
after  the  control  valve  at  the  bottom  of  the  conduit  was 
opened,  four  consecutive  samples  of  the  efflux  were  taken 
for  timed  intervals  of  roughly  ten  seconds  each.  A 
photograph  of  the  birefringent  pattern  was  taken  before 
each  of  these  samplings.  At  the  end  of  this  sequence, 
about  two  minutes  in  all  after  the  valve  was  opened,  the 
flow  was  halted  and  the  efflux  for  each  of  the  timed 
intervals  was  weighed  and  its  temperature  recorded.  These 
data  are  given  in  Table  C-I, 
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TABLE  C-I 
CHANGE  IN  FLOW  RATE  AS  OVERHEAD  TANK  EMPTIES 


Order 
sampled 

Efflux, 
gm 

1831.9 

Time, 
sec 

10.1 

Flow  Rate, 
gm/sec 

181.4 

Temp, 

Order,  temp 
measurement 

1. 

21.2 

2 

2 

1699.8 

9.5 

178.9 

21.4 

1 

3 

2010.9 

11.1 

181.2 

21.2 

3 

4 

1110.1 

6.3 

176.2 

20.8 

4 

From  these  data  the  flow  rate  during  the  test  was 
G  =  179.4  ±  2.4  gm/sec.  The  standard  deviation  In  the 
flow  rate  is  only  slightly  higher  than  the  1  percent 
variation  to  be  expected  from  the  ±0.1-sec  precision 
of  the  timer.  The  photographs  taken  during  the  run  are 
indistinguishable.  The  temperature  variation  appears  to 
have  resulted  from  cooling  of  the  samples  which  remained 
on  the  floor  between  weighings.  The  floor  tempera t^lre 
is  noticeably  cooler  than  the  ambient  temperature  in  the 
vicinity  of  the  overhead  tank. 

On  the  basis  of  this  test  it  was  concluded  that 
variations  in  the  fluid  level  of  the  overhead  tank  have 
a  negligible  effect  upon  the  flow  rate  in  the  channels. 


APPENDIX  D 
ALIGNMENT  OF  POLARIZING  ARRAYS 

The  polarizing  material  used  in  this  dissertation 
was  commercial  Polaroid  sheeting  secured  to  the  inside 
surfaces  of  two  parallel  plates  of  transparent  plexiglass. 
One  of  the  two  longer  edges  of  each  plate  was  milled  to 
provide  a  reference  surface,  after  which  the  two  plates 
were  clamped  together  and  drilled.  When  the  clamp  was 
removed  threaded  rods  were  passed  through  the  holes  and 
provided  with  nuts  which  could  be  tightened  to  secure 
the  plates  at  any  desired  location  on  the  rectangular 
conduit  through  which  the  Milling  Yellow  was  observed  in 
flow.  The  individual  elements  of  the  polarizing  arrays 
shown  in  Figure  12  were  aligned  with  the  milled  edge  of 
the  plastic  plates  in  the  following  manner. 

On  a  piece  of  graph  paper  were  drawn  a  series  of 
30  mm  by  5  nmi  rectangles.  These  rectangles  were  oriented 
so  that  some  were  horizontal,  some  vertical,  some  were 
sloped  at  45°  angles  to  the  right,  and  some  at  4-5°  to  the 
left.  The  Polaroid  sheeting  was  laid  on  top  of  the  graph 
paper  and  secured  to  it  by  a  thin  coating  of  rubber  cement, 
the  optical  axis  of  the  sheeting  lying  parallel  with  one 
axis  of  the  graph  paper.  Rectangular  strips  were  then  cut 
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from  the  Polaroid  sheet  along  the  lines  already  drawn. 

One  of  the  plastic  plates  was  laid  on  a  second  sheet 
of  graph  paper  having  the  same  scale,  the  milled  edge  of 
the  plate  coinciding  with  one  axis  of  the  paper.  One  of 
the  rectangular  bits  of  sheeting  was  selected  by  observing 
the  orientation  of  its  optic  axis  by  means  of  the  graph 
paper  still  cemented  to  its  under  side.  This  rectangle 
was  slid  carefully  beneath  the  plate  and  maneuvered  until, 
looking  through  both  the  transparent  plastic  and  the  polar- 
izing material,  it  could  be  seen  that  the  grid  on  the 
rectangle  was  aligned  with  the  grid  on  the  imderlying 
sheet.  When  this  was  accomplished  the  small  strip  was 
taped  in  place  and  its  alignment  reconfirmed,  after  which 
the  graph-paper  backing  was  carefully  peeled  away.  The 
process  was  repeated  until  all  of  the  elements  shown  in 
Figure  12  had  been  secured  in  place.  For  strips  having 
45*'  orientations  it  was  necessary  to  align  the  milled  edge 
with  a  line  drawn  at  45"  on  the  underlying  graph  paper,  but 
this  posed  no  difficulties. 

After  the  plexiglass  plates  were  secured  to  the 
conduit,  the  alignment  was  completed  by  aligning  the 
milled  edges  with  each  other  outside  the  wall  of  the 
conduit  and  parallel  with  that  wall,  or,  alternately,  by 
aligning  each  milled  edge  in  turn  to  a  common  vertical 
reference,  usually  the  centerline  of  the  bolts  which  hold 
the  conduit  together. 


APPENDIX  E 
CALIBRATION  OF  THE  H(3PPLER  RHE0-VI3C0METER 

The  instr-ument  used  in  this  investigation  v/as  the 
"Rheo-Viscometer  nach  H5ppler,"  serial  number  002770, 
manufacturer:  VEB  Priifgerate-Vferk  Medingen,  Sitz  Freital, 
D.  D.  R.  The  serial  number  of  the  cylinder  v/as  31871. 

According  to  the  instruction  booklet  provided  by 
the  manufacturer,  the  viscosity  of  a  Nev/tonian  fluid  is 
obtained  by  use  of  the  simple  expression 

Pl<l     =     KPt 
XvThere  K  is  a  calibration  constant,  P  is  the  force  applied 
to  the  ball  divided  by  the  cross-sectional  area  of  the 
ball,  and  t  is  the  time  of  fall  for  a  displacement  of 
3.00  cm.  Although  both  cylinder  and  ball  have  dimensions 
which  are  specified  to  a  few  thousandths  of  a  centimeter, 
the  closeness  of  fit  requires  calibration  for  each  com.bin- 
ation  of  cylinder  and  ball.  Further,  it  is  necessary 
that  the  orientation  of  both  ball  and  cylinder  be  noted 
so  that  there  is  no  rotation  of  the  one  relative  to  the 
other  each  time  the  cylinder  is  removed  and  replaced. 

The  ball  and  cylinder  used  in  this  dissertation  were 
calibrated  on  May  3,  1973,  using  a  "normal  fluid"  supplied 
by  the  manufacturer.  The  viscosity  of  the  norm.al  fluid 
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was  certified  as  11.06  cP  at  20«»  C  if  used  within  three 
months  of  April  2,  19735  the  date  at  which  it  was  standard- 
ized. Temperature  control  at  20*»  C  was  provided  by  a 
Haake  Constant  Temperature  Circulator,  Model  F423,  serial 
71011.  A  partially  insulated  container  of  ice  water  was 
immersed  in  the  bath  of  the  Haake  circulator  during  the 
calibration  to  provide  ambient  cooling.  The  ice  water 
was  replaced  at  about  ten-minute  intervals. 

The  calibration  data  are  provided  in  Table  E-I.  The 
temperature  recorded  here  are  from  the  gas- jacketed, 
precision  thermometer  located  adjacent  to  the  cylinder  in 
the  constant  temperature  bath  of  the  viscometer  (see 
Figure  21) ,  From  these  calibration  data  a  separate  cali- 
bration constant  K  was  calculated  for  each  value  of  P,  the 
so-called  "average  stress"  on  the  sliding  ball.  These 
data  are  plotted  in  Figure  E-1.  It  will  be  noted  that  K 
varies  very  little  along  the  ordinate  of  the  graph. 

The  form  of  this  curve  suggests  the  presence  of  a 
frictional  resistance  in  the  viscometer  bearings.  This 
supposition  led  to  a  corrected  viscometer  equation  of 
the  form 

jun  =  K(P-Pf)t 

where  P^  is  the  frictional  resistance.  Rearrangement  of 
this  equation  into  the  form 
Pt  =  P^t  +  ;ijj/K 

shows  that  the  slope  d(Pt)/dt  would  have  a  value  P„  when 
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P,  gm/cra2  t,  sec   T,  °C 


TABLE  E-I 
CALIBRATION  DATA 

P,   gm/cm^     t,   sec       T,   *»C 


250 

8.4 

20.00 

40 

51.9 

20.00 

8.4 

20.00 

51.9 

20.00 

8.4 

20.00 

51.9 

20.00 

8.4 

20.00 

51.9 

20.00 

8.4 

20.00 

51.8 

20.02 

K250 

=  0.00527 

K40 

=  0.00533 

200 

10.5 

20.00 

25 

83.3 

20.02 

10.5 

20.00 

83.3 

20.00 

10.4 

20.00 

83.4 

19.96 

10.4 

20.00 

83.4 

19.94 

10.5 

20.00 

83.5 
83.3 

20.02 
20.02 

K200 

=  0.00529 

83.3 

20.00 

160 

13.0 
13.1 

20.00 
20.00 

K25 

=  0.00531 

13.0 

20.00 

20 

104.5 

20.00 

13.0 

20.00 

104.8 

20.00 

13.0 

20.00 

104.8 
104.5 

20.02 
20.04 

Ki  60 

=  0.00531 

104.9 

20.00 

100 

20.9 
20.8 

20.00 
20.00 

K20 

=  0.00528 

20.8 

20.00 

15 

139.5 

20.00 

20.7 

20.00 

139.9 

20.00 

20.7 

20.00 

139.7 
139.7 

20.04 
20.02 

Ki  00 

=  0.00532 

140.4 

19.98 

80 

25.9 
26.0 

20.00 
20.00 

K15 

=   0.00527 

26.0 

20.00 

10 

211.6 

19.96 

26.0 

20.00 

212.2 

20.00 

25.8 

20.00 

212.8 
212.0 

20.00 
20.00 

K80  = 

=  0.00533 

212.2 

20.00 

50 

41.5 
41.6 

20.00 
20.00 

Ki  0 

=  0.00521 

41.5 

20.00 

5 

434.9 

20.04 

41.6 

20.00 

435.8 

20.02 

41.6 

19.98 

435.9 

20.00 

K50  =  0.00532 


K5  =  0.00508 
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the  product  Pt  was  plotted  versus  t.  When  this  was  done 
(not  shown)  the  result  was  P^  =  0.28  gm/cm^.  With  this 
correction  the  calibration  constant  could  be  expressed 
as 

K  =  )ij^/(P-0.28)t. 

Values  of  K  calculated  in  this  manner  are  plotted 
at  the  top  of  Figure  E-1  and  show  a  nearly  linear  vari- 
ation in  K  with  P.  For  the  line  drawn  in  Figure  E-1, 

K  =  ;pj^/(P-0.28)t  =  0.00538  -  O.OOOOOO388P. 
For  most  applications  it  may  be  stated  with  sufficient 
accuracy  that 

K  =  0.00534  cP-cmVgm-sec, 
the  units  to  be  used  for  >ijj  in  cP,  P  in  gm/cm^,  and  t 
in  seconds.  In  consistent  units  (;i„  in  poise,  P  in  dyne/cm^ 
t  in  seconds)  K  =  5.44  x  lO"^. 

As  a  check  of  the  calibration,  a  distilled  water 
sample  {yi^  -   1.002  cP  at  20°  C,  according  to  Swindells, 
Coe,  and  Godfrey,  1966)  was  tested  in  the  instrument. 
Because  of  a  restriction  on  the  use  of  the  viscometer  at 
lower  viscosities  (see  Appendix  J),  only  three  points 
could  be  obtained,  but  these  lay  in  the  range  (P  :$  15 
gm/cm^)  where  the  frictional  correction  is  most  signif- 
icant. As  Table  E-II  and  Figure  E-1  confirm,  the  cali- 
bration previously  performed  is  also  valid  for  distilled 
water . 
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TABLE  E-II 
CALIBRATION  CHECK  (WATER) 


P,  gin/cm2 

t,  sec 

T,  °C 

15 

12.6 

20.02 

12.7 

20.02 

12.7 

20.00 

12.7 

20.00 

12.7 

20.02 

Ki  5  =  0.00527 
10 


19.3 

20.02 

19.2 

20.02 

19.3 

20.04 

19.2 

20.04 

19.2 

20.04 

19.2 

20.02 

19.3 

20.02 

19.2 

20.00 

19.2 

20.00 

19.2 

20.00 

K10  =  0.00521 


39.7 

20.00 

39.5 

20.00 

39.6 

20.00 

39.5 

20.00 

39.5 

20.00 

K5  =  0.005065 
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Analysis  of  the  viscometer  in  Chapter  Five  resulted 
in  the  following  equation  for  the  viscosity  of  a  Nev/tonian 
fluid: 


56m^  ^Ji^l^ 


^N  -  95ilL  ^* 


where  6  is  the  maximum  width  of  the  eccentric  annulus, 

m 

a  is  the  radius  of  the  ball,  ft   is  the  distance  moved  by 
the  ball  in  time  t,  and  L  is  the  effective  length  to  be 
obtained  by  calibration.  Comparison  of  this  result  with 
the  calibration  equation  yields 

L  =  5&m^/96alK 
and  on  substitution, 

L  =  (5)(.007)V(96)(.798)(3.00)(5.44  X  10"^) 

=  0.1365  cm. 
Following  the  correction  for  friction,  the  variation 
in  the  calibration  constant  with  P  is  less  than  1  percent. 
A  much  more  serious  uncertainty  is  that  due  to  b^   for 
^^^hich  the  value  given  above  was  based  upon  the  subtraction 
of  a  micrometer-read  diameter  for  the  ball  from  manu- 
facturer 's  specification  for  the  inside  diameter  of  the 
cylinder : 

6jjj  =  1.603  -  1.596  =  0.007  cm. 
This  value  is  accurate  to  only  one  significant  figure.  To 
illustrate  the  effect  of  this  uncertainty,  consider  the 
expression  for  the  calibration  constant  proposed  by 
Hubbard  and  Brown  (1943a)  for  the  rolling  ball  viscometer: 
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where  Jj^  was  calculated  by  Hubbard  and  Brown  to  be  0,398 
and  K  has  been  corrected  to  compensate  for  the  fact  that 
the  force  in  the  Hubbard-Brown  analysis  was  due  to  the 
weight  of  the  ball,  while  in  the  present  case  a  force 
was  applied  along  a  rigid  shaft.  This  equation  may  be 
rearranged  to  solve  for  6^^  when,  as  in  the  present  case, 
the  value  of  K  is  known  from  calibration.  When  this  is . 
done  it  is  found  that 

6jjj  =  0.00552  cm, 
in  which  case  the  corresponding  effective  length  becomes 
L  =  0.0670  cm. 

Since  the  Hubbard-Brown  analysis  was  for  the  rolling 
ball  instrument  a  difference  in  6^^  is  not  surprising,  but 
the  variation  in  L  with  bj^   is  dramatized  by  the  comparison. 

In  the  present  study  the  values  &jjj  =  O.OO7  and 
L  =  0,1365  were  accepted  as  the  most  likely  ones. 


APPENDIX  F 
BILINEAR  MATERIALS 

An  important  class  of  empirical  constitutive  equations 
are  those  which  postulate  two  separate  relationships,  often 
linear,  with  a  critical  value  of  the  shear  stress  at  which 
the  transition  occurs.  Figure  F-1  is  an  example. 

The  integration  of  the  viscometer  equation  in  this 
case  requires  the  separate  treatment  of  two  regions  of 
flow,  one  in  which  all  the  stresses  are  below  the  critical 
value,  Tq,  and  the  other  in  which  the  transition  occurs 
within  the  annulus.  As  shown  in  Figure  F-2,  corresponding 
to  each  value  of  t^  there  will  exist  an  angle  9^  repre- 
senting the  maximum  value  of  0  for  which  a  single  region 
of  flov;  may  be  assumed. 

In  this  case  the  equation  for  the  viscometer  is  no 
longer  equation  (5»1)>  but  its  analogue: 

4L2 


1/t  = 


TTaJ?P2 


\^'  r^  gi  (t)  dT  dG  + 

-'o  Jo 


rf't^  gi(T)  dT  de  +  ^''f  T  g2(T)  dr  de 


where 

g(T)  = 


137 
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FIGURE  F-1.  Bi- functional  material. 


FIGURE  F-2,  Flow  field  within  viscometer  showing  boundary 
between  regions  obeying  separate  constitutive  equations. 
Along  the  boundary,  t  =  t^. 


189 


and  0  is  related  to  t  by 


Of  special  importance  is  the  bilinear  case: 


g  = 


r  '^/-^f 


T  <  T. 


In  this  instance, 


i/t  = 


TTaSP2  L^o  Jo  P 


TT  f-Tw 


n    ~  d-r  de     +    f    1    (tJ-  -  kT)  dt  dG 


Integrating  once  and  substituting 
-c^     =     P6ni(l-cos  e)/4L 

yields 


Vt    = 


II^S>— «3de   . 


4L^ 

TraJ?P2  Ll9^o'L^  J^ 


P3&m3    r. 


Jfi  3uo  192uoo  L3  Jft_ 


\ 


\  3uoo 


de      +     ^~|^rCl-COS    0)2    d0      + 
^2    h^    Jfi 


Integrating  again  and  collecting  terms,  a  final  expression 
for  1/t  is  obtained.  The  result  is  given  on  the  next  page. 
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^     P6a!£o 


56, 


8  sin  Gq     + 


3  sin  29f 


+     2_sin3e£    +^I„     +     (1-cos  0c)2(^-ej 
3  cro)icD  ^  ° 


l6TraJl 


3(w-ep)     +     A-  sin  Gc     - 


sin  2Gc 
2 


-     k(l-cos  Gjj)2(tt-G(,) 


The  substitution 

Co  =  (1/jiio)  -  (l/;ioo ) 
has  been  made  to  reduce  slightly  the  number  of  terms, 
and  where  t^  appeared  it  was  replaced  by  6^^(1-003  G(,)/4L. 


APPENDIX  G 
SOME  CHANNEL  CONSTANTS 


262(5i  +  52)Sop/6i       A2/4-6i2 

Channel     61,   cm     62?   cm     So,   cm2  gm  (x  108) 


1 

.672 

.663 

.127 

.336 

.192* 

2 

.672 

.620 

.126 

.303 

.192* 

6 

.671 

.596 

.126 

.28^ 

.193* 

*A     =     5892  Angstrom 


191 


APPENDIX  H 
VARIATION  IN  MILLING  YELLOW  SOLUTIONS  WITH  TIME 

As  alluded  to  in  Chapter  Two,  other  investigators, 
beginning  with  Prados  (1957)?  have  reported  difficulty 
with  the  stability  of  Milling  Yellow  in  solution.  The 
experience  of  the  present  study  confirms  these  earlier 
findings. 

Two  factors  contribute  especially  to  the  variability 
of  Milling  Yellow  solutions:  evaporation  and  sedimentation, 

If  evaporation  is  not  rigidly  controlled,  the  concen- 
tration of  the  dye  increases  with  a  disproportionate  rise 
in  both  birefringence  and  intrinsic  viscosity.  This  was 
evident  in  the  current  work  as  shown  in  Figure  l8,  where 
the  birefringence  rises  markedly  between  runs  4l8  and  4-19 
due  to  the  accidental  displacement  for  about  ten  hours  of 
the  lid  on  the  lower  storage  tank.  Evaporation  also 
occurred  during  the  testing  when  the  lid  was  set  aside  to 
permit  the  efflux  from  the  conduits  to  be  caught  and 
measured.  A  deliberate  concentration  by  evaporation  was 
carried  out  on  the  original  stock  solution  in  January  1973 
just  prior  to  initial  testing. 

A  more  insidious  change  in  Milling  Yellow  properties 
occurs  as  the  result  of  slow  but  continuous  sedimentation 
to  which  the  organic  dye,  with  its  limited  solubility, 
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contributes  the  chief  constituent.  This  is  particularly- 
true  when  the  solution  is  diluted  to  concentrations  near 
the  limit  of  birefringence.  The  original  stock  solution 
set  in  a  tightly  sealed  container  from  May  to  December, 
1972.  During  this  time,  as  mentioned  in  Chapter  Two,  a 
continuous  process  of  evaporation  and  recondensation 
occiured  beneath  the  lid.  The  result  is  shown  in  Figure 
H-1.  The  apparent  viscosity  as  measured  with  the  Hoppler 
Rheo-Viscometer  according  to  the  manufacturer's  - 
instructions  ( i . e . ^  as  opposed  to  the  method  of  analysis 
in  Chapter  Five)  is  plotted  as  a  function  of  the  recip- 
rocal of  the  average  shear  stress  at  four  different  times 
during  1972.  The  reduction  in  viscosity  with  time, 
particularly  at  lower  shear  rates,  is  very  evident. 

During  the  testing  of  both  the  original  and  the 
freshly  made  stock  solutions  described  in  Appendix  B 
there  was  evaporation  and  sedimentation  occurring,  the 
former  having  the  dominant  role.  The  results  in  Chapters 
Four  and  Five  confirm  that  the  optical  and  Theological 
properties  of  Milling  Yellow  vary  significantly  with  time 
and  must  be  measured  at  the  time  of  use  if  quantitative 
results  are  to  be  obtained. 
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APPENDIX  I 

VARIATION  IN  MILLING  YELLOW  APPARENT  VISCOSITY 
WITH  CONCENTRATION  AND  TEMPERATURE 

When  measured  with  the  Hoppler  Rheo-Viscometer 
using  the  manufacturer's  instructions,  the  apparent 
viscosity  of  Milling  Yellow  varies  considerably  with 
the  concentration  of  the  fluid,  with  the  temperature, 
and  with  the  average  shear  stress.  Prior  to  the 
inception  of  the  testing  reported  in  Chapters  Four 
and  Five  of  this  dissertation  and  before  the  analysis 
of  the  viscometer  in  the  latter  chapter  was  completed, 
a  series  of  measurements  were  made  to  determine  the 
temperature  dependence  of  Milling  Yellow  solutions  of 
different  concentrations.  These  measurements  are 
described  in  this  appendix. 

Apparatus 

The  instrument  used  was  the  Hoppler  Rheo-Viscometer 
which  has  been  adequately  described  in  Chapter  Five  and 
Appendix  Ej  however,  the  calibration  described  in  Appendix 
E  had  not  been  performed  at  the  time  these  temperature 
studies  were  conducted  and  the  apparent  viscosities 
reported  herein  are  based  on  an  older  calibration  using 
the  simple  expression: 
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ji^    =  0.00543  Pt 
where  ^^  is  the  apparent  viscosity  in  cP,  P  is  the 
"average  shear  stress"  (force  on  ball  divided  by  cross- 
sectional  area  of  ball)  in  gm/cm^  [sis] ,  and  t  is  the 
time  required  for  the  ball  to  fall  30  mm. 

With  a  single  exception  (run  103-20  at  19**  C,  P  =  15) 
the  apparent  viscosity  was  measured  at  least  three  and 
oftener  five  times  for  each  sample,  each  temperature,  and 
each  value  of  P  indicated  in  the  data  which  follow.  The 
variations  in  the  fall  times  during  these  runs  were  compar- 
able to  those  in  Tables  XIX  to  XXVII,  although  the  temper- 
ature variations  in  the  earlier  runs  were  somewhat  larger. 
In  run  103-20  at  16"  C,  for  example,  the  measured  temper- 
ature ranged  from  15.70*'  to  15.88  °   C.  This  range  may  be 
regarded  as  typical  for  the  earlier  runs  as  may  the 
departure  of  the  actual  temperatures  from  the  nominal 
value.  Toward  the  end  of  the  testing  period  the  temper- 
ature variations  were  no  greater  than  those  reported  in 
Chapter  Five.  All  of  the  tests  were  conducted  in  October, 
1972,  Run  1012,  for  instance,  was  prepared  on  October  12. 

Preparation  of  Samples 

With  the  exception  of  one  sample  taken  from  the 
original  stock  solution  whose  preparation  has  been  detailed 
in  Appendix  B,  all  of  the  samples  were  prepared  by  weighing 
out  a  prescribed  amount  of  Milling  Yellow  and  adding  it  to 
a  sufficient  abundance  of  distilled  water  so  that  the 
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initial  concentration  was  less  than  1.3  percent.  The 
mixture  was  then  boiled  until  the  total  volume  was 
slightly  greater  than  a  predetermined  value.  At  this 
time  the  solution  v;as  removed  from  the  hot  plate  and 
placed  on  one  scale  of  a  triple  beam  balance  previously 
set  to  the  exact  mass  desired.  When  the  mass  was  attained 
by  evaporation,  the  solution  was  immediately  bottled  and 
capped.  A  check  of  the  solid  content  of  these  samples 
was  later  made  by  evaporating  to  dryness  as  described  in 
Appendix  L,  The  value  of  the  concentration  cited  in  the 
data  which  follow  is  the  value  given  in  Appendix  L, 

Experimental  Data 

The  variation  in  apparent  viscosity  with  temperature 
is  given  in  Figures  I-l  through  1-6.   It  will  be  noted 
that  there  are  no  data  on  certain  of  these  plots  above 
a  certain  temperature.  This  reflects  a  restriction  on 
the  use  of  the  viscometer  which  is  discussed  at  greater 
length  in  Appendix  J. 

Computation  of  Temperature  Coefficients 

The  replication  of  the  data  from  separately  prepared 
samples  is  quite  good  at  higher  concentrations  (q  >  2.3 
percent)  as  shown  in  Figures  I-l  and  1-2.  At  very  low 
concentrations  (q  =  1.3  percent)  the  agreement  between 
samples  104-14  and  1016-14  was  so  close  that  the  data 
coincide,  and  only  the  former  are  plotted  in  Figure  I-5« 
At  intermediate  concentrations,  which  unluckily  include 
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50 
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100 
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200 
^200        - 

©]Vcm2 


Temperature,  °C 

FIGURE  I-l.  Temperature  variation  of  samples  103-30, 
concentration  q  =  2.94-  (solid  lines),  and  1012-30,  concen- 
tration q  =  2.04  (dashed  lines).  The  load  on  the  ball  is 
given  in  the  right  hand  margin. 


Apparent 
viscosity, 

cP 


^^   25 


^    35 


60 
65 

*^  90 

gm/cm2 


Temperatiire ,   ®C 

FIGURE  1-2.  Temperature  variation  of  samples  107-25, 
concentration  q  =  2,33  (solid  lines),  and  1011-25,  concen- 
tration q  =  2.34  (dashed  lines).  The  load  on  the  ball  is 
given  in  the  margin.  Note  that  for  the  third  pair  of  lines 
the  load  on  the  ball  differs  in  the  t;vo  samples. 
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FIGURE  1-3.  Temperature  variation  of  samples  103-20, 
concentration  q  =  1,89  (solid  lines),  and  109-20,  concen- 
tration a  =  1,92  (dashed  lines).  The  load  on  the  ball  is 
given  in  the  right  hand  margin. 


4  - 


Apparent 
viscosity, 

cP        2 


Temperature,  ''C 

FIGURE  1-4.  Temperature  variation  of  samples  104-17, 
concentration  q  =  1.59  (solid  lines),  and  1011-17,  concen- 
tration q  =  1.57  (dashed  lines).  Load  on  ball:  e  10  gm/mS; 
O  15  gm/cm2;  ©  25  gm/cm^;  (D  35  gm/cm2. 
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Apparent 
viscosity, 

cP     ' 


Temperature,  ''C 


FIGURE  1-5.  Temperature  variation  of  sample  104-14-, 
concentration;  a'=  1.31.  Load  on  ball,  P:  •  5  gm/cm^; 
©  10  gm/cm2;  O  15^  gm/cm2. 


Apparent 
viscosity, 

cP 


Temperature,  °C 


FIGURE  1-6.  Temperature  variation  of  original  stock 
solution,  nominal  concentration:  q  =  1.5«   Load  on  ball, 
P:  O  15  gin/cm2;  ©  25  gm/cm^;  •  40  gm/cm^. 
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the  range  of  concentrations  studied  in  this  dissertation, 
the  results  obtained  from  separately  prepared  samples 
were  disparate,  as  shown  in  Figures  1-3  and  I-4-.  The 
disparity  was  also  evident  in  the  appearance  of  the 
sediment  resulting  when  the  samples  were  evaporated  to 
dryness  to  determine  their  concentration.  The  data  in 
Figures  1-3  and  1-4  form  the  principal  basis  for  the 
decision  not  to  attempt  the  preparation  of  standardized 
solutions  of  the  commercial  dye. 

Besides  a  failure  to  replicate,  there  is  evidence 
of  an  alternation  in  the  character  of  the  temperature 
dependence  as  the  concentration  decreases.  At  high 
concentrations  the  data  plot  as  parallel  lines  indicating 
that  the  temperature  dependence  is  independent  of  the 
shear  stress.  At  intermediate  concentrations  the  vari- 
ations due  to  temperature  increase  in  magnitude  as  the 
shearing  stress  decreases.  Although  straight  lines  have 
been  drawn  through  the  data  in  Figure  1-2  to  facilitate 
comparisons  from  figure  to  figure,  close  inspection 
indicates  a  distinct  curvature  is  present.  Since  each 
value  of  the  apparent  viscosity  reflects  a  wide  range 
of  local  shear  stresses,  the  curvature  may  be  due  to  the 
integration  of  two  characteristic  responses,  only  one  of 
which  is  dependent  upon  the  shear  stress.  In  Figure  1-2 
these  shear  rate  dependencies  are  nearly  negligible,  but 
at  the  lower  concentrations  plotted  in  Figures  1-3  and  1-4 
the  shear-rate  dependent  terms  dominate  the  relationship. 
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At  the  lov/est  concentration,  plotted  in  Figure  1-5)  the 
shear-rate  dependence  is  still  present,  but  its  magnitude 
is  negligible. 

The  straight  lines  dra^-m  in  Figures  I-l  through  1-6 

have  the  general  form 

-kT 
Pa  =  Pao  e 

from  which 

where  k  has  an  empirical  value  which  depends  upon  both 

the  concentration  and  the  average  shear  stress,  P.  Typical 

values  of  k  are  given  in  Table  I- I, 

The  data  for  the  original  stock  solution  v/hich  are 
plotted  in  Figure  1-6  are  accurately  represented  by  the 
eiapirical  expression 

p^     =  exp  [0.252(30. 3-T)P" -377 
from  vvhich 

For  comparison,  run  104-17  yields  hy  the  method  of 
least  squares: 

yi^     =     e:q3  [0. 335(26. 8-T)P" -^2^] 

and 

-  .4-29 
k  =  0.335  P 

The  temperatures  in  these  empirical  expressions  are  in 

degrees  Celsius, 
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TA3T,-F, 

I- 1 

TEMPER. 

A.TURH   ( 

COEFFICIENTS 

Run 

Concentrat 

\'it  % 

ion, 

she 

Average 
iaring  stress, 

P,    gryca^ 

^  ^    1     ^ua 
>^a  6T 

103-30 

2.94 

200 

100 

50 

35 

0.040 

0.039 
0.040 

0.044 

107-25 

2.33 

90 
65 
35 
25 

0.051 
0.050 

0.050 
0.050 

103-20 

1.89 

50 

35 
20 

15 

0.058 
0.061 
0.076 
0.035 

104-17 

1.59 

35 

25 

15 
10 

0.071 

0.086 

0.109 
0.121 

104-14 

1.31 

15 
10 

5 

0.026 
0.025 

0.024 

<CU' 


It  appears  that  at  higher  concentrations  (q>2.3 
percent)  and  at  low  concentrations  (q  ^^  1,3  percent)  the 
variations  in  the  apparent  viscosity  can  be  computed  from 
the  empirical  constants  in  Table  I- I.   At  intermediate 
concentrations  the  temperature  dependence  has  the  form 

p^   5T 

but  the  constants  ko  and  n  must  be  obtained  on  a  case  by 
case  basis. 


APPENDIX  J 

RESTRICTION  OF  VISGOI-IETER  TO  LIQUIDS  liAVDIG- 
VISCOSITY  ABOVE  4-  CENT  IPO  ISE 


The  instrument  booklet  furnished  v/ith  the  Hoppler 
Rheo-Visconeter  gives  the  useful  range  of  the  instrument 
as  4  to  280  centipoise  for  the  ball  and  cylinder  which 
were  used  in  this  dissertation.  V/hen  queried  about  the 
use  of  the  instrument  at  lower  viscosities,  the  manu- 
facturer replied: 

Ue  are  sorry  to  inform  you  that  the  extension 
of  the  application  range  of  the  Rheo-Viscometer  to 
liquids  with  a  viscosity  below  4  cP  is  not  possible 
in  general.  The  measurement  of  such  viscosities  may 
be  carried  out  only  v;ith  Newtonian  liquids.* 

The  difficulty  v;hich  arises  when  the  viscometer  is 

used  indiscriminately  with  liquids  having  viscosities 

less  than  4  cP  is  sho\m  in  Figure  J-1.  Here  the  apparent 

viscosity  of  a  Killing  Yellow  sample  at  21.8°  C  has  been 

plotted  against  the  reciprocal  of  the  so-called  average 

shear  stress,  P.   It  is  evident  that  the  character  of  the 

flow  undergoes  an  abrupt  transition  on  or  about  P  =  60 

gm/cm2.  A  similar  transition  is  observed  in  water  at  aboTit 

P  =  15  gm/cm2. 


*i'Ieyer,  Conuaercial  Director,  VEB  MLW  Priifgerate-VJerk, 
Medingen,  Sits  Freital,  D.  D.  R,  Letter  dated  October  19, 
1972 
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FIGURE  J-1,  Transition  in  measurement  of  apparent 
viscosity  when  P  is  too  large.  Fluid  sampled  from  orig- 
inal stock  solution  in  May,  1972.  Temperature:  21.8°  C. 
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FIG'jP.E  J-2.  Transition  in  measurement  of  apparent 
viscosity  when  temDerature  is  too  high,  sanDle  104—14. 
Load  on  ball:  O  10  gm/cm2;  ©  15  gm/cm2;  •    20  gm/cm2, 
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Similar  transitions  are  observed  if  P  is  held 
constant  and  the  temperature  of  the  fluid  is  raised.  As 
shovm  in  Figure  J-2,  above  a  certain  critical  temperature 
the  fall  time  of  ball  (and  hence  the  apparent  viscosity) 
remains  nearly  constant. 

All  of  the  data  observed  during  the  current  study, 
for  both  Newtonian  and  non-Nev;tonian  liquids,  indicate 
that  the  transition  in  the  viscometer  occurs  at  about 

P/^a^  =  15.5  gm/cm2-cP2. 
If  P/;ia^  is  significantly  greater  than  15.5,  reliable 
results  cannot  be  obtained  with  the  Hoppler  instrument. 
Since  the  largest  mass  supplied  by  the  manufacturer  (i.e., 
the  sum  of  the  masses  provided)  corresponds  to  P  =  200 
gm/cm2,  the  lowest  viscosity  v;hich  can  be  measured  without 
concern  for  the  transition  is 

;^a  =  v^200/l5.5  =  3.6  cP, 
a  value  in  good  agreement  with  the  manufacturer's 
restrictioni  however,  if  only  those  masses  less  than 
^^•^  ?a.^   are  used  v^ith  the  less  viscous  liquids,  the  diffi- 
culty v;ith  the  transition  is  circumvented,  .^ji  example  is 
the  use  of  distilled  water,  as  described  in  Appendl:c  E,  to 
check  the  calibration  of  the  instrument. 

It  is  not  certain  what  phenomenon  is  responsible  for 
the  change  in  the  response  of  the  viscometer  at  the 
transition  point.  Three  possibilities  are  turbulence, 
flow  separation,  and  movement  of  the  ball  away  from  the 
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wall  of  the  cylinder.  Since  no  "use  of  the  viscometer  has 
been  made  for  values  of  F/pa^  greater  than  l5.5,  it  is 
not  necessary  to  discr ruminate  between  these  alternatives; 
however,  the  question  is  not  without  interest  and  a  brief 
discussion  of  each  possibility  follows. 

Turbulence .  -  It  is  commonly  knoim  (e.g.  Schlichting, 
1968,  pp.  15-19)  that  drag  becomes  substantially  greater 
during  the  transition  from  laminar  to  turbulent  flow.  The 
Reynolds  number  for  the  viscometer  ball  is 

Re  =  2pJia/)it 
where  p  is  the  density,  i/t  is  the  velocity  of  the  ball, 
a  is  the  radius  of  the  ball,  and  ji   is  the  viscosity  of 
the  fluid.  For  water,  Re  =  38O  when  transition  occurs  in 
■  the  viscometer.  Calculation  of  the  ReiTiolds  number  based 
^lpon  the  average  width  of  the  annulus  and  the  average  flow 
velocity  in  the  annulus  yields  the  value  Re  =  95.  Soth  of 
these  estimates  are  less  than  those  usually  associated 
with  turbulence,  although  Eckert  and  Irvine  (1956)  have 
shown  that  in  narrow  channels  of  varying  width,  turbulence 
can  occur  in  the  wider  portion  even  when  flow  in  the 
narrovrer  part  of  the  cross-section  remair-S  laminar. 
Occurrence  of  the  Eckert- Irvine  phenomenon,  iiowever,  would 
be  expected  to  yield  a  gradual,  rather  than  an  abrupt 
transition,  so  this  possibility  appears  unlikely. 

A  Reynolds-number  dependence  is  evident  in  the 
transition  whether  turbulence  is  present  or  not.  Since 
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the  apparent  viscosity  is  defined  by 

>^a  =     ^^-» 
it  follows  that 

t  =     >i^/!T. 

Substituting  this  value  for  t  into  the  expression  given 

previously  for  the  Reynolds  nunber  yields 
Re  =  2Kp^P/a;^g^2 

V7hich  differs  only  by  a  proportionality  constant  from 
the  criterion  for  transition: 
P/j^g2  =  constant. 
Flow  Separation.  -  Separation,  if  it  occurs,  must 
take  place  on  the  dovmstream  side  of  the  ball,  that  is, 
on  its  upper  surface.  .'In  analysis  of  this  possibility 
requires  specification  of  flow  around  the  sphere  in  the 
viscometer .   In  this  dissertation  the  flow  has  been 
specified  only  along  the  "effective  length"  of  the  annulus, 
a  narrow  region  in  the  vicinity  of  the  point  of  closest 
clearance  betv;een  ball  and  cylinder.   If  the  flow  along 
the  effective  length  is  compared  with  flow  in  a  convergent- 
divergent  channel,  then  the  criterion  for  flow  in  the 
latter  case  may  give  some  indication  of  the  lii:elihood  of 
flow  seiparaticn  in  the  viscometer. 

Millsaps  and  Pohlhausen  (1953)  found  that  backflow 
occurs  in  divergent  charjnels  having  5°  half-angles  vrhen 
the  Reynolds  number  exceeds  684.   In  this  c-se  the  Reynolds 
number  was  defijied  as: 
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where  u^^^  is  the  speed  at  the  center  of  the  divergent 

channel  and  r^^  is  the  radial  distance  of  the  flov;  from 
the  origin  of  the  divergence.   In  the  Hoppler  viscometer 
the  na:>:iaum  flow  rate  for  a  Nev;tonian  fluid  is  about 
2,4  times  the  average  rate,  so  one  can  estimate 
^ax  =  2.4II  =  2.4  (^/t)(a/6j 

ia 


=   2.4 


t6 


'm 


Perhaps  the  least  unreasonable  estimate  of  r^  is  that  it 

c 

is  half  the  previously  calculated  effective  length: 

^c  =  ^/2- 

With  these  very  crude  estimates  the  corresponding  Resmolds 
number  may  be  calculated.  When  the  viscomieter  undergoes 
its  transition  v;ith  vra.ter  (p  =  1  gm/cm^,  yx   =  0,01  poise, 
t  =  12,7  sec)  the  corresponding  Reynolds  number  is 
Re  ^  500. 
Since  separation  is  more  apt  to  occur  at  larger  half- 
angles  than  at  sm.aller  ones,  and  since  the  half-ajngle*  in 
the  Rheo-Viscometer  at  the  upper  end  of  the  effective 
length  is  about  2,5°)  compared  to  the  5°  for  which  Re  =  684 
implied  separation,  the  present  analogy,  if  valid,  seems 
to  indication  that  flow  separa.tion  is  not  responsible  for 
the  transition  point  observed  in  the  viscometer. 


*In  the  Hoppler  Rheo-Viscometer  the  half-angle  is  defined 
as  half  the  angle  between  the  v;all  and  a  tangent  to  the  ball 
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Movement  av;ay  from  the  wall.  —  Irving  (1972)  has 
sho\m  that  cylinders  in  inclined  tubes  fall  as  much  as 
twice  as  fast  as  those  in  vertical  tubes.  This  result 
leads  to  the  expectation  that  movement  of  the  ball  away 
from  the  v/all  in  the  Hoppler  instrument  would  increase 
the  fall  time  vrith  a  corresponding  increase  in  the 
apparent  viscosity.  In  the  original  design  of  the  Rheo- 
Viscometer,  Hoppler  (194-2)  insured  that  the  ball  would 
remain  against  the  wall  by  inclining  the  cylinder  to  a 
significant  angle  from  the  vertical.   In  the  modern 
design  a  coimterweight  attached  to  the  side  of  the  ball's 
shaft  accomplishes  the  same  objective.  Mien  the  pressure 
drop  in  the  eccentric  annulus  exceeds  a  certain  value, 
the  ball  may  be  expected  to  swing  av/ay  from  the  wall  into 
the  flow  much  as  cylindrical  shell  viscometers  center 
themselves  as  shown  by  Lindgren  (1972).   It  is  not  unlikely 
that  this  is  what  is  happening  at  the  transition  point  in 
the  Hoppler  viscometer,  but  it  cannot  be  stated  v;ith 
certainty  that  this  is  the  case. 


appeI'Idd:  k 
viscometer  response  at  very  slow  fall  times 

Among  the  assumptions  implicit  in  the   analysis  of 
the  Hopplor  Rheo-Visconeter  in  Chapter  Five  vras  the 
presLiiiption  that  the  response  of  the  viscometer  was 
independent  of  the  location  of  the  ball  in  the  axial 
direction.  This  assujnption  vras  tested  by  observing 
the  descent  of  the  ball  at  speeds  sufficiently  slci;  so 
that  the  fall  time  could  be  separately  measured  at 
each  l-rjr.  interval  along  the  path.   As  an  additional 
check,  the  point  at  which  the  fall  was  initiated  was 
also  varied. 

It  was  found  that  fall  tim.es  for  l-ivm.   increments 
vary  from  point  to  point  along  the  cylinder  axis.  The 
speed  ap:oeared  greater  (times  were  smaller)  during  the 
fir.st  millimeter  of  travel.  As  a  result,  it  appears 
that  calibrations  of  the  Hoppler  Rheo-Viscomieter  are 
valid  only  if  the  starting  pnint  and  fall  distance  for 
each  measurement  are  identical  with  tliose  used  for  cali- 
bration of  the  instrur-'icnt.  Further  details  are  provided 
in  th.e  paragraphs  which  follov;. 

Tests  were  conducted  at  21.3°  C  using  a  1.89  percent 
solution  of  ICilling  Yellow  and  an  "average  shearing  stress" 
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of  P  =  10  gn/cEi^.  Under  those  conditions  the  total  fall 
tine  is  about  340  seconds,  so  tiiat  the  average  fall  time 
for  each  of  the  30  1-ran  intervals  along  the  axis  exceeds 

11  seconds,  an  ample  tine  for  sea_uential  measurements. 

The  procedure  was  as  folloirs.  Starting  points  at 

12  different  locations  within  the  cylinder  were  selected 
\/ith  the  aid  of  a  table  of  random  numbers.  The  ball  was 
positioned  opposite  the  first  starting  point  and  the 
temperature  was  recorded.   Simultaneously,  the  ball  was 
released  and  a  timer  started.  As  the  ball  passed  the 
first  millimeter  marking,  easily  observed  on  the  Hoppler 
instruxient,  the  first  timer  vras  stopped  and  a  second  \ias 
started.  The  time  on  the  first  was  read  into  a  recorder, 
the  tiffxer  was  reset,  and  as  the  next  millimeter  marking 
was  passed,  restarted  as  the  other  timer  vras  stopped. 
This  process  v;as  repeated  until  the  ball  had  fallen  3  mm 
belovr  the  terminal  point  of  the  calibration  rims  reported 
in  Appendix  E. 

Individual  fall  tines  for  each  l-irm  interval  are 
provided  in  Table  K-I.   Also  tabulated  are  temperatures 
measured  betiv'een  runs  and  the  average  fall  time  for  each 
location  along  the  viscor.eter  axis. 

Inspection  of  the  data  indicates  that  the  1-mm  fall 
times  vaTv  significantly  along  the  viscometer  cylinder 
and  that  the  variations  are  roughly  reproducible  from 
one  run  to  the  next.   An  exception  is  the  first  millimeter 
of  travel  for  each  run.  Tiiies  for  this  interval  tend  to 
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be  significantly  lower  than  for  the  same  location  \ihen 
the  ball  was  started  further  up  the  axis.   Because  of 
this  disparity,  fall  times  for  "first  r.illirr.eters"  are 
not  included  in  the  computation  of  average  fall  times 
for  their  respective  locations. 

An  attempt  to  observe  time-dependent  properties  of 
the  fluid  was  made  by  correcting  the  fall  time  for  each 
l-mm  interval  by  multiplying  the  measured  fall  time  at 
each  location  by  the  average  fall  time  for  all  measurements 
divided  by  the  average  fall  time  at  the  axial  location  in 
question.  Corrected  fall  times  are  compared  in  Table  E-II 
for  the  first  10  mm.   of  travel,  whatever  the  starting  point. 
Those  runs  v;hich  terminate  before  the  ball  has  fallen  10  mm 
are  those  which  iiere   begun  at  axial  locations  near  the 
bottom  of  the  range  tested.  No  significant  dependence 
upon  time  is  apparent  after  the  first  millimeter  of  travel. 
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APPENDIX  L 
DETSRiH'INATIOil  OF  SA2'1PLE  CONCEOTRATIOIIS 

Because  of  the  variations  in  concentration  which 
occurred  during  testing,  the  solid  content  was  deter- 
mined for  each  of  the  samples  taken  for  viscometric 
analysis.  The  value  of  these  measurements  is  limited 
by  the  un-lmown  inorganic  salt  content  of  commercial 
Milling  Yello^^r  which  averages  30  percent  of  the  total 
weight  according  to  Swanson  and  Green  (I969)  and  the 
supplier. 

The  samples  used  for  concentration  m.easurenents 
averaged  ^0   ml  and  were  vreighed  out  to  the  nearest 
milligram  on  a  chain-linl:  balance.  Following  evaporation 
to  dryness  in  a  cle3.n  container  maintained  at  about  42<^  C, 
the  residue  v;as  again  vreighed,  this  time  to  -  0,1  milligram, 
a  variation  easily  v;ithin  the  calibrated  precision  of  the 
instrument.  A  typical  calculation  follov;s,  principally  to 
indicate  the  relative  magnitudes  of  the  tare,  the  solution, 
and  the  residue. 

Sample  calculation,  run  42,  — 

SamBle  &  tare:   86.535  gm      Residue  &  tare:   34.l64'^  gra 
-  Tare:   -33.353  gm      -  Tare:  -33.352o  gm 

Scjnple  mass:   53.232  gm       Residue  mass:    0,3ll7  gm 
VJeight  fraction  =0.3117/53.232  =  O.OI525 
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The  results  of  these  calculations  have  been  sunnar- 
ized  in  Table  L-I.  Two  sets  of  data  are  tabulated.   One 
set  corresponds  to  the  tests  run  in  the  rectangular 
conduit,  the  other  to  the  study  of  temperature  variability 
reported  in  Appendix  I. 

TABLE  L-I 
COriCENTRATIOITS  OF  SA.MPLES 


Sample 

Percentage 
solids 

SamDle 
104-14 

Percentage 
solids 

4-2 

1.525 

1.309 

44 

1.534 

1016-14 

1.316 

418 

1.695 

104-17 

1.586 

419 

1.716 

1011-17 

1.572 

420 

1.753 

103-20 

1.886 

420B 

1.663 

109-20 

1.922 

424 

1.719 

107-25 

2.329 

4243 

1.717 

1011-25 

2.344 

425 

1.697 

103-30 

2.939 

1012-30 

2.333 

The  dried  samples  showed  a  considerable  variation  in 
color  and  texture,  tvo  distinct  constituents predoainating: 
coarse,  deeply  tinted,  orange  crystals  and  a  fine,  nearly 
J anorphous  powder.  The  relative  proportions  of  these 
varied  widely.  For  example,  the  original  stock  solution, 
after  a  year's  sedimentation  of  fine  yellov/  powder,  yielded 
dried  residues  in  which  the  orange  crystals  heavily 
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predominated,  while  samples  from  the  fresh  stock  solution 
were  as  predominately  yellow.   Considering  the  difference 
in  weight  percentages  cited  in  Table  L-I  and  the  fact 
that  the  optical  and  rheological  properties  are  quite 
similar,  one  is  dra;m  to  the  conclusion  that  there  are 
two  conditions  present  within  Milling  Yellow  solutions: 
one,  the  deep  orange  which  characterizes  birefringent 
solutions,  is  active;  the  other,  the  light  yellow  of 
dilute  solutions,  is  inactive.  A  further  investigation 
of  the  nature  of  Milling  Yellow' s  special  properties  may 
require  the  specialized  attention  of  physical  chemists. 
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